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"When you are solving a problem, don't worry. Now, after you have 
solved the problem, then that's the time to worry." 



— Richard Feynman 



Abstract 



This dissertation is devoted to the study of thermodynamics for quantum gauge the- 
ories. The poor convergence of quantum field theory at finite temperature has been 
the main obstacle in the practical applications of thermal QCD for decades. In this 
dissertation I apply hard-thermal-loop perturbation theory, which is a gauge-invariant 
reorganization of the conventional perturbative expansion for quantum gauge theories 
to the thermodynamics of QED and Yang-Mills theory to three-loop order. For the 
Abelian case, I present a calculation of the free energy of a hot gas of electrons and 
photons by expanding in a power series in triu/T, nij/T and e^, where mo and nif are 
the photon and electron thermal masses, respectively, and e is the coupling constant. 
I demonstrate that the hard-thermal-loop perturbation reorganization improves the 
convergence of the successive approximations to the QED free energy at large coupling, 
e ~ 2. For the non- Abelian case, 1 present a calculation of the free energy of a hot gas of 
gluons by expanding in a power series in mo/T and g^, where is the gluon thermal 
mass and g is the coupling constant. I show that at three-loop order hard-thermal-loop 
perturbation theory is compatible with lattice results for the pressure, energy density, 
and entropy down to temperatures T ~ 2 — 3 Tc. The results suggest that HTLpt 
provides a systematic framework that can be used to calculate static and dynamic 
quantities for temperatures relevant at LHC. 
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Chapter 1 

Introduction 



The beginning of experiments at the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL) in 1999 marked the beginning of a new era in ultrarelativis- 
tic heavy-ion collisions. One of the primary goals of the RHIC program is to discover 
and study the quark-gluon plasma (QGP) whose existence is predicted by quantum 
chromodynamics (QCD). In addition, looking forward, ultrarelativistic heavy-ion colli- 
sion experiments are part of the Large Hadron Collider (LHC) program at European 
Organization for Nuclear Research (CERN). The LHC experiments, for which full beam 
runs are scheduled in 2011, will provide data on heavy ion collisions at center of mass 
energies of 5.5 TeV per nucleon pair collision and will open a new chapter in the study 
of partonic matter under extreme conditions. 

For the RHIC and LHC experiments to have the greatest possible impact on science, 
it is essential to make as close a connection to the fundamental theory of QCD as 
possible. There is an urgent need for theoretical analysis that is based rigorously 
on QCD but which can also make contact with more phenomenological approaches, 
particularly in the area of equilibrium and non-equilibrium dynamics of QCD at 
intermediate coupling, ^ ~ 2, or equivalently dg = / (47r) ~ 0.3. 

We have to be extremely careful when dealing with this intermediately coupled 
region. Naively, ~ 2 seems to suggest the breakdown of perturbation theory in this 
region. This is also in line with the observations from the early RHIC data that the 
state of matter created there behaved more like a strongly coupled fluid than a weakly 
coupled plasma (ll. As a result, the term "quark-gluon plasma" might need to be 
modified to "quark-gluon liquid", and a description in terms of hydrodynamics or 
AdS/CFT correspondence might be more appropriate. However on the other hand, 
^ ~ 2 is not huge especially when considering that as = / (47r) ~ 0.3 is still a 
small number. So people have not yet totally lost faith in perturbation theory and 
as a payback observables like jet quenching (21 and elliptic flow O have been able 
to be described using a perturbative formalism. Therefore it seems that a complete 
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understanding of QGP would require knowledge from both strong-coupling and weak- 
coupling formalisms, and in this dissertation I focus on the latter approach. 

Thermodynamics describes the bulk properties of matter in or near equilibrium 
which are theoretically clean and well defined. The calculation of thermodynamic 
functions for finite temperature field theories has a long history. In the early 1990s the 
free energy was calculated to order for massless scalar (p'^ theory (H 13, quantum 
electrodynamics (QED) ©IHl and QCD (51, respectively. The corresponding calcula- 
tions to order were obtained soon afterwards ll7HT2ll. Recent results have extended 
the calculation of the QCD free energy by determining the coefficient of the g log g 
contribution I1T3| . For massless scalar theories the perturbative free energy is now 
known to order g^ im and g^ log g llTSl . 

Unfortunately, for all the above-mentioned theories the resulting weak-coupling 
approximations, truncated order-by-order in the coupling constant, are poorly con- 
vergent unless the coupling constant is tiny. Therefore a straightforward perturbative 
expansion in powers of Us for QCD does not seem to be of any quantitative use even 
at temperatures many orders of magnitude higher than those achievable in heavy-ion 
collisions. 

The poor convergence of finite-temperature perturbative expansions of thermo- 
dynamic functions stems from the fact that at high temperature the classical solution 
is not described by massless gluonic states. Instead one must include plasma effects 
such as the screening of electric fields and Landau damping via a self-consistent hard- 
thermal-loop (HTL) resummation IIT6ll . The inclusion of plasma effects can be achieved 
by reorganizing perturbation theory. 

There are several ways of systematically reorganizing the finite-temperature per- 
turbative expansion IIT7l - IT9ll . In this dissertation I will focus on the hard- thermal-loop 
perturbation theory (HTLpt) method Il20ti24l . The HTLpt method is inspired by varia- 
tional perturbation theory ||25] - |27| . HTLpt is a gauge-invariant extension of screened 
perturbation theory (SPT) ||28] - |32| , which is a perturbative reorganization for finite- 
temperature massless scalar field theory. In the SPT approach, one introduces a single 
variational parameter which has a simple interpretation as a thermal mass. In SPT 
a mass term is added to and subtracted from the scalar Lagrangian, with the added 
piece kept as part of the free Lagrangian and the subtracted piece associated with the 
interactions. The mass parameter is then required to satisfy a variational equation 
which is obtained by a principle of minimal sensitivity. This naturally led to the idea 
that one could apply a similar technique to gauge theories by adding and subtracting 
a mass in the Lagrangian. However, in gauge theories, one cannot simply add and 
subtract a local mass term since this would violate gauge invariance. Instead, one adds 
and subtracts an HTL improvement term which modifies the propagators and vertices 
self-consistently so that the reorganization is manifestly gauge invariant ||33) . 
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This dissertation focuses on the study of thermodynamics for gauge theories. In 
the rest of this chapter, a brief introduction to statistical physics and thermal QCD is 
provided. In Chapter |2| we show the emergence of infrared divergences in thermal 
field theory and how the weak-coupling expansion treats them systematically. HTLpt 
is introduced in Chapter |3j where we discuss its formalism as well as techniques that 
make HTLpt calculations tractable. Chapters |4] and |5] are devoted to the study of 
thermodynamics to three-loop order using HTLpt for QED and Yang-Mills theory, 
respectively. We summarize in Chapter |6] together with a brief outlook for the real-time 
application of HTLpt. 



1.1 Statistical physics and quantum partition function 

For a relativistic system which can freely exchange energy and particles with its sur- 
roundings, the most important function in thermodynamics is the grand canonical 
partition function 

Z = ^ e-f,/T _ ^ = Tre-«/^ (1.1.1) 

states states 

Here £, is the energy of the state \Si) and "H is the Hamiltonian of the system. The 
temperature of the system is denoted by T, and since this dissertation only concerns 
with high temperature physics, the chemical potential of the particles in the system 
is set to zero for simplicity. All of the thermodynamic properties can be determined 
from ( 1.1. 1| . For example, the pressure, entropy and energy are given by 



V = "^^^f^ , (1.1.2) 

S=^M. (1.1.3) 
g = -pv + TS, (1.1.4) 

where V is the volume of the system. Typically, the width L of a system is much larger 
than the inverse temperature, (i.e. L ^ In/ T), such that one can use the infinite 
volume limit to describe the thermodynamics of a finite volume to good approximation. 
The advantage of the infinite volume limit is that field theoretic calculations simplify. 
In all calculations performed in this thesis, this infinite volume limit is taken. Then it 
turns out that log Z becomes proportional to V, such that the pressure becomes 

P = . (1.1.5) 

The extension to field theory is straightforward. If is the Hamiltonian of a 
quantum field theory in d-dimensional space and hence {d + 1) -dimensional spacetime. 
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then the partition function ( 1.1.1 1 is 

Z = Tre-^/^ = J P<pe-/o^'''^^/'^''-^^W , (1.1.6) 

with C the Lagrangian density of the theory and periodic boundary conditions 

<p(0,x) = <p(l/r,x) . (1.1.7) 

for bosonic fields (p. For fermionic fields, it turns out that to implement Pauli statistics 
one must impose anti-periodic boundary conditions 

<p(ax) = -<p(l/T,x) . (1.1.8) 



1.2 QCD at finite temperature 

Quantum chromodynamics is a gauge theory for the strong interaction describing the 
interactions between quarks and gluons. The QCD Lagrangian density in Minkowski 
space reads 

^QCD = -^Tr [G^.yG^''] +J^ipi [i^W^ - rrn] xPi + C^t + £ghost ■ (1-2.1) 

The gluon field strength is G^y = d^Ai, — dyA^ — ig[A^, Ay\. The gluon field is A^ = 
A'^t", with generators of the fundamental representation of SU(3) normalized so that 
Tr t"t^ = S"^ /2. In the quark sector there is an explicit sum over the Nf quark flavors 
with masses m, and D^^ = 3^, — igA^ is the covariant derivative in the fundamental 



representation. The Lagrangian ( 1.2.1 1 is mathematically simple and beautiful, however 
in order to carry out a physical calculation with it, a gauge fixing is needed to remove 
unphysical degrees of freedom. The ghost term >Cghost depends on the choice of the 
gauge-fixing term Cgf. One popular choice for the gauge-fixing term that depends on 
an arbitrary gauge parameter ^ is the general covariant gauge: 

1, 



{H'Ai.Yl . (1.2.2) 
The corresponding ghost term in the general covariant gauge reads 



Vost - -rdy+gf^'fdt'iAli^'), (1.2.3) 

where t] and f] are anti-commuting ghosts and anti-ghosts respectively and /"'"^ is 
structure constant of SU(3). 

The finite temperature QCD partition function is obtained by a Wick rotation of 
the theory from Minkowski space to Euclidean space. It is achieved by the substitution 
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t = ZT with t being the Minkowski time and t being the Euclidean one. The resulting 
Euclidean partition function is 



/ 



VAfiVxpVtpVi]Vt] exp 



l/T 

dr / d^xCn 







QCD 



(1.2.4) 



with Cqqjj the Wick-rotated Lagrangian denity. Feynman rules are exactly the same 
as in zero-temperature field theory except that the imaginary time T is now compact 
with extent l/T. To go from t to frequency space, we should perform a Fourier 
series decomposition rather than a Fourier transform. The only difference with zero- 
temperature Feynman rules will then be that loop frequency integrals are replaced by 
loop frequency sums: 

with the sum over discrete imaginary-time frequencies known as Matsubara frequencies 

cOn = 2nnT bosons, (1.2.6) 

co„ = {2n + l)nT fermions. (1.2.7) 

to implement the periodic or anti-periodic boundary conditions. A detailed explanation 



of the imaginary-time formalism is given in Appendix A.12 



1.3 Beta function and asymptotic freedom 

The beta function of a quantum field theory encodes the dependence of a coupling 
parameter g on the energy scale /y of a given physical process. It is defined by the 
relation: 

= f I . (1.3.1) 

This dependence on the energy scale is known as the running of the coupling parameter, 
and theory of this kind of scale-dependence in quantum field theory is described by 
the renormalization group which refers to a mathematical apparatus that allows one to 
investigate the changes of a physical system as one views it at different distance scales. 

To lowest order in the coupling constant a beta function is either positive indicating 
the growth of charge at short distance or negative indicating the decrease of charge at 
short distance. Until 1973, only examples of the former were known *. The discovery 
that only non-Abelian gauge theories allow for a negative beta function is usually 

* 't Hooft reported a similar discovery at the Marseille conference on renormalization of Yang-Mills 
fields and applications to particle physics in 1972 without publishing it l34l . 
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credited to Gross and Wilczek ||35l, and to Politzer ||36l. The solution to ( |1.3.1| for QCD 
reads 

ccsiH) = ^ = - , (1.3.2) 

471 (ll-?N^)log(?//AQCD) ^ ' 

which clearly shows asymptotic freedom Il35ll36l , i.e. a — > as — > oo. The parameter 
-^QCD is a scale above which the theory works "chosen" by the world in which we 
live. It is well known that QCD exhibits confinement at large distances or low energies 
which terminates the validity of perturbation theory due to the infrared growth of the 
coupling. However it is precisely the asymptotic freedom that ensures the possibility 
of a perturbative treatment for the ultraviolet sector of the theory which sets the stage 
to study the high temperature phase of non-Abelian theory in this dissertation. 
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In this chapter, and in the rest of the dissertation, we consider thermal field theories at 
high temperatures, which means temperatures much higher than all zero-temperature 
masses or any mass scales generated at zero temperature. 

It has been known for many years that naive perturbation theory, or the loop 
expansion breaks down at high temperature due to infrared divergences. Diagrams 
which are nominally of higher order in the coupling constant contribute to leading 
order in g. A consistent perturbative expansion requires the resummation of an infinite 
subset of diagrams from all orders of perturbation theory. We discuss these issues next. 



2.1 Scalar field theory 

We start our discussion by considering the simplest interacting thermal field theory, 
namely that of a single massless scalar field with a (p^ interaction. The Euclidean 
Lagrangian is 

^ = 1(3^^)' + ^^'- (2-1-1) 

Perturbative calculations at zero temperature proceed by dividing the Lagrangian into 
a free part and an iateractiag part according to 

Aree = , (2.1.2) 

Ant = ^(p'' . (2.1.3) 

Radiative corrections are then calculated in a loop expansion which is equivalent to a 
power series in g^. We shall see that the perturbative expansion breaks down at finite 
temperature and the weak-coupling expansion becomes an expansion in g rather than 
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We will first calculate the self-energy by evaluating the relevant diagrams. The 
Feynman diagrams that contribute to the self-energy up to two loops are shown in 
Fig.lZT] 



1-loop : 2-loop : 



Figure 2.1: One- and two-loop scalar self-energy graphs. 

The one-loop diagram is independent of the external momentum and the resulting 
integral expression is 



n(i) = -a^i — 

2^fpP2' 



- *_t2 
24 ' 

= , (2.1.4) 

where the superscript indicates the number of loops. The notation P = {Po,p) rep- 
resents the Euclidean four-momentum. The Euclidean energy Pq has discrete values: 
Pq = InnT for bosons and Pq = (2w -|- l)7rT for fermions, where n is an integer. 
Eq. ( 2.1.4[ | represents the leading order thermal mass of our scalar field. The sum- 



integral ^p, which is defined in Eq. (B.0.1 1, represents a summation over Matsubara 



frequencies and integration of spatial momenta md = 3 — 2e dimensions * . The above 
sum-integral has ultraviolet power divergences that are set to zero in dimensional 
regularization. We are then left with the finite result ( 2.1.4[ |, which shows that ther- 



mal fluctuations generate a mass for the scalar field of order gT. This thermal mass 
is analogous to the Debye mass which is well-known from the nonrelativistic QED 
plasma. 

We next focus on the two-loop diagrams and first consider the double-bubble in 
Fig.|22](b). 



(a) (b) (c) 

Figure 2.2: Bubble diagrams contributing to the scalar self -energy. 



This diagram is also independent of the external momentum and gives the following 
*For an introduction to thermal field theory and the imaginary time formalism see Refs. l37l and l38l . 
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sum-integral 

n(2^) = -l^'i 44. (2.1.5) 



4^ TpqP^Q^ ■ 

This integral is infrared divergent. The problem stems from the middle loop with two 
propagators. In order to isolate the source of the divergence, we look at the contribution 
from the zeroth Matsubara mode to the Q integration 



with defined in Eq. (C.0.1 The integral J^l/q^ behaves like 1/ q, as a result Eq. (2.1.6 
is linearly infrared divergent as q ^ 0. This infrared divergence indicates that naive 
perturbation theory breaks down at finite temperature. However, in practice this 
infrared divergence is screened by a thermally generated mass and we must somehow 
take this into account. The thermal mass can be incorporated by using an effective 
propagator: 

MP) = , (2.1.7) 

with m ^ gT <^T. 

If the momenta of the propagator is of order T or hard, clearly the thermal mass is 
a perturbation and can be omitted. However, if the momenta of the propagator is of 
order ^T or soft, the thermal mass is as large as the bare inverse propagator and carmot 
be omitted. The mass term in the propagator ( 2.1. 7| provides an infrared cutoff of order 



gT. The contribution from ( 2.1.6| would then be 

Since m ~ gT , this shows that the double-bubble contributes at order ^'^T^ to the 
self -energy and not at order g^T"^ as one might have expected. Similarly, one can show 



that the diagrams with any number of bubbles like Fig. 2.2 : are all of order ^ . Clearly, 



naive perturbation theory breaks down since the order-^ correction to the thermal 
mass receives contributions from all loop orders. On the other hand, the three-loop 



diagram shown in Fig. 2.3 is of order ^*T^ and thus subleading. Therefore, we only 
need to resum a subset of all possible Feynman graphs in order to obtain a consistent 
expansion in g. 

If we use the effective propagator to recalculate the one-loop self -energy, we obtain 
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Figure 2.3: Subleading three-loop self -energy diagram. 



r 



24 



p p2 + 
An 



(2.1.9) 



where here, and in the following, the prime on the sum-integral indicates that we have 
excluded the n = mode from the sum over the Matsubara frequencies. The order 
corresponds to the summation of the bubble diagrams in Fig. 



2.2 



which can be 



verified by expanding the effective propagator ( 2.1.7[ | around m = 0. Thus by taking 
the thermal mass into account, one is resumming an infinite set of diagrams from all 
orders of perturbation theory. 

The self-energy ( 2.1.4[ | is the first example of a hard thermal loop (HTL). Hard 
thermal loops are loop corrections which are g^T'^/ times the corresponding tree- 
level amplitude, where P is a momentum that characterizes the external lines. From 
this definition, we see that, whenever P is hard, the loop correction is suppressed by 
and is thus a perturbative correction. However, for soft P, the hard thermal loop 
is 0{\) and is therefore as important as the tree-level contribution to the amplitude. 
These loop corrections are called "hard" because the relevant integrals are dominated 
by momenta of order T. Also note that the hard thermal loop in the two-point function 
is finite since it is exclusively due to thermal fluctuations. Quantum fluctuations do not 
enter. Both properties are shared by all hard thermal loops. 

What about higher-order n-point functions in scalar thermal field theory? One can 
show that within scalar theory the one-loop correction to the four-point function for 
high temperature behaves as Il39ll 



(2.1.10) 



oc ^4iog(T/p) , 



where p is the external momentum. Thus the loop correction to the four-point function 
increases logarithmically with temperature. It is therefore always down by g^ log(l/^), 
and it suffices to use a bare vertex. More generally, it can be shown that the only hard 



thermal loop in scalar field theory is the tadpole diagram in Fig. 2.1 and resummation 



is taken care of by including the thermal mass in the propagator. In gauge theories, the 
situation is much more complicated as we shall see in the next section. 
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2.2 Gauge theories 

In the previous section, we demonstrated the need for resummation in a hot scalar 
theory. For scalar theories, resummation simply amounts to including the thermal 
mass in the propagator and since the running coupling depends logarithmically on the 
temperature, corrections to the bare vertex are always down by powers of ^ log 1 /^. In 
gauge theories, the situation is more complicated. The equivalent HTL self -energies are 
no longer local, but depend in a nontrivial way on the external momentum. In addition, 
it is also necessary to use effective vertices that also depend on the external momentum. 
It turns out that all hard thermal loops are gauge-fixing independent 111611401 - 1441 . This 
was shown explicitly in covariant gauges. Coulomb gauges, and axial gauges. They 
also satisfy tree-level like Ward identities. Furthermore, there exists a gauge invariant 
effective Lagrangian, found independently by Braaten and Pisarski II33I and by Taylor 
and Wong Il43l , that generates all of the hard thermal loop n-point functions. From a 
renormalization group point of view this is an effective Lagrangian for the soft scale gT 
that is obtained by integrating out the hard scale T. We return to the HTL Lagrangian 
in Chapter |3] 



2.2,1 Polarization tensor 



We next discuss in some detail the hard thermal loop for the vacuum polarization tensor 
YW^ . For simplicity, we focus on QED here. The Feynman diagram for the one-loop 
photon self-energy is shown in Fig. 2.4 and results in the following sum-integral 

"K7''(K-f)7 



Tr 



K2(]<C-P)2 



(2.2.1) 



where Tr denotes the trace over Dirac indices. After taking the trace, the self-energy 
becomes 



2^ w¥y _ _^^,.^2^ 



2 j; K^'P^ + K^p^ 



{K}K^{K-P) 



4e- 



{K} K^K-P)^ ' 



(2.2.2) 



where we have assumed, for now, that d = 3. Since we are interested in the high- 
temperature limit, we may assume that K ^ P because the leading contribution in T to 
the loop integral is given by the region K ~ T. With this assumption, the self-energy 
simplifies to 



- 4^f"'e^/ — 



(2.2.3) 



We first consider the spatial components of n''^'(P). The sum over Matsubara frequen- 
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Figure 2.4: One-loop photon self-energy diagram. 



cies can be evaluated using 

TV ' 



l-npik) -nf(|p-k| 



Ak\p-k\ 

- nf{k) +nF(|p - k|) 



iPo - k-\p-k\ ~'~ fPo + ^ + |p - k| 
1 



+ 



_/Po-'^+ |p-k| iPo + k-\p-k 



, (2.2.4) 



which is derived from a contour integral in the complex energy plane. The second term 
in Eq. ( 2.2.3| l is rather simple. We obtain 



n'/{P) = -2A.-/J(i-2„,W) + 2.^X^ 



X i^(l-nf{k) -nf{\k-p 
- (-nF(/c) + nH|k-p|) 



zPo-^-|k-p| iPo+k+\k-p 
-1 1 " 



fPo-^+|k-p| iPo + k-\k-p\ 



(2.2.5) 



where nf{x) = 1/ (exp(^x) + 1) is the Fermi-Dirac distribution function. The zero- 
temperature part of Eq. ( 2.2.5| is logarithmically divergent in the ultraviolet. This term 
depends on the external momentum and is cancelled by standard zero-temperature 
wavefunction renormalization. We next consider the terms that depend on temperature. 
In the case that the loop momentum is soft, the Fermi-Dirac distribution functions can 
be approximated by a constant. The contribution from the integral over the magnitude 
of k is then of order and subleading. When the loop momentum is hard, one can 
expand the terms in the integrand in powers of the external momentum. We can then 
make the following approximations 



nK|k-p|) ^ nrik)-^^p-k, 



|k — p| k — p k , 



(2.2.6) 
(2.2.7) 
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where k = k/A; is a unit vector. Thus the angular integration decouples from the 
integral over the magnitude k. This implies 

^ ^ Jo dk J An -fPg + p-k 



3 J 47r-zPo + p-k 

The other components of the self-energy tensor H'"'(P) are derived in the same 
manner or obtained using the transversality of polarization tensor: 

pFuf''{P) = 0. (2.2.9) 

One finds EH 

n«(P) ^ ill (f€l^f^^A, a.2.10) 



3 47r-zTo + p-k 

hO^(P) = g^/^ Al. (2.2.11) 

In d dimensions, we can compactly write the self-energy tensor as 

nf"'(P) = mf,[rf'^(P,-P)-N^N^] , (2.2.12) 

where N specifies the thermal rest frame is canonically given by N = (— f, 0). We have 
defined 

n^l = -m-De^f^^^^, = (2.2.13) 
and the tensor T^^(P, Q), which is defined only for momenta that satisfy P + Q = 0, is 

V'^^P,-P) = (vn'^y . (2.2.14) 

The angular brackets indicate averaging over the spatial directions of the light-like 
vector Y = (— f,y). The tensor T'^" is symmetric in ji and v. Because of transverality 
and the rotational symmetry around the p-axis, one can express the self -energy in terms 
of two independent functions, H7(P) and Hl(P): 

Hf-(P) = H,(P) ^ " ^ 

nr(P) + ^HL(P) 

13 
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where the functions YIt{P) and YIl{P) are 

1 



Hr(P) 



(3'' - ff)n^P) , 



noo(p) . 



(2.2.16) 
(2.2.17) 



In three dimensions, the self-energies n7(P) and rLiiP) reduce to 



UriP) 
ndP) 



m 



P2 

_D 10 

2 p2 



P^ , iPo + P 
log? 



m7 



2fPop ° iPo - V 

fPo , iPo + p' 

— loe 

2p ^iPo-p 



(2.2.18) 
(2.2.19) 



The hard thermal loop in the photon propagator was first calculated by Silin more than 
forty years ago 1451 . The hard thermal loop in the gluon self -energy was first calculated 
by Klimov and Wei don ||46ll47| . It has the same form as in QED, but where the Debye 
mass mu is replaced by 



ml 



(2.2.20) 



where Ca = Nc is the number of colors and Nf is the number of flavors. When d = 3 
the QCD gluon Debye mass becomes 



(2.2.21) 



2.2.2 Fermionic self-energy 

The electron self-energy is given by 



E(P) = mjJiP) , 



where 



V\P) 



and mj\s the thermal electron mass 



P^ I ~ ' 



= — 3e^ 



In QCD, the quark mass is given by 



(2.2.22) 



(2.2.23) 



(2.2.24) 



(2.2.25) 
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2,2.3 Higher n-point functions 

In gauge theories, there are also hard thermal loops involving vertices. For instance, 
the one-loop correction to the three-point function in QED, can compactly be written as 

Tf" (P, Q, R) = 7^' - mjf (P, Q, R) , (2.2.26) 

where the tensor in the HTL correction term is only defined for P 

The quark-gluon vertex satisfies the Ward identity 

P''r''(P,Q,R) = S-\Q) - S-^R) , (2.2.28) 

where S{q) is the resummed effective fermion propagator. 

In QED there are, in fact, infinitely many amplitudes with hard thermal loops. 
To be precise, there are hard thermal loops in all n-point functions with two fermion 
lines and n — 2 photon lines. In non-Abelian gauge theories such as QCD, there are in 
addition hard thermal loops in amplitudes with n gluon lines llT6l . 



-Q + R = 0: 

(2.2.27) 



2.3 Weak-coupling expansion 

The Braaten-Pisarski resummation program has been used to calculate the thermody- 
namic functions as a weak-coupling expansion in g. They have now been calculated 
explicitly through order log ^ for massless (p^ theory ||4l|5l|7l|8l[l4l[l5), through order 
for QED IHllEISHni, and through order g^ \ogg for QCD ElTlHEl. In this section, 
we review these calculations in some detail. 



2.3.1 Scalar field theory 

The simplest way of dealing with the infrared divergences in scalar field theory is to 
reorganize perturbation theory in such a way that it incorporates the effects of the 
thermally generated mass m into the free part of the Lagrangian. One possibility is to 
divide the Lagrangian ( 2.1.1| into free and interacting parts according to 



Aree = ^ + V , (2.3.1) 

Cin, = ^4>'-lmY- (2.3.2) 
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-^0 -^la -^Ib ^2a ^2b ^Ic ^2d 

Figure 2.5: Diagrams which contribute up to three-loop order in scalar perturbation theory. 
A boldfaced x indicates an insertion of m^. 



Both terms in Eq. (2.3.2 1 are treated as interaction terms of the same order, namely g^. 
However, the resummation implied by the above is rather cumbersome when it comes 
to calculating Green's function with zero external energy. Static Green's functions 
can always be calculated directly in imaginary time without having to analytically 
continue them back to real time. This implies that we can use Euclidean propagators 
with discrete energies when analyzing infrared divergences which greatly simplifies the 
treatment. In particular, since only propagators with zero Matsubara frequency have no 
infrared cutoff of order T, only for these modes is the thermal mass of order gT relevant 
as an IR cutoff. Thus, another possibility is to add and subtract a mass term only for the 
zero-frequency mode. This approach has the advantage that we do not need to expand 
the sum-integrals in powers of /T'^ in order to obtain the contribution from a given 
term in powers of g^. We will follow this path in the remainder of this section and write 

Aree = \{Mf + \^^'P^^P^,0 r (2-3.3) 

^int = ^^(p" - ImY^Pofl ■ (2.3.4) 
The free propagator then takes the form 

A(P) = + . (2.3.5) 

This way of resumming is referred to as static resummation ||48| . It is important to point 
out that this simplified resummation scheme can only be used to calculate static quan- 
tities such as the pressure or screening masses. Calculation of dynamical quantities 
requires the full Braaten-Pisarski resummation program. The problem is that the calcu- 
lation of correlation functions with zero external frequencies cannot unambiguously be 
analytically continued to real time Il49ll . 

Perturbative expansion 

Following the decomposition in (2.3.3 1 and (2.3.4 1, the partition function for (p^ theory 



reads 

P^g-CSfree+Sint) ^ (2.3.6) 
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with the free actior\ Sf^ee and the action due to interactions Smt defined by 

nl/T 



Sfree ~ / ^-T d X /Ifree r 

Jo J 

dz J d^X dnt . 



Sint — 



(2.3.7) 
(2.3.8) 



Expanding in powers of Sint, ( 2.3.6| becomes 



(-Sint)" 



«=0 



nl 



(2.3.9) 



Taking logarithm for both sides, we get 



logZ = log 



V(pe 



log 



= logZo + logZi 



(2.3.10) 



In this way, the partition function has been divided into two parts. The free part log Zq 
describes the classical limit of the theory, i.e. a gas of non-interacting scalar particles, 
which can be evaluated analytically due to the fact that it is a quadratic or Gaussian 
function in the field (p. The interacting part log Zj contains quantum corrections to the 
classical theory, which is accessed perturbatively as a power series in Sint- Using the 
relation log det A = Tr log A it is easy to show that 



logZo = O = -^Trlogl-S^ 



(2.3.11) 



where the single closed loop denotes the corresponding Feynman diagram contributing 
to log Zq. As an example of the perturbative expansion of static resummation, we next 
consider the calculation of the free energy which is the negative pressure T = —V, 
through order in scalar field theory. This involves the evaluation of vacuum graphs 
up to three-loop order shown in Fig. |2.5[ 



One loop 

The one-loop contribution to the free energy is 

-^0 = ^T|^log(p2 + m2) + l^^ogp2 
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Using the integrals and sum-integrals contained in the appendices the result for this 
diagram in the limit e — >^ is 



90 



j,4 



Tnr" 

1271 



Two loops 

The two-loop contribution to the free-energy is given by 



with 



T\ — J-"ia + ^ib / 

1 2^ / 1 

--m^T / y . 

2 JvP + m 



The result for these diagrams in the limit e — >^ is 



-^la — 



90 1^4 \An 

_ 5\/6 /^\3 
2 Utt/ 



1 + e (^4 + 4 
1 + e ( 4 + 2 



a-i) 



V47rTy 



4e 



(2.3.13) 



(2.3.14) 



(2.3.15) 
(2.3.16) 



2e 



AtzTJ \2mJ 



J^lb — 



~90 1~ Vin) ' 



(2.3.17) 
(2.3.18) 



where we have kept all terms that contribute through order e, because they are needed 
for the counterterm diagrams in the three-loop free energy. 



Three loops 

The three-loop contribution is given by 



-^2 = -F2a + J^2b + ^2c + J^2d + ^^ig" , (2.3.19) 



where the expressions for the diagrams are 



= -^^'(^/p^^+^pf) (^/p(p2 + „2)2+^, 



-F2b = -^^gf^ 



PQrP^ Q2 K2 (P + Q + R)2 
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(2.3.21) 



■ (p + q + r)2 + ' 
^2d = I . (2.3.23) 



p (p2 + OT2)2 



Note that the basketball diagram J^2b is the only diagram through three loops that 
cannot be decomposed to solely a set of one loop integrals, which therefore needs 
special treatment II50I . The result for these diagrams in the limit e — > is 



J"2a = 



n^T^ j 5^6 / g x3 5 
90 I 8 UttJ 8 V47r 



+ 27-8 + 4 



5^6 /^x5 
4 Utt 



1 C'(-i) 



4e 



90 \ 4 Utt/ 

. ^r^^'Ti 

2 Utt. 



2b 



1 1 o r(-i) , r(-3) , 91 
e^^a-i) a-3)^i5 





( 


a-i)J 






\2mJ 




( ^ ^ 




V47rT. 





6€ 



(2.3.24) 



8-41og2 



\2m) 



J-2c 



•^2d — 



90 



5v/6/gx3 15/gx4 
4 Utt/ 2 Utt 



"90 8^ \An 



(2.3.25) 
(2.3.26) 
(2.3.27) 



Free energy through 

Combining the one-, two-, and three-loop contributions given by Eqs. ( |2.3.12| , (2.3.14 1, 
and (2.3.19 1, respectively, gives the free energy through order 



T_ 
-^0 



4 V47r/ 

2 r(-3) 
3 a-3) 



5v/6 



3 V47r 



3 15 



log 



47rT 



59 1 

h -7 ■ 

45 3 ^ 



4 r(-l) 
3^-1) 



47r 



_ 15^6 
2 

2r(-i) 



log 



2 2 1 

3 log - + -7 3^^_-^ 



V47r/ 



47rr 
5 



V47r/ 6 



(2.3.28) 



where J-q = — 7r2T^/90 is the free energy of an ideal gas of non-interacting scalar 
bosons. The pressure through order ^ was first calculated using resummation by 
Parwani and Singh [Zl and later by Braaten and Nieto using effective field theory fHl. 
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With a = [g/ An)^, the renormalization group equation for the coupHng reads 

da 



d}i 



(2.3.29) 



Using Eq. ( 2.3.29[ l^ one can verify that the free energy ( 2.3.28[ l is RG-invariant up to 
corrections of order g^ log g. 




Figure 2.6: Weak-coupling expansion for the free energy of massless theory normalized 
to that of an ideal gas as a function of g{2nT) to orders g^, g^, g"^, and g^ are shown as 
bands that correspond to varying the renormalization scale }i by a factor of 2 around 2nT. 
This figure is adapted from Ref. Ill9ll . 



In Fig. 2.6 we show the successive perturbative approximations to the free energy 
as a function of g{2nT). The bands are obtained by varying the renormalization scale 
^ by a factor of 2 around the central value }i = 27rT. The lack of convergence of the 
weak-coupling expansion is evident from this figure. The band obtained by varying 
ji hy a factor of 2 is not necessarily a good measure of the error, but it is certainly a 
lower bound on the theoretical error. Another indicator of the theoretical error is the 



deviation between successive approximations. We can infer from Fig. 2.6 that the error 
grows rapidly when ^(27rT) exceeds 1.5. 



2.3.2 Gauge theories 

hi this section, we discuss the application of weak-coupling expansion to gauge theo- 
ries. The EucUdean Lagrangian for an SU(Nc) gauge theory with Nf fermions in the 
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fundamental representation is 

c = ^g;,g;, + ipj,,d^ip, (2.3.30) 

where G^y = d^A^ — dyA^^ + gf'^'^A'j^A^ is the field strength, g is the gauge coupling 
and /"'"^ are the structure constants. The covariant derivative is = — igA^^T", 
where T" are the generators in the fundamental representation. 

The constants and Ca are the dimension and quadratic Casimir invariant of the 
adjoint representation, with 

S"" = dA , fabcjrdbc ^ (-^^ad _ (2.3.31) 

c?f is the dimension of the total fermion representation , and Sf and S2F are defined in 
terms of the generators T" for the total fermion representation as 

Sf = ;^tr(r2) , S2F = ;^tr[(T2)2] , (2.3.32) 

O-A UA 

where = TT". For SU(Nc) with Nf fermions in the fundamental representation, the 
standard normalization of the coupling gives 

1 _ 1 

dA=N^-l, Ca=Nc, df = NcNf, SF = -Nf, S2F = ^^N;. 

(2.3.33) 

For U(l) theory, relabel ^ as e and let the charges of the Nf fermions be qte. Then 

dA = l, Ca=0, dF = Nf, SF = Y2^h S2F = Y2qf. (2.3.34) 

If we are only interested in static quantities, we can apply the same simplified 
resummation scheme also to gauge theories Q [TT]| . Thus we are interested in the 
static limit of the polarization tensor H^,i,. In that limit H^ vanishes and TLi = rn^. 
In analogy with the scalar field theory, we rewrite the Lagrangian by adding and 
subtracting a mass term \m^AQAQ5p^^fl. One of the mass terms is then absorbed into 
the propagator for the timelike component of the gauge field Aq, while the other is 
treated as a perturbation. 

The free energy through g^ requires the evaluation of diagrams up to three loops. 
The strategy is the same in the scalar case, where one distinguishes between hard and 
soft loop momenta. The result reads 
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22 

y 
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38rM) _148rM) _ 64 
AnT^ 3 a-3) 3 ^(-1) ^ 5 



C.S J? log ^ 
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ir(-3) 74r(-i; 
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87 



47 
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60 
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105 



+ 24 log 2 



(471) 



Ca + Sf 



Q 176 log 



47rT 



1767 - 247r^ - 494 + 264 log 2 



CaSf ( 112 log ^ + 1127 + 72 - 128 log 2 



S^(-641og 



144S2F 




647 + 32 -128 log 2 



(2.3.35) 



The one-loop beta function for an SU(Nc) gauge theory with Nf fermions reads 
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In ' 



(2.3.36) 



written in terms of as = g /An. Using Eq. (2.3.36 \ one can verify that the free energy 



(2.3.35 I is RG-invariant up to corrections of order g° log^. 

The free energy for QCD through order was first derived by Arnold and Zhai IS]. 
Later it was extended to order by Zhai and Kastening lUTl using the above resum- 
mation techniques, and by Braaten and Nieto using effective field theory IIT2II . The 
order-^^ contribution is the last contribution that can be calculated using perturbation 
theory. At order g^, although the electric g^ log g contribution is still perturbatively 
accessible ||T3ll , however perturbation theory breaks down due to infrared divergences 
in the magnetic sector ||5n[52l . 

the free energy of QCD {Nc = 3) with Nj = 2is shown as a function 



In Fig. 



2.7 



of the temperature T/Tc, where Tc is the critical temperature for the deconfinement 
transition. In the plot we have scaled the free energy by the free energy of an ideal gas 
of quarks and gluons which for arbitrary Nc and Nj is 



ideal 



n^ 
45 



:r N. 



NcNf 



(2.3.37) 
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Figure 2.7: Weak-coupling expansion for the free energy of SU(3) gauge theory with Nf = 2 
normalized to that of an ideal gas as a function of T/Tc to orders g^, g'^, g*, and g^ are 
shown as bands that correspond to varying the renormalization scale }i by a factor of 
2 around InT ||2TI . Also shown is the lattice estimate by Karsch et al. Il53ll for the free 
energy. The band indicates the estimated systematic error of their result which is reported 
as (15 ± 5)%. This figure is adapted from Ref. ll2TI . 

The weak-coupling expansions through orders g^, g^, g'^, and g^ are shown as bands 
that correspond to varying the renormalization scale, }i, by a factor of 2 around the 
central value }i = InT. As successive terms in the weak-coupling expansion are added, 
the predictions change wildly and the sensitivity to the renormalization scale grows. 
It is clear that a reorganization of the perturbation series is essential if perturbative 
calculations are to be of any quantitative use at temperatures accessible in heavy-ion 
collisions. 
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Hard-Thermal-Loop Perturbation 
Theory 

In this chapter, we introduce the hard-thermal-loop perturbation theory, which is 
a gauge-invariant resummation scheme for thermal gauge theories. We discuss its 
formalism and properties, as well as some technicalities which make the evaluation of 
loop diagrams tractable analytically. This chapter forms the basis for the discussions in 
the rest of the dissertation. 



3.1 Introduction 

One possible conclusion from the bad convergence of the weak-coupling expansion 



of QCD free energy in Section 2.3 is that the quark-gluon plasma is completely non- 
perturbative, and that it can only be studied by nonperturbative methods like lattice 
gauge theory. This would be a very unfortunate conclusion from the perspective of the 
search for the quark-gluon plasma. Real-time processes can serve as the signatures for 
a quark-gluon plasma at intermediate coupling. While lattice gauge theory can be used 
to calculate thermodynamic properties Il54l[53l , its application to dynamical quantities 
currently suffers from large systematic errors Il55ll56l . 

There is another possible interpretation of the failure of the conventional perturba- 
tion series. It could simply be a signal that we are using the wrong degrees of freedom. 
Naive perturbation theory is an expansion around an ideal gas of massless quarks 
and gluons. This generates infrared divergences that must be rendered finite either by 
resumming infinite classes of diagrams or by nonperturbative methods. While such a 
procedure gives a well-defined weak-coupling expansion, in practice the coefficients 
seem to be too large for the expansion to be of any use. It is possible that another choice 
for the degrees of freedom would generate diagrams with better infrared behavior and 
successive approximations with better convergence properties. 
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The high-temperature limit of QCD provides a clue as to what those degrees of 
freedom might be. In this limit, quarks and gluons are quasiparticles with temperature- 
dependent masses Il57ll58ll47l . Furthermore, quarks and gluons acquire additional 
propagating degrees of freedom: in addition to the two usual transverse polarization 
modes of the gluon, there is a collective mode with longitudinal polarization called 
the plasmon; in addition to the two usual spin states of a quark, there is a collective 
mode with two spin states called the plasmino. The quasiparticle mass of the gluon is 
also intimately tied to the screening properties of the plasma. Chromoelectric fields are 
screened by the Debye mechanism beyond a screening length of ~ 1/wJd where m^, 
is the gluon quasiparticle mass. Oscillating chromomagnetic fields are also screened, 
with a screening length that scales like {mjjU))^^^'^, where co is the frequency. At 
very low frequencies (a; of order oc^T), nonperturbative effects take over, so that static 
chromomagnetic fields have a screening length of order 1/ {ocsT). 

Quasiparticle masses, collective modes, and screening are all tied together by 
gauge invariance. The problem is therefore how to incorporate plasma effects into the 
perturbation expansion for QCD while preserving gauge invariance. This problem was 
solved at leading order in ^ by Braaten and Pisarski I1T6| . They developed a method 
called hard-thermal-loop (HTL) resummation for summing all Feynman diagrams that 
are leading order in g for amplitudes involving soft external momenta of order gT. 
This method can be used to systematically calculate higher order corrections as an 
expansion in powers of g. 

As one step further for resumming graphs, Andersen, Braaten and Strickland 
introduced hard-thermal-loop perturbation theory (HTLpt) ||20| , which is essentially a 
reorganization of the conventional perturbation expansion for QCD that selectively 
resums higher order effects related to quasiparticles and screening. 

3.2 Formalism 

HTLpt is formulated in Minkowski space, therefore it applies to both thermodynamics 
and real-time dynamics straightforwardly. We use pure-glue QCD next as an example 
to show the formalism of the theory. Note that all the discussions here apply equally 
to the case with quarks. The Minkowskian Lagrangian density that generates the 
perturbative expansion for pure-glue QCD can be expressed in the form 

I 

^QCD = "2^1 (GffyG^") + £gf + £ghost + A-^QCD/ (3.2.1) 

where G^y = dp, Ay — dyA^f — ig[Ap, Ay] is the gluon field strength and A,, is the gluon 
field expressed as a matrix in the SU(Nc) algebra. The ghost term >Cghost depends on 
the choice of the gauge-fixing term £gf . 
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The perturbative expansion in powers of g generates ultraviolet divergences. The 
renormalizability of perturbative QCD guarantees that all divergences in physical 
quantities can be removed by renormalization of the coupling constant ocg = g^/ATZ. 
There is no need for wavefunction renormalization, because physical quantities are 
independent of the normalization of the field. There is also no need for renormalization 
of the gauge parameter, because physical quantities are independent of the gauge 
parameter. 

Hard-thermal-loop perturbation theory is a reorganization of the perturbation 
series for thermal gauge theories with the Lagrangian density written as 



^ = (-Cqcd + -^htl) 



AW . (3.2.2) 



The HTL-improvement term appearing above is (1 — S) times the isotropic HTL effec- 
tive action which generates all HTL n-point functions Il33ll 



1- - .2^Jn /ryL\ 



CnjL = "2(1- ^)"^DTr I^G,, ^ j^-^ G^^^j , (3.2.3) 

where D^, is the covariant derivative in the adjoint representation, y'' = (1, y) is a 
light-like four- vector, and {■ ■ ■)y represents the average over the directions of y. The 



term (3.2.31 has the form of the effective Lagrangian that would be induced by a 
rotationally-invariant ensemble of charged sources with infinitely high momentum. 
Note that the covariant derivatives in the denominators make the HTL-improvement 
terms gauge invariant by modifying all n-point functions self-consistently. The param- 
eter mo can be identified with the Debye screening mass. 

HTLpt is defined by treating ^ as a formal expansion parameter Il59ll *. Physical 
observables are calculated in HTLpt by expanding them in powers of S, truncating 
at some specified order, and then setting S = 1. This defines a reorganization of the 
perturbation series in which the effects of the term in (3.2.31 are included to all 



orders but then systematically subtracted out at higher orders in perturbation theory 



by the Sm^ term in (3.2.3 L If we set S = 1, the Lagrangian (3.2.2 ) reduces to the QCD 



Lagrangian ( 3.2.1[ |; while the free Lagrangian, which reads £qcd + -^htl/ is obtained 



by setting 3 = and describes gluon quasiparticles with screening masses mo. 

We stress here that the S expansion is equivalent to loop expansion, i.e. leading 
order (LO) S expansion is one loop, next-to-leading order (NLO) is two loops, next-to- 
next-to-leading order (NNLO) is three loops, and so on. If the expansion in 3 could be 
calculated to all orders, the final result would not depend on otd when we set 3 = 1. 
However, any truncation of the expansion in 3 produces results that depend on m^. 



* For applications of this so-called linear delta expansion other than HTLpt, please see Ref . l60l for a 
broad but far from complete list. 
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Some prescription is required to determine ntu as a function of T and a. In the next two 
chapters we will discuss different mass prescriptions. 

The HTL perturbation expansion generates ultraviolet divergences. In QCD pertur- 
bation theory, renormalizability constrains the ultraviolet divergences to have a form 
that can be cancelled by the counterterm Lagrangian ACqcb- We will demonstrate 
that renormalized perturbation theory can be implemented by including a counterterm 



Lagrangian A£htl among the interaction terms in (3.2.21. There is no proof that the 



HTL perturbation expansion is renormalizable, so the general structure of the ultra- 
violet divergences is not known; however, it was shown in previous papers ||2T] that 
it was possible to renormalize the NLO HTLpt prediction for the free energy of QCD 
using only a vacuum counterterm, a Debye mass counterterm, and a fermion mass 
counterterm. In this dissertation we will show that renormalization is also possible at 
NNLO. 

The free Lagrangian in general covariant gauge is obtained by setting S = 
in( |322l l: 

1, 



C 



free 



-Tr {dj.Aydf'A^ - 3j,A,9Mf') 
1 



^Tr 



(3.2.4) 



The resulting propagator is the HTL gluon propagator and the remaining terms in ( 3.2.2| 
are treated as perturbations. The propagator can be decomposed into transverse and 
longitudinal pieces which in Minkowski space are given by 



Ar(p) 
Al(p) 



Hr(p) ' 
1 



nL(p) ' 



(3.2.5) 
(3.2.6) 



where rip 



nf — ■pf'{n--p/p^) withn = (1,0) being the vector that specifies the thermal 
rest frame, and Ylj and YIl are the transverse and longitudinal self -energies, respec- 
tively, and read 



m 



D 



id 



l)nl 



m 



Mi-'^°(P'-p)] ' 



(3.2.7) 
(3.2.8) 



with T^^{p, —p) defined in ( A.1.15 1. Note that there are also HTL vertex corrections 
which are given by similar but somewhat more complicated expressions which can be 
found in Refs. lISllZllllll. 
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4 




Figure 3.1: Longitudinal and transverse dispersion relations. This figure is adapted from 
Ref. ESI. 

As mentioned above, HTLpt is a systematic framework for performing calculations 
in thermal gauge theories which is gauge invariant by construction. It systematically 
includes several physical effects of the plasma such as the propagation of massive 
particles, screening of interactions, and Landau damping. We briefly comment on these 
issues next. 



3.2.1 Massive quasiparticles 

The HTL self-energies are included in the zeroth order propagators which results in the 
resummed HTL propagator for gluons. The dispersion relations for the transverse and 
longitudinal gluonic degrees of freedom are determined by locating the zeros of the 
inverse propagator which gives the following two equations 

cv^-k^ -nT{coT,k) = 0, (3.2.9) 
k^ + nL{cvi,k) = 0. (3.2.10) 



The dispersion relations for transverse and longitudinal gluons are shown in Fig. 3.1 
As can been seen from this Figure both modes approach a constant in the limit of small 
momentum and approach the light-cone in the limit of large momentum. 
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3.2.2 Screening 

HTLpt also includes screening of interactions which can be seen by examining the 
static limit of the HTL propagators. For instance, the inclusion of the longitudinal 
self-energy changes the Coulomb potential of two static charges in the plasma to a 
Yukawa potential: 

\imAL{cv,k) = r^^, (3.2.11) 

This result shows that chromoelectric fields are screened on a scale r ~ wj^^- Like- 
wise, the screening of long wavelength chromomagnetic fields is determined by the 
transverse propagator for small frequencies 

~ p2^-n 2 (3.2.12) 

This shows that there is no screening of static magnetic fields meaning that the magnetic 
mass problem of non-Abelian gauge theories at high temperature is not solved by HTL 
resummation. However, HTL resummation does give dynamical screening at a scale 
r ~ (^ni^cv) ' . Note that the divergences associated with the absence of static magnetic 
screening do not pose a problem until four-loop order; however, at four loops the lack 
of static magnetic screening gives rise to infrared divergences that cause perturbation 
theory to break down. 



3.2.3 Landau damping 

The transverse and longitudinal HTL self-energies also contain the physics of Landau 
damping. Landau damping represents a transfer of energy from the soft modes to the 
hard modes for spacelike momentum. This can be seen from the analytic structure of 



the self-energies given by Eqs. (2.2.18 1 and (2.2.19 1. Because of the logarithms appearing 



in these functions there is an imaginary contribution to the self -energies for —k<co<k 
which gives the rate of energy transfer from the soft to hard modes. Note that ignoring 
this contribution leads to gauge variant and unrenormalizable results. 



3.3 Technicalities 

Calculations in HTLpt are much more difficult than in ordinary perturbative QCD, 
because the Feynman rules are more complicated. In spite of the complexity of the 
Feynman rules, calculations do appear to be tractable with the help of the following 
techniques. Due to the fact that only thermodynamic quantities are considered in this 
dissertation, from now on we switch the discussions to EucUdean space for convenience. 
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3.3.1 Mass expansion 

The calculation of the free energy in HTLpt involves the evaluation of vacuum diagrams. 
In Refs. Il20[|2n , the free energy was reduced to scalar sum-integrals. The one-loop free 
energy were evaluated exactly by replacing the sums by contour integrals, extracting 
the poles in e, and then reducing the momentum integrals to integrals that were at most 
two-dimensional and could therefore be easily evaluated numerically. Evaluating two- 
loop free energy exactly would involve the evaluation of five-dimensional numerical 
integrals which turned out to be intractable. Therefore attacking the third loop in this 
way is hopeless. 

The fact that niu ~ gT suggests that niu/T can be treated as an expansion param- 
eter of order g in terms of which the sum-integrals can be further expanded IISTl . It 
was shown that the first few terms in the wjd/T expansion of the sum-integrals gave 
a surprisingly accurate approximation to the exact result ||20l|3ll. We will adopt this 
mass expansion trick in the calculation of three-loop HTL free energy in the next two 
chapters. We will carry out the mu/T expansion to high enough order to include all 
terms through order g^ if niD / T is taken to be of order g. The two-loop approximation 
will be perturbatively accurate to order g'^ and the three-loop approximation accurate 
to order g^. We demonstrate next how the mass expansion works by using the simplest 



example of one-loop photon diagram (la) in Fig. 4.2 



The expression of the one-loop photon diagram (la) in Fig. 4.2 after taking into 
account the ghost contribution is. 



^1 



la 



l^^iid - 1) log [-At(P)] +log A,(P)} . 



(3.3.1) 



After plugging in (3.2.51 and (3.2.61 for Aj/i with Ylj/i defined in (3.2.71 and ( 3.2.8[ l 
and expanding to second order in m^,, the hard contribution from (la) reads. 



(h) 



d-14 



14 



^^^log(P^) + ^^^log(,^) + ^^^i,-^ 
r 1 d 2 



1 d 

+ 



p4 p4 p2p2 p4 



^.d ^ 2 _ d ,^,7 



p2p2 



(3.3.2) 



Note that the integrands log(p^) and have no scale, so the corresponding sum- 
integrals vanish in dimensional regularization. Finally, the hard contribution from (la) 
becomes 



■J 1 



la 



A{d-l)fp 



2"fpP2 
1 2 2d 



pi p2p2 p4 p2p: 



^ 2 _ d , _ , 9 
Tp + -:^rp + ^{Tpr 



, (3.3.3) 
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with sum-integrals listed in App.|B] All the other diagrams are to be evaluated in the 
same spirit. 

3.3.2 Simplified 5 expansion 



We have introduced the \ — 5 description in ( 3.2.3) . The purpose of doing so is to 



distinguish interactions from the free part in the Lagrangian by associating every 



interaction term with a label 5. The subtracted term in (3.2.21 generates self- 
energy and vertex insertions that systematically eliminate the effects of the added rn^ 
term from lower orders. The number of the diagrams with self-energy and vertex 
insertions grows exponentially as we go to higher and higher orders in the 5 expansion. 
Therefore evaluating each diagram individually would become hopeless at higher loop 
order. Since all the self -energy and vertex insertions originate from the (1 — S)m'^ 
term in ( 3.2.3| , the diagrams with self -energy and vertex insertions can be obtained by 
substituting rn^ — ?> (1 — 5)nv^ in the original diagrams and expanding to appropriate 
order in 5. The 5 expansion for hard contributions is trivial since the mo dependence in 
hard modes only enters as multiplicative factors which are of even order in rriD. The 5 
expansion for soft contributions are much more involved due to the fact that nij^ also 
appears in denominators for the soft contributions. We use again the one-loop photon 



diagram (la) in Fig. 4.2 next to show how to carry out the 5 expansion. 



The soft contribution of the one-loop photon diagram (la) in Fig. 4.2 reads. 



•^1? = It I \og{p^ + ml) . (3.3.4) 



P 



After substituting — ?► (1 — ^)m^ and expanding to order S'^ to include all terms 
through g^, we obtain 



l^log [p^ + {l-S)ml] 



]t [ log (p^ + ml) - ^-mlT [ - ^\%T 



2 Jp ' u, ^ u J^pi + rnl 4 " yp(p2 + ^2^)2 

from which the one-loop photon diagram with one and two self -energy insertions. 



i.e. (2c) in Fig. 4.2 and (3h) in Fig. |4.3[ are generated systematically. In the following 
chapters we will show that with the help of 3 expansion, the evaluation of the diagrams 
with self-energy and vertex insertions becomes incredibly simple and straightforward. 
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Chapter 4 

QED Thermodynamics to Three 
Loops 

The thermodynamics of QED is studied in this chapter using the hard-thermal-loop 
perturbation theory reorganization of finite-temperature gauge theory. We calculate 
the free energy through three loops by a dual expansion mmi)/T,mf/T and e^, where 
niD and mjr are thermal masses and e is the coupling constant. The results demonstrate 
that the hard-thermal-loop perturbation reorganization improves the convergence of 
the successive approximations to the QED free energy at large coupling, e ~ 2. The 
reorganization is gauge invariant by construction, and due to cancellation among 
various contributions during renormalization, we obtain a completely analytic result 
for the resummed thermodynamic potential at three loops. This chapter is based on: 
Three-loop HTLfree energy for QED, J. O. Andersen, M. Strickland and N. Su, Phys. Rev. 
D 80, 085015 (2009). 



4.1 Introduction 

The weak-coupling expansion of the QED free energy is known to order e^ Il5ll6ll9l-[TT1. 



In Fig. 4.1 we show the successive perturbative approximations to the QED free energy. 
As can be seen from this figure, at couplings larger than e ~ 1 the QED weak-coupling 
approximations also exhibit poor convergence which is as bad as its counterpart in 
QCD. 

In spite of the complexity of the Feynman rules, calculations with HTLpt do appear 
to be tractable. Andersen, Braaten, Petitgirard, and Strickland have demonstrated this 
by calculating the next-to-leading order (NLO) free energy for QCD ||2T1 . Although 
their results showed striking improvement of convergence comparing to the naive 
weak-coupling expansion, there were still problems remained at two-loop order. First, 
both the leading order (LO) and NLO free energies have a wrong curvature below 2Tc. 



32 



CHAPTER 4. QED THERMODYNAMICS TO THREE LOOPS 




0.5 1 1.5 2 

e(2jtT) 

Figure 4.1: Successive perturbative approximations to the QED pressure (negative of the 
free energy). Each band corresponds to a truncated weak-coupling expansion accurate 
to order e^, e^, e*, and e^, respectively. Shaded bands correspond to variation of the 
renormalization scale }i between tzT and AtzT. 

Instead of going down the results rise up towards Tc. This is due to the fact that the 
truncation order in the dual expansion was and the NLO approximation is only 
perturbatively accurate to order g^. The missing g^ and g^ terms will enter at three loops. 
Second, in the NLO renormalization only vacuum and mass counterterms were needed, 
therefore the self -consistent running coupling could not be derived systematically from 
the calculation and it had to be added by hand in the results. The coupling constant 
renormalization also enters at three-loop order. Therefore it is clear that in order to 
complete the calculation, we need to attack the third loop. However, comparing to 
Abelian case, a direct three-loop non-Abelian calculation might cause unnecessary 
complications which should not be the main concern, we therefore decided to use QED 
as a test case to develop the necessary techniques for attacking QCD. 



4.2 HTL perturbation theory 

The Lagrangian density for massless QED in Minkowski space is 

^QED = -^F^,vF^"' + iipjt'Df,xp + Cgi + Cgi, + ACQED. (4.2.1) 

Here the field strength is F^^' = 9''A^ — d^'A^ and the covariant derivative is D'' = 
3'' + zeA^. The ghost term Cg^ depends on the gauge-fixing term Cg(. In this chapter 
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we choose the class of covariant gauges where the gauge-fixing term is 

^gf = -^i^f^'f ' (4-2.2) 

with ^ being the gauge-fixing parameter. In this class of gauges, the ghost term decou- 
ples from the other fields. 

The HTLpt Lagrangian density for QED is written as 

L = (/:qed + Cutl) ^ + A/:htl , (4.2.3) 
where the HTL-improvement term reads 

W = -4(1 - ^)^hF,, ( (f^y ^'p + - '^^ ™/'^'^" 1" ' ^^-^-^^ 

with the parameter wjd identified with the Debye screening mass and the parameter 
m f identified as the induced finite-temperature electron mass. 

Although the renormalizability of the HTL perturbation expansion has not yet 
been proven, the renormalization was achieved at NLO for the free energy of QCD 
using only a vacuum energy counterterm, a Debye mass counterterm, and a quark 
mass coimterterm Il2l1l . In this chapter we will show that this is also possible at next- 
to-next-to-leading order (NNLO) with the introduction of a new coupling constant 
counterterm which coincides with its perturbative value at zero temperature giving 
rise to the standard one-loop running. The necessary counterterms are 



2 



3Aa = Nf^S^ , (4.2.5) 
Ami = (-^ + O{3^oc^)) (1 - S)ml , (4.2.6) 



Am} = (-^ + 0{3^cc^) ) (1 - 3)m} , (4.2.7) 



As discussed in Section 3.2 a prescription is required to determine mu and as a 
function of T and a when truncating at a finite order in ^. As one possibility we will treat 
both as variational parameters that should be determined by minimizing the free energy. 
If we denote the free energy truncated at some order in (5 by Q(r, ci,mD,mf,ii,5), our 
prescription is 

—n{T,ci,mD,mf,^,5 = l) = 0, (4.2.9) 
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(la) 



(2b) 




Figure 4.2: Diagrams contributing through NLO in HTLpt. The undulating lines are 
photon propagators and the solid lines are fermion propagators. A circle with a n indicates 
a photon self -energy insertion and a circle with a Z indicates a fermion self-energy insertion. 
All propagators and vertices shown are HTL-resummed propagators and vertices. 



dm 



-Q(r, OL,mi:,,mf, }i,5 



. 



(4.2.10) 



Since Cl{T,(x,mD,mj,}i,S = 1) is a function of the variational parameters niD and 
m f, we will refer to it as the thermodynamic potential. We will refer to the variational 



equations (4.2.91 and (4.2.101 as the gap equations. The free energy J-" is obtained by 
evaluating the thermodynamic potential at the solution to the gap equations ( |4.2.9| and 



(4.2.10 1. Other thermodynamic functions can then be obtained by taking appropriate 
derivatives of T with respect to T. 



4.3 Diagrams for the thermodynamic potential 



In this section, we list the expressions for the diagrams that contribute to the ther- 
modynamic potential through order 5^, aka NNLO, in HTL perturbation theory. The 
diagrams are shown in Figs. 4.2[ 4.3 and 4.4 Because of our dual truncation in m^, 
m f, and e the diagrams listed in Fig. 4.4 do not contribute to our final expression so we 
will not explicitly list their integral representations. The expressions here will be given 
in Euclidean space; however, in Appendix |A] we present the HTLpt Feynman rules in 
Minkowski space. 

The thermodjmamic potential at leading order in HTL perturbation theory for QED 
with Nf massless electrons is 



(4.3.1) 



Here, J-'ia is the contribution from the photons 



l^^iid - 1) log [-At(P)] +log A,(P)} 



(4.3.2) 
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(3g) 



(3h) 



QD ^ ^ cxy 



(3c) 



(3d) 



(3e) 



(3f) 



Oooocx) 



(3i) 



(3j) 



(3k) 



(31) 



Figure 4.3: Diagrams contributing to NNLO in HTLpt through order e^. The undulating 
lines are photon propagators and the solid lines are fermion propagators. A circle with a Tl 
indicates a photon self-energy insertion and a circle with a Z indicates a fermion self -energy 
insertion. The propagators are HTL-resummed propagators and the black dots indicate 
HTL-resummed vertices. The lettered vertices indicate that only the HTL correction is 
included. 




Figure 4.4: Diagrams contributing to NNLO in HTLpt beyond order e^. The diagrams 
in the first line above first contribute at order and the second line at order e^. The 
undulating lines are photon propagators and the solid lines are fermion propagators. All 
propagators and vertices shown are HTL-resummed propagators and vertices. 



The transverse and longitudinal HTL propagators At{P) and (P) are given in ( A.12.2 
and ( A.12.3| |. The electron contribution is 



-^ib = ^ logdet[P-E(P)] . 



{P} 

The leading-order vacuum energy counterterm Aq^^o is given by 

1 



(4.3.3) 



(4.3.4) 



The thermodynamic potential at NLO in HTL perturbation theory can be w^ritten 



as 



OnLO = nLO + N/(J'2a + -7^2b+-^2d)+-^2c + Al<fo 
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+ Aimfj^^nLo + Aim^^OLo, (4.3.5) 

where AiSq, Aim^, and Aim^ are the terms of order 3 in the vacuum energy and mass 
counterterms: 

= "64^"^'^' ^^-^-^^ 
Aiml = Nf£^ml, (4.3.7) 

Aim^ = --^m^. (4.3.8) 

The contributions from the two-loop diagrams with electron-photon three- and four- 
point vertices are 

^2a = Tx\S:^{V,Q,R)S{Q)Y\V,Q,R)S{R)\A^\V) , (4.3.9) 

-7"2b = ^e'f Tr[r?'^(P,-P,Q,Q)S(Q)]A?'''(P), (4.3.10) 
2 rp{Q} 

where R = Q — P. The contribution from the HTL photon counterterm diagram with a 
single photon self-energy insertion is 

^2c = ^^^n'''^(P)Af^(P) . (4.3.11) 

The contribution from the HTL electron counterterm diagram with a single electron 
self-energy insertion is 

^2d = -t Tr [E(P)S(P)] . (4.3.12) 

r {p} 

The role of the counterterm diagrams (2c) and (2d) is to avoid overcounting of diagrams 
when using effective propagators in (la) and (lb). Similarly, the role of counterterm 
diagram (3k) is to avoid overcounting when using effective vertices in (2a). 

The thermodynamic potential at NNLO in HTL perturbation theory can be written 

as 

f^NNLO = i^NLO + i J^3c + -^Sj ) + Nf ( + Jsb + J^3d + -^3e + -^3f + J^3g + ^3i 

+ -^sk + -^ai) + J^3h + M£o + A2m|)^^OLo + A2mJ^OLo 

+2-((4pH('^'"^)'^^'^^- ^'-'"'^ 
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where I^iEq, A2m^, l^iw^j, and Aia are terms of order 5^ in the vacuum energy, mass 
and coupling constant counterterms: 



^2ml = -Nj^ml, (4.3.15) 
la 

Aam^. = - — ml. (4.3.16) 
47re ^ 

The contributions from the three-loop diagrams are given by 

-Fsa = Tr[r''(-P,Q-P,Q)S(Q)r^(Q-R,Q,R)S(]^)r^(P,i?,K-P) 

xS{R - P)r^(-Q + R - P, Q - P)S(Q - P)] A''''(P)A'*^(Q - R) ,(4.3.17) 

J^3b = Tr [r?^(P,P + Q,Q)S(Q)r^(-R + Q,Q,R)S(R)r'*(R-Q,R,Q) 

xS(Q)r(-P,Q,P + Q)S(P + Q)] A?"^(P)A«^(R - Q) , (4.3.18) 

J'3c = -]e^i Tr [r?'(P,P + Q,Q)S(Q)r^(-P,Q,P + Q)S(P + Q)" 

xTr [r (-P, R,P + R)S(P + R)r'*(P,P + R,R)S{R)] A''^(P)A'^^(P) , (4.3.19) 

^31 = -^e^^ Tr [r'^^(P,-P,R,R)S(R)l A'^'^(P)A^^(P) 
2 Tp{QR} L J 

xTr P(P,P + Q,Q)S(Q)r(-P,Q,P + Q)S(P + Q)] . (4.3.20) 

The contributions from the two-loop diagrams with electron-photon three- and 
four-point vertices with an insertion of a photon self-energy 

^3d = -le^i Trr(P,Q,R)S(Q)r'^(P,Q,R)S(R)]A?'«(P)n«^(P)A^^(P), 
2 JP{Q} 

(4.3.21) 

-Fsf = -le^i Tr[r'"'(P,-P,Q,Q)S(Q)]A?'«(P)n«^(P)A^''(P), (4.3.22) 

2 rp{Q} 

where R = Q-P. 

The contributions from the two-loop diagrams with the electron-photon three and 
four-point vertices with an insertion of an electron self-energy are 



-Fse = -e'^^^^^A«^(P)Tr [r«(P,Q,R)S(Q)Z(Q)S(Q)r^(P,Q,R)S(R)] , 

(4.3.23) 

-^3g = -h^i A^^(P)Tr[P"'(P,-P,Q,Q)S(Q)E(Q)S(Q)] , (4.3.24) 



P{Q} 

where R = Q-P. 
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The contribution from the HTL photon counterterm diagram with two photon 
self-energy insertions is 

-^3h = -^i ^''''(P)A''^(P)^''^^(P)A^^'(P) . (4.3.25) 

The contribution from HTL electron counterterm with two electron self-energy 
insertions is 

-^3i = It Tr[S(P)E(P)S(P)E(P)] . (4.3.26) 
2T{p} 

The remaining three-loop diagrams involving HTL-corrected vertex terms are 
given by 



-^3k = eW./ Tr[r''(P,Q,R)S(Q)r^(P,Q,R)S(R)] A'"^(P), (4.3.27) 
^Tp{q} 

-^31 = -h^t Tr[P^^P,-P,Q,Q)S{Q)]A^^P) , (4.3.28) 
2 TpiQ} 



where is the HTL correction term given in Eq. ( A.8.2 '. Note also that diagram (31) is 
the same as (2b) since there is no tree-level electron-photon four-vertex. 

In the remainder of this chapter, we work in Landau gauge (^ = 0), but we em- 
phasize that the HTL perturbation theory method of reorganization is gauge-fixing 
independent to all orders in S (loop expansion) by construction. 



4.4 Expansion in the mass parameters 

In this section we carry out the mass expansion for all the diagrams listed in the last 
section to high enough order to include all terms through order if mo and tUf are 
taken to be of order e. The NLO approximation will be perturbatively accurate to order 
and the NNLO approximation accurate to order e^. 

The free energy can be divided into contributions from hard and soft momenta. 
In the one-loop diagrams, the contributions are either hard (h) or soft (s), while at 
two-loop level, there are hard-hard (hh) and hard-soft (hs) contributions. There are no 
soft-soft (ss) contributions since one of the loop momenta is fermionic and always hard. 
At three loops there are hard-hard-hard (hhh), hard-hard-soft [hhs), and hard-soft-soft 
(hss) contributions. There are no soft-soft-soft (sss) contributions, again due to the 
hard fermionic lines. 

In the process of the calculation we will see that there are many cancellations 
between the lower-order HTL-improved diagrams and the higher-order HTL-improved 
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counterterm diagrams. This is by construction and is part of the systematic way in 
which HTLpt converges to the known perturbative expansion. For example, one can see 
that diagrams (2c) and (3h) subtract out the modification of the hard gluon propagator 
due to the HTL-improvement of the propagator in diagram (la). Likewise, one expects 
cancellations to occur between diagrams (lb), (2d) and (3i); (2a), (3d), (3e) and (3k); and 
(2b), (3f), (3g), and (31). Below we wiU explicitly demonstrate how these cancellations 
occur. 

4.4.1 One-loop sum-integrals 
Hard contribution 

For hard momenta, the self -energies are suppressed hy nio/T and nif/T relative to the 
inverse free propagators, so we can expand in powers of T1t{P), TLi{P), and E(P). 
For the one-loop graph (la), we need to expand to second order in m^: 



la 



= T^H 

45 24 



l+(2 + 2||-^ 16 

.2 



V47rT/ 



2e 



2 T-'2 



(4.4.1) 



The one-loop graph with a photon self -energy insertion (2c) has an explicit factor of 
and so we need only to expand the sum-integral to first order in m^: 



2c 



1 

24 



+ 



1 

1 /I 



le 



--7 + 1^^ 



AtzT) 



le 



(4.4.2) 



The one-loop graph with two photon self-energy insertions (3h) must be expanded to 



zeroth order in m|,: 



3h 4(d-l)' 



>2P 
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1 ^ . 27r2 
--7 + 2, + — 



128712 

The sum of Eqs. (4.4.1 |-( 4.4.3| l is very simple 
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(4.4.3) 



(4.4.4) 



This is the free er\ergy of an ideal gas of photons. 

The one-loop graph (lb) needs to be expanded to second order in mjr-. 
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(4.4.5) 



The one-loop fermion loop with a fermion self -energy insertion (2d) must be expanded 
to first order in m ^ : 
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(4.4.6) 



The one-loop fermion loop with two self -energy insertions (3i) must be expanded to 
zeroth order in 
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3i 
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(4.4.7) 



The sum of Eqs. (4.4.5 1-( 4.4.7| | is particularly simple: 
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(4.4.8) 
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This is the free energy of an ideal gas of a single massless fermion. 
Soft contribution 

The soft contributions in the diagrams (la), (2c), and (3h) arise from the Pq = term 
in the sum-integral. At soft momentum P = (0, p), the HTL self-energy functions 
reduce to Ht(P) = and nL(P) — The transverse term vanishes in dimensional 
regularization because there is no momentum scale in the integral over p. Thus the soft 
contributions come from the longitudinal term only and read 
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It J log(p2 + m|,) 
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(4.4.9) 



p (p2 + m2,)2 



32 TT 



(4.4.10) 



(4.4.11) 



Note that we have kept the terms through order e in Eqs. ( |4.4.9[ l and ( 4.4.10[ | as they 
are required in the calculation of the contributions due to counterterms. There is 
no soft contribution from the leading-order fermion term ( 4.3.3[ l or from the HTL 
counterterms ( 4.3.12| | and ( 4.3.26| |. 



4.4.2 Two-loop sum-integrals 

For hard momenta, the self-energies are suppressed by mo/^ arid mj IT relative to the 
inverse free propagators, so we can expand in powers of Hj, H^, and E. 



{hh) contribution 



The (Ml) contribution from (4.3.91 and (4.3.10 1 was calculated in Ref. ||2T1 and reads 

1 
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(4.4.12) 



Consider next the ihh) contribution from ( 4.3.21[ | and (4.3.221. The easiest way to 
calculate this term, is to expand the two-loop diagrams (2a) and (2b) to first order in 
m^. This yields 



(hh) 
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(4.4.13) 



We also need the {hh) contributions from the diagrams (3e), (3g), (3k), and (31) The first 
two diagrams are given by ( |4.3.23| , ( 4.3.24| , while the last remaining ones are given by 
(4.3.271 and ( 4.3.28[ |. The easiest way to calculate these contributions is to expand the 
two-loop diagrams (2a) and (2b) to first order in mr. This yields 
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+ 8.97544 
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The sum of the terms in ( 4.4.12| |-( 4.4.14| is very simple 
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(4.4.14) 



(4.4.15) 



This is the two-loop contribution from the perturbative expansion of the free energy in 
QED. 



(hs) contribution 

In the (hs) region, the momentum P is soft. The momenta Q and R are always hard. 
The function that multiplies the soft propagator A7-(0, p), Al(0, p), or Ax(0, p) can be 
expanded in powers of the soft momentum p. The soft propagators A7-(0, p), Al(0, p), 
and Ax(0, p) are defined in Eqs. ( A.2.4 1, ( A.2.5| and ( A.2.10 1, respectively. In the case of 
Ar(0, p), the resulting integrals over p have no scale and they vanish in dimensional 
regularization. The integration measure scales like m^, the soft propagators Ai{0, p) 
and Ax(0, p) scale like and every power of p in the numerator scales like mo. 



The terms that contribute through order e^m^T and e^mjmuT from (|4.3.9|) and (|4.3.10|l 
were calculated in Ref. Il2l1 and read 
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(4.4.16) 



The (hs) contribution from (4.3.211 and (4.3.221 can again be calculated from the 
diagrams (2a) and (2b) by Taylor expanding their contribution to first order in m^. This 
yields 
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(4.4.17) 



We also need the {hs) contributions from the diagrams (3e), (3g), (3k), and (31). Again 
we calculate their contributions by expanding the two-loop diagrams (2a) and (2b) to 
first order in mjr. This yields 
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(4.4.18) 



(ss) contribution 

There are no contributions from the (ss) sector since fermionic momenta are always 
hard. 



4.4.3 Three-loop sum-integrals 
{hhh) contribution 

If all three loop momenta are hard, we can expand the propagators in powers of H^^i/ (P) 
and S(P). Through order e^, we can use bare propagators and vertices. The diagrams 
(3a), (3b), and (3c) were calculated in Refs. ||6l|5l and their contribution is 
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Using the expression for the sum-integrals in Appendix |B| we obtain 
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(4.4.19) 



(4.4.20) 



(hhs) contribution 

The diagrams (3a) and (3b) are both infrared finite in the limit niD — > 0. This implies 
that the contribution is given by using a dressed longitudinal propagator and bare 
vertices. The ring diagram (3c) is infrared divergent in that limit. The contribution 
through is obtained by expanding in powers of self-energies and vertices. Finally, the 
diagram (3j) also gives a contribution of order e^. Since the electron-photon four-vertex 
is already of order e^mjr, we can use a dressed longitudinal propagator and bare fermion 
propagators as well as bare electron-photon three-vertices. Note that both (3c) and (3j) 
are proportional to Nj and so it is more convenient to calculate their sum. One finds 
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Using the expressions for the integrals and sum-integrals listed in Appendixes [B| and [C} 
we obtain 



J + ^+27 + 21og2 



+ 2^ 



AnTJ \2mDJ ^ An 



+ Nf 



N? mlT^ . (4.4.24) 
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{hss) contribution 



The {hss) modes first start to contribute at order e^, and therefore at our truncation 
order the (hss) contributions vanish. 

(sss) contribution 

There are no contributions from the (sss) sector since fermionic momenta are always 
hard. 

4.5 Thermodynamic potentials 

In this section we present the final renormalized thermodynamic potential explicitly 
through order S'^, aka NNLO. The final NNLO expression is completely analytic; 
however, there are some numerically determined constants which remain in the final 
expressions at NLO. 



4.5.1 Leading order 



The complete expression for the leading order thermodynamic potential is given by 
the sum of Eqs. ( 4.4.1[ |, (4.4.5 ), and ( 4.4.9[ l plus the leading vacuum energy countert- 
erm ( |4.3.4| |: 
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where thu, nip and p. are dimensionless variables: 
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(4.5.1) 



(4.5.2) 
(4.5.3) 
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iTZT 



(4.5.4) 



4.5.2 Next-to-leading order 

The renormalization contributions at first order in 3 are 
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(4.5.5) 



Using the results listed in Eqs. ( 4.3. 6[ l, ( 4.3. 7[ |, and ( 4.3.8[ l, the complete contribution 
from the counterterm at first order in 3 is 
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(4.5.6) 



Adding the NLO counterterms ( 4.5. 6[ l to the contributions from the various NLO 
diagrams, we obtain the renormalized NLO thermodynamic potential 
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(4.5.7) 



4.5.3 Next-to-next-to-leading order 

The renormalization contributions at second order in 3 are 
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Using the results listed in Eqs. (4.3.141, ( 4.3. 15| , and ( 4.3. 16| , the complete contribution 
from the counterterms at second order in (5 is 
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Adding the NNLO counterterms (4.5.9 1 to the contributions from the various NNLO 
diagrams, we obtain the renormalized NNLO thermodynamic potential. We note that 
at NNLO all numerically determined coefficients of order drop out and we are left 
with a final result which is completely analytic. The resulting NNLO thermodynamic 
potential is 
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(4.5.10) 



We note that the coupling constant counterterm listed in Eq. ( 4.2.5| coincides with 
the known one-loop running of the QED coupling constant 



d}i 



Nfe^ 



(4.5.11) 



Below we will present results as a function of e evaluated at the renormalization scale 
2tcT. Note that when the free energy is evaluated at a scale different than ji = 2tcT we 
use Eq. ( 4.5.11| to determine the value of the coupling at = 27rr. 
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We have already seen that there are several cancellations that take place alge- 
braically, irrespective of the values of mo and m^. For example the (hh) contribution 
from the two-loop diagrams (2a) and (2b) cancel against the {hh) contribution from the 
diagrams (3d), (3e), (3f), and (3g). As long as only hard momenta are involved, these 
cancellations will always take place once the relevant sum-integrals are expanded in 
powers of m D / T and nif /T. This is no longer the case when soft momenta are involved. 
However, further cancellations do take place if one chooses the weak-coupling values 
for the mass parameters. For example, if one uses the weak-coupling value for the 
Debye mass. 



^{Q} 



1 2q^ 



= ^NfCcT\ (4.5.12) 

the terms proportional to mjr in Onnlo cancel algebraically and HTLpt reduces to the 
weak-coupling result for the free energy through e^. This reduction is by construction 
in HTLpt which also provides a consistency check that in the weak-coupling limit, 
HTLpt coincidences with weak-coupling expansion. 



4.6 Free energy 



The mass parameters mo and nif m hard-thermal-loop perturbation theory are in prin- 
ciple completely arbitrary. To complete a calculation, it is necessary to specify mo 
and as functions of e and T. In this section we will consider two possible mass 
prescriptions in order to see how much the results vary given the two different assump- 
tions. First we will consider the variational solution resulting from the thermodynamic 



potential, Eqs. (4.2.9 ' and (4.2.10 1, and second we will consider using the perturbative 
expansion of the Debye mass ll62l[T0l and the perturbative expansion of the fermion 



mass 



4.6.1 Variational Debye mass 



The NLO and NNLO variational Debye mass is determined by solving Eqs. ( 4.2.9[ l 
and ( |4.2.10[ ) using the NLO and NNLO expressions for the thermodynamic potential, 
respectively. The free energy is then obtained by evaluating the NLO and NNLO 
thermodynamic potentials, ( 4.5. 7[ | and ( 4.5.10[ l, at the solution to the gap equations 
( |4.2.9| and ( 4.2.10| . Note that at NNLO the gap equation for the fermion mass is trivial 
and gives mj = 0. The NNLO gap equation for mu reads 
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(4.6.1) 



In Figs. [43} |4^ and [47| we plot the NLO and NNLO HTLpt predictions for the 

the renormalization scale 



4.7 



free energy of QED with Nj = 1. As can be seen in Fig. 
variation of the results decreases as one goes from NLO to NNLO. This is in contrast to 
weak-coupling expansions for which the scale variation can increase as the truncation 
order is increased. 



-0.95 




1 1.5 

e(2n;T) 



Figure 4.5: NLO HTLpt predictions for the free energy of QED with = 1 and the 
variational Debye mass. Different curves correspond to varying the renormalization scale 
j,i by a factor of 2 around fi — 2nT. 



One troublesome issue with the variational Debye mass is that at NNLO the 
solutions to ( 4.6.1[ l have a small imaginary part. The NNLO gap equation ( 4.6.1[ l is 
quadratic in m^, so it has two complex conjugate solutions for m^. Although the 
solutions are real for large couplings, they become complex when the coupling is 
smaller than a critical value, e.g. the critical value for Nf = 1 and }i = 2nT is 
e{2nT) = 3.38946. For small coupling and finite Nf, the variational Debye masses can 
be expanded as follow: 
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nir 
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n 



v/6 



(4.6.2) 



which reproduces the weak-coupling result (4.6.3 > at leading order, however the term 
becomes imaginary which starts to deviate from the weak-coupling result. It should 
be mentioned that as — )■ oo, the critical value below which mo becomes complex 
goes to zero. We plot the imaginary part of the free energy which results from these 
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Figure 4.6: Real (top panel) and imaginary (bottom panel) parts of NNLO HTLpt predic- 
tions for the free energy of QED with Njr = 1 and the variational Debye mass. Different 
curves correspond to varying the renormalization scale ^ by a factor of 2 around = InT. 



imaginary contributions to the variational Debye mass in Fig. 4.6 (bottom panel). The 
imaginary contributions to the variational Debye mass come with both a positive and 
negative sign corresponding to the two possible solutions to the quadratic variational 
gap equation. The positive sign would indicate an unstable solution while the negative 
sign would indicate a damped solution. These imaginary parts are most likely an 
artifact of the dual truncation at order e^; however, without extending the truncation to 
higher order, it is difficult to say. They do not occur at NLO in HTLpt in either QED or 
QCD. We note that a similar effect has also been observed in NNLO scalar theory 
[301 and Yang-Mills theory ||23l|24l. Because of this complication, in the next subsection 
we will discuss a different mass prescription in order to assess the impact of these small 
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0.5 1 1.5 2 

e(2jtT) 

Figure 4.7: A comparison of the renormalization scale variations between NLO and NNLO 
HTLpt predictions for the free energy of QED with Nf = 1 and the variational Debye 
mass. The bands correspond to varying the renormalization scale ^ by a factor of 2 aroimd 
}i = 2nT. 



imaginary parts. 



4.6.2 Perturbative Debye and fermion masses 

The perturbative Debye and fermion masses for QED have been calculated through 
order ll62l[T0l and Il63l , respectively: 
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(4.6.3) 
(4.6.4) 



Plugging ( |4.6.3[ ) and ( |4.6.4[ ) into the NLO and NNLO thermodynamic potentials, ( |4.5.7| 
and ( |4.5.10| , we obtain the results shown in Fig. |4.8[ The renormalization scale variation 
is quite small in the NNLO result. 



4.6.3 Comparison with the O-derivable approach 

Having obtained the NNLO HTLpt result for the free energy we can now compare the 
results obtained using this reorganization with results obtained within the O-derivable 



approach. In Fig. 4.9 we show a comparison of our NNLO HTLpt results with a 
three-loop calculation obtained previously using a truncated three-loop O-derivable 
approximation ||64l . For the NNLO HTLpt prediction we show the results obtained 
using both the variational and perturbative mass prescriptions. As can be seen from 
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Figure 4.8: A comparison of the renormalization scale variations between NLO and NNLO 
HTLpt predictions for the free energy of QED with = 1 using the perturbative thermal 
masses given in Eqs. 14.6.3 1 and | 4.6.4| . The bands correspond to varying the renormaliza- 
tion scale by a factor of 2 around ji = 2nT. 




0.5 1 1.5 2 

e(2jTT) 

Figure 4.9: A comparison of the predictions for the free energy of QED with Nf = 1 
between three4oop O-derivable approximation 1641 and NNLO HTLpt at = 27rT. 



this figure, there is very good agreement between the NNLO O-derivable and HTLpt 
approaches out to large coupling. The difference between these two predictions at 
e = 2.4 is merely 0.6%. In all cases we have chosen the renormalization scale to be 

]i = iTlT. 



As a further consistency check, in Fig. 4.10 we show a comparison between the 
untruncated two-loop numerical O-derivable approach calculation of Ref. Il65l and 
our NLO HTLpt result using the variational mass. In both cases we have chosen the 
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Figure 4.10: A comparison of the predictions for the free energy of QED with Njr = 1 
between the two-loop 2P1 approximation in Landau gauge Il65l and NLO HTLpt at }i = 
2tzT. 



renormalization scale to be }i = 2nT. From this figure we see that there is a reasonable 
agreement between the NLO numerical O-derivable and NLO HTLpt results; however. 



the agreement is not as good as the corresponding NNLO results shown in Fig. 4.9 
We note that the results of 



were computed in the Landau gauge = 0). As 
detailed in their paper, their result is gauge dependent. Such gauge dependence is 
unavoidable in the 2PI O-derivable approach since it only uses dressed propagators. In 
Ref. Il64l it was explicitly shown that the two-loop O-derivable Debye mass is gauge 
independent only up to order e^, resulting in gauge variation of the free energy at order 
e^. This is in agreement with general theorems stating that the gauge variance appears 
at one order higher than the truncation ||66| . 



4.6.4 QCD free energy at large Nj 

The large Nj limit is achieved by taking Nj to be large while holding e^Nf fixed. The 
large Nf coupling for QED is defined by ^eff = ^\J^f- power counting, it is easily 
to see that in perturbation theory only ring diagrams survive in the large Nf limit, 
which indicates the equivalence of QED and QCD at large Nf. In the large Nf limit, it 
is possible to solve for the 0{N^) contribution to the free energy exactly Il67ll68ll . In 



Fig. 4.11 we plot the NLO and NNLO HTLpt predictions for the free energy at large 
Nf along with the numerical result of Ref. II68I , as well as the perturbative ^^g, and 
newly obtained g^^ Il69l predictions at fi = e^^TzT. The NLO HTLpt result seem to 
diverge from the exact result around ^eff = 2, while the NNLO one from ^gff = 2.8, 
however both of their large coupling behaviors qualitatively fit that of the ntimerical 
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result. 
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Figure 4.11: A comparison of the predictions for the large Ny^ free energy of QED between 
the numerical result from ||68l , NLO and NNLO HTLpt, and perturbative through 
results at p = e^i'/rT. 
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4.7 Conclusions 

In this chapter we calculated the three-loop HTLpt thermodynamic potential in QED. 
Having obtained this we applied two mass prescriptions, variational and perturbative, 
to fix the a priori undetermined parameters and m j that appear in the HTL-improved 
Lagrangian. We found that the resulting expressions for the free energy were the same 
to an accuracy of 0.6% at e = 2.4 giving us confidence in the prediction. We also 
compared the HTLpt three-loop result with a three-loop O-derivable approach 1641 and 
found agreement at the subpercentage level at large coupling. Besides, the large Nf 
HTLpt three-loop result is in reasonable agreement with the exact numerical one Il68l . 

In addition, we showed that the HTLpt NLO and NNLO approximations have 
improved convergence at large coupling compared to the naively truncated weak- 
coupling expansion and that the renormalization scale variation at NNLO using both 
the variational and perturbative mass prescriptions was quite small. Therefore, the 
NNLO HTLpt method result seems to be quite reliable. This is important since, unlike 
the O-derivable approach, the HTLpt reorganization is gauge invariant by construction 
and is formulated directly in Minkowski space allowing it to, in principle, also be 
applied to the calculation of dynamical quantities. 

The renormalization of the three-loop thermodynamic potential required only 
known vacuum energy, mass, and coupling constant counterterms, and the resulting 
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running coupling was found to coincide with the canonical QED one-loop running. 
This provides further evidence that the HTLpt framework is renormalizable despite 
the new divergences which are introduced during HTL improvement. 

Finally, we note that at three loops we could obtain an entirely analytic expression 
for the renormalized NNLO thermodynamic potential. There were a number of can- 
cellations that took place during renormalization which resulted in an expression that 
was independent of any numerically determined subleading coefficients in the sum- 
integrals. With the confidence in the techniques, we are ready to step into non-Abelian 
theories. 
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Chapter 5 



Yang-Mills Thermodynamics to 
Three Loops 

In this chapter, we study the thermodynamics of Yang-Mills theory using the hard- 
thermal-loop perturbation theory in the same spirit of Chapter |4] We show that at three- 
loop order hard-thermal-loop perturbation theory is compatible with lattice results for 
the pressure, energy density, and entropy down to temperatures T ~ 2 — 3 T^. This 
chapter is based on: Gluon Thermodynamics at Intermediate Coupling, J. O. Andersen, 
M. Strickland and N. Su, Phys. Rev. Lett. 104, 122003 (2010), and Three-loop HTL gluon 
thermodynamics at intermediate coupling, J. O. Andersen, M. Strickland and N. Su, JHEP 
1008, 113 (2010). 

5.1 Introduction 

The goal of ultrarelativistic heavy-ion collision experiments is to generate energy 
densities and temperatures high enough to create a quark-gluon plasma. One of 
the chief theoretical questions which has emerged in this area is whether it is more 
appropriate to describe the state of matter generated during these collisions using 
weakly-coupled quantum field theory or a strong-coupling approach based on the 
AdS/ CFT correspondence. Early data from the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Labs indicated that the state of matter created there behaved 
more like a fluid than a plasma and that this "quark-gluon fluid" is strongly coupled 

m. 

In the intervening years, however, work on the perturbative side has shown that 
observables like jet quenching ||2l and elliptic flow ||3| can also be described using 
a perturbative formalism. Since in phenomenological applications predictions are 
complicated by the modeling required to describe, for example, initial-state effects, 
the space-time evolution of the plasma, and hadronization of the plasma, there are 
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Figure 5.1: Weak-coupling expansion for the scaled pressure of pure-glue QCD. Shaded 
bands show the result of varying the renormalization scale by a factor of 2 around 
ji = 2nT. 

significant theoretical uncertainties remaining. Therefore, one is hard put to conclude 
whether the plasma is strongly or weakly coupled based solely on RHIC data. To have 
a cleaner testing ground one can compare theoretical calculations with results from 
lattice quantum chromodynamics (QCD). 

Looking forward to the upcoming heavy-ion experiments scheduled to take place 
at the Large Hadron Collider (LHC) at the European Organization for Nuclear Research 
(CERN) it is important to know if, at the higher temperatures generated, one expects a 
strongly-coupled (liquid) or weakly-coupled (plasma) description to be more appropri- 
ate. At RHIC, initial temperatures on the order of one to two times the QCD critical 
temperature, Tc ~ 190 MeV, were obtained. At LHC, initial temperatures on the order 
of 4 — 5 Tc are expected. The key question is, will the generated matter behave more 
like a plasma of quasiparticles at these higher temperatures. 

As is well known, the weak-coupling expansion for the free energy of SU(3) pure- 
glue QCD fails to converge at phenomenologically relevant temperatures that are being 
created in the colliders HSl lTTI - ITSll . In Fig. |5.1| we show the weak-coupling expansion 
results through order a^''^ from which one sees severe oscillations as higher orders are 
included in the expansion. Equipped with the techniques as well as the confidence 
from the QED calculation in Chapter |4| we are ready to generalize HTLpt to Yang-Mills 
theory. 
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5.2 HTL perturbation theory 

The Lagrangian density for SU(Nc) Yang-Mills theory in Minkowski space is 

Cym = - ^Tr [G^,G'"'] + £gf + + ACyu • (5.2.1) 

Here the field strength is Gf" = 3^'^" - d^A^" - ig[At',A"], with A^' an element of the 
SU(Nc) gauge group. The ghost term Cg^ depends on the gauge-fixing term Cgf. In this 
chapter we choose the class of covariant gauges where the gauge-fixing term is 

1 



(5.2.2) 



with ^ being the gauge-fixing parameter. 

The HTLpt Lagrangian density for Yang-Mills theory is written as 

C = (Cym + Chjl) ^ + A£htl • (5.2.3) 

g^VSg 

The HTL-improvement term is 

W = - ^ (1 - S)mlTr (^G,. ^ ^ G"^^ , (5.2.4) 

where the covariant derivative reads = d^' — igA^. 

Similar to the QED case, as we will show HTLpt for Yang-Mills theory is also renor- 
malizable at NNLO with a coupling constant counterterm, a Debye mass counterterm 
and a vacuum energy counterterm which read 

SAocs = -^<x^,S^ + 0{S^4), (5.2.5) 
iZrce 

Ami = ("I^'^^^ + ^('^'«')) (1 - ^Wd , (5.2.6) 

In the following we will first obtain the thermodynamic potential Q which is a 
function of T, cUs and mo. In order to get the free energy T , some prescription is required 
to determine mo as a function of T and dls. We will discuss several prescriptions in 
Section 1531 
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Figure 5.2: Diagrams contributing through NLO in HTLpt. The spiral lines are gluon 
propagators and the dotted lines are ghost propagators. A circle with a Tl indicates a gluon 
self-energy insertion. All propagators and vertices shown are HTL-resummed propagators 
and vertices. 





3c 





3g 




Figure 5.3: Diagrams contributing to NNLO in HTLpt through order g^. The spiral lines 
are gluon propagators and the dotted lines are ghost propagators. A circle with a Tl 
indicates a gluon self -energy insertion. The propagators are HTL-resummed propagators 
and the black dots indicate HTL-resummed vertices. The lettered vertices indicate that 
only the HTL correction is included. The yellow box in (3f) denotes the insertion of the 



one-loop self -energy defined in Fig. 5.4 



5.3 Diagrams for the thermodynamic potential 

In this section, we list the expressions for the diagrams that contribute to the ther- 
modynamic potential through order S^, aka NNLO, in HTL perturbation theory. The 



diagrams are shown in Figs. |5.2[ and 5.3 A key to the diagrams is given in Fig. 5.4 The 
expressions here will be given in Euclidean space; however, in Appendix |A] we present 
the HTLpt Feynman rules in Minkowski space. 

The thermodynamic potential at leading order in HTL perturbation theory for 
QCD 

Hlo = (N^ - 1) J-ia+ib + Ao^o • (5.3.1) 



61 



CHAPTER 5. YANG-MILLS THERMODYNAMICS TO THREE LOOPS 



1^66 = Dressed Propagator (nj = Self-Energy Insertion 
• = Dressed Vertex © = Vertex Insertion 



□ 




+ my. 



Figure 5.4: Key to the diagrams in Figs. 5.2 and 5.3 



Here, J-'u+ib is the contribution from the gluons and ghost shown on the first line of 



Fig. 5.2 



^la+ib = -^^^{(d-l)log[-Ar(P)]+logAL(P)}. (5.3.2) 



The transverse and longitudinal HTL propagators Aj (P) and (P) are given in ( A.12.2 



and ( A.12.3| |. The leading-order vacuum energy counterterm AqSq was determined in 
Ref. 11201 : 



7\r2 _ 1 



The thermodynamic potential at NLO in HTL perturbation theory can be written 

as 



OnLO = nLO + (N2-l)[J2a + -^2b + -^2c + -?='2d] 

+ Ai£o + Aifnf,^QLo , (5.3.4) 

where AiSq and Aim^ are the terms of order S in the vacuum energy density and mass 
counterterms: 

M2 _ 1 

Aimf, = -^^^s^^D- (5.3.6) 

The contributions from the two-loop diagrams with the three-gluon and four-gluon 
vertices are 



^2a = ^g^i r^'P{P,Q,R)r^^P,Q,R)A>''{P)A^''{Q)AP^{R), (5.3.7) 
12 Tpq 

J'lb = ^g^i rf-'^-(P,-P,Q,-Q)A'-(P)A^'^(Q), (5.3.8) 
8 Tpq 
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where R = —Q — P. The contribution from the ghost diagram is 



The contribution from the HTL gluon counterterm diagram with a single gluon self- 
energy insertion is 

^2d = l^n>^'iP)A^'iP) . (5.3.10) 
The thermodynamic potential at NNLO in HTL perturbation theory can be written 

as 

f^NNLO = f^NLO + (^c ~ 1) [-^Sa + ^3b + J^3c + ^3d + J^3e + ^3f + -^Sg + -^Sh + ^3i 
+ J"3j + J"3k + ^31 + ^3m] + A2£^o + A2m|) ^^Olo + Aim|) ^^Onlo 

+ I ( + ^""^ - D^^^^i'^^ ' (5.3.11) 

where Aiccg, A2inj), and A2^o are terms of order 3^ in the coupling constant, mass, and 
vacuum energy counterterms: 



UN, 



9 llNc 7 

Az^D = — oistnf) , (5.3.13) 



The contributions from the three-loop diagrams are given by 

N2 



N r 

'sa = -^g% T''fHP,Q,-P - Q)A'''{P)At'nQW{P + Q) 
24 TpQR 

xr?'^^(-Q, -R, Q + R)A^''(R)A^^(Q + R) 



^ptrAp^p _^Q_Q_jljl_ p)AP</'(R - P)r^'^"^(-P, P - R) , (5.3.15) 

^3b - f pq^^2(ptwTqW ^-^'-Q'^ + Q^ 

X A«^(P)A^''(P + Q)AT^(Q) , (5.3.16) 

-7^3d = -^E r«^'?'^(P,Q,R,S)rT"^''^^(P,Q,R,S)A«T(P)A^^(Q)A?''^(R)A''^(S), 
48 Tpqr 

(5.3.18) 
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-^3e = T-^''^''{P,Q,R,S)A^f'{P)A>^''{Q)A^'{R)A''f{S)A'\P + Q) 

4 TPQR 

xTf^''^{-P, -Q,P + Q)T^^'^{-P - Q, -R, -S) , (5.3.19) 
-^3f = ^ nf^(P)A^'^(P)n'^^(P)A^?'(P) , (5.3.20) 

= fpQR P4(p + Q)2(p + R)2 ^' (Q)^ '^(K) ' (5.3.21) 

where S = — (P + Q + R) and nf'^(P) is the one-loop gluon self-energy with HTL- 



resummed propagators and vertices defined by the yellow box in Fig. 5.4 



U^-(P) = \Ncg'f P''''^{P,-P,Q,-Q)A^^{Q) + \Ncg'f 1^^^{P,Q,-R 



Q z TQ 

xA-^{Q)r^\P,Q,-R)Af'\R) + Ncg^^^^ , (5.3.22) 

where R = P -|- Q. The contributions from the two-loop diagrams with a single self- 
energy insertion are 



•^3h = -^g^i r-f''{P,Q,-Q-P)r^^\P,Q,-Q-P) 

4 Tpo 



PQ 

xA'''"(P)n'"^(P)A^^(P)A^T(Q)A^^(-Q - P) , (5.3.23) 



^3i = -^/^^^r''^''-(P,-P,Q,-Q)A«^(P)m^(P)A^^(P)A^-(Q) ,(5.3.24) 

^3j = - Y^'fpQ P2 (P + Q)2 ^^HQWiQW^Q) ■ (5.3.25) 
The two-loop diagrams with a subtracted vertex is 

-^3k = ^g^mli V'^P{P,Q,-Q-P)r"'^P,Q,-Q-P) 
6 Tpo 



PQ 

xAf'^P)A^''{Q)AP^{-Q - P) , (5.3.26) 



-^31 = ^8^<i r'"^"'"(P,-P,Q,-Q)A'-(P)A^'^(Q). (5.3.27) 
8 TpQ 



The contribution from the HTL gluon counterterm diagram with two gluon self-energy 
insertions is 

-^3m = -H n^'"{P)A'"'(P)n"NP)A^^'(P) . (5.3.28) 
4~p 

5.4 Expansion in the mass parameter 

In this section we carry out the mass expansion for all the diagrams listed in the last 
section to high enough order to include all terms through order if nio is taken to be 
of order g. The NLO approximation will be perturbatively accurate to order g^ and the 
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NNLO approximation accurate to order g^. 

The free energy can be divided into contributions from hard and soft momenta. In 
the one-loop diagrams, the contributions are either hard (h) or soft (s), while at the 
two-loop level, there are hard-hard {hh), hard-soft (hs), and soft-soft (ss) contributions. 
At three loops there are hard-hard-hard (hhh), hard-hard-soft (hhs), hard-soft-soft 
{hss), and soft-soft-soft (sss) contributions. 

5.4,1 Leading order 
Hard contribution 



For hard momenta, the self-energies are suppressed by mo/T relative to the propaga- 
tors, so we can expand in powers of YIt{P) and Hl(P). 

For the one-loop graphs (la) and (lb), we need to expand to second order in m^: 



(h) 

la+lb 



A{d-1) 



1 

P4 



>2P2 



45 



24 'n^+^fFT'^ 



1 /I ^ . 



(JL_) 



AttT 

2e 

m 



D ■ 



(5.4.1) 



Soft contribution 

The soft contribution in the diagrams (la) and (lb) arises from the Pq = term in the 
sum-integral. At soft momentum P = (0, p), the HTL self -energy functions reduce 
toHT(P) = OandHL(P) — Tit^. The transverse term vanishes in dimensional regu- 
larization because there is no momentum scale in the integral over p. Thus the soft 
contributions come from the longitudinal term only and read 



7r(s) 
•^la+lb 



log (p2 



m 



2e 



12n \2mJ 



(5.4.2) 



We have kept the order e terms in the and m|, terms, respectively in Eqs. (5.4.1 
and( |5A2| since they contribute in the counterterms at next-to-leading order. 
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5.4.2 Next-to-leading order 
Hard contribution 

The one-loop graph with a gluon self-energy insertion (2d) has an explicit factor of 
and so we need only to expand the sum-integal to first order in m^: 



(h) 



2d 



1 iJ 1 



+ 



2{d 
1 

24 



1+ 2 + 2 



p4 p2p2 



1 1 

2d-^Tp + 2-^Tp 



+ 



647r2 



- - 7 + 27e 



\4:TTTJ 



(5.4.3) 



Soft contribution 



The soft contribution from (2d) arises from the Pq = term in the sum-integral. Only 
the longitudinal part TLl{P) of the self-energy contributes and reads 



is) 



2d 



1 2r 



P 



Stt V 2m 



Id 



) [l+2e] 



(5.4.4) 



(hh) contribution 

For hard momenta, the self-energies are suppressed by niu/T relative to the propa- 
gators, so we can expand in powers of H^ and H^. The two-loop contribution was 
calculated in Ref. 11211 and reads 



(hh) 



2a+2b+2c 



PQ 



1 1 



Nc 2 2 

4 J PQ 



2{d-2)^-^ + {d + 2) ^ 



P2^2Q2 



2d 



P-Q 



4p2^-2(rf-l)]!2^rs-(d + l) 

4d^jL^rR + 2dZ^rR 



-Ad 



p2Q2f. 



p2Q2^4 
Tr 



•p2Q2^4 ' P2Q2^4 

Using the expressions for the sum-integrals listed in Appendix |B| we obtain 



(5.4.5) 



(hh) 



2a+2b+2c 



Neds 



1+ (2 + 4|^ \e 



1^ \ 



4e 



47rT 



) 



p4 
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7 
96 



4.621 



NcOis 
371 



47rTy ° 



(5.4.6) 



{hs) contribution 

In the {hs) region, the momentum P is soft. The momenta Q and R are always hard. 
The function that multiplies the soft propagator A7(0, p), Al(0,p) or Ax(0, p) can be 
expanded in powers of the soft momentum p. In the case of A7(0, p), the resulting 
integrals over p have no scale and they vanish in dimensional regularization. The 
integration measure scales like w^,, the soft propagators Al(0,p) and Ax(0,p) scale 
like l/rn^, and every power of p in the numerator scales like mo. The two-loop 
contribution was calculated in Ref. Il2l1l and reads 



(hs) 



2a+2b+2c 



1 



-{d 



2{d 



+ Ncg^mlT [ -Y—^i 
Jp p2 + mfjTQ 



■(^-4)o4 



1 



{d-l){d + 2) 
d 



4(^-1) q4 



(5.4.7) 



In order to facilitate the calculations, it proves useful to isolate the terms that are 
speficic to HTL perturbation theory. After integrating by parts and using the results 
from Zhai and Kastening llTTjl , we can write 



(hs) 



2a+2b+2c 



2 



.^-5. + 16]/T.?,/^-^^^i, 



12 

Nc 
2 



Using the expressions for the integrals and sum-integrals in Appendices |B] and |Cj we 
obtain 



■j-(hs) 

'^2a+2b+2c 



2 in 
11 



1+ 2 + 2 



AnTJ \2mD 

le 



le 



11/1 27 ^ \ NcOis / U \2e/ u ,^ 

-H +27 ^^ M— M— mlT . 5.4.9 

327r 11 J 37r V47rT/ \2mD J ° 



(ss) contribution 

The (ss) contribution was obtained by Braaten and Nieto by a two-loop calculation in 
electrostatic QCD (EQCD) in three dimensions fl2\ . Alternatively, one can isolate the 
(ss) contributions from the two-loop diagrams which were calculated by Arnold and 
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Zhai in Ref. ||5l. In Ref. Il2l1 , this contribution was calculated and agrees with earlier 
results. One finds 



(ss) 



2a+2b+2c 



4 

3 
16 



pq p^{q^ + ml){r^ + ml) 

4e 



+ 3 



3n \2mu 



2 T-2 



(5.4.10) 



We have kept the order e in Eqs. ( 5.4.3| , ( 5.4.4| , ( 5.4.6| l, and ( 5.4.9| since they contribute 
in the counterterms at NNLO. 



5.4.3 Next-to-next-to-leading order 
Hard contribution 

The one-loop graph with two gluon self-energy insertions (3m) is proportional to mjj 
and so must be expanded to zeroth order in m^. 



(h) 



3m 



4(£f-l) -jp 

1 /I ^ ^ 2n^ 



(5.4.11) 



Soft contribution 

The soft contribution from (3m) arises from the Pq = term in the sum-integral. Only 
the longitudinal part Yli{P) of the self-energy contributes and reads 



(s) 



3m 



' 327r 



(p2 



m 



2 )2 

d) 



(5.4.12) 



(hh) contribution 

We also need the {hh) contribution from the diagrams (3h)-(31). We calculate their 
contributions by expanding the two-loop diagrams (2a)-(2c) to first order in m^^. This 
yields 



(hh) 



3h-31 



j^'-it,^ [-2(rf - 1 ) J, ij + 2(rf - 2) ^ ^ + (rf + 2) ^ 

2 



^ 'P^q^Q 

P^Q^r^R^ ^ >P^q^Q^'^ P^Q^r^ P^Q^r^ 



{d + 1] 



^ TR+4d-^-^rR+2dJr^rR 



P^Q^r 
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(hs) contribution 



4.621 



n \4tiTJ ° 



(5.4.13) 



We also need the (hs) contribution from the diagrams (3h)-(31). Again we calculate 
their contributions by expanding the two-loop diagrams (2a)-(2c) to first order in m^. 
This yields 



•^3h-31 



Nc 
2 



12 



^mlT[(f-5d + l6] [ ^ 
■ 'v 

g\d-5)mlT j 

Jo 



2- - ''"■""^^(p2 + m|,)2fQQ4 

Using the expressions in Appendices |B] and |C| we obtain 



(5.4.14) 



[hs) 



3h-31 



TzNcOc, 33/1 59 



4 3n 



6471 \e 33 



+ — + 27 



NcOCg f H \ 
~3n 



(4717) ( 



IniD 



2€ 



mlT. 
(5.4.15) 



(ss) contribution 

The (ss) contribution from the two-loop diagrams with a single self -energy insertion 
can be easily obtained by expanding the two-loop result in powers of m^. This yields 



3h-31 



pq 



2{p^+4ml) 





NcUs / 


e 


3n \ 



u 

^ ^ ml)T^ . (5.4.16) 



4 

3 
16 



We have verified this by explicitly calculating the relevant diagrams. 
(hhh) contribution 

If all the three loop momenta are hard, we can obtain the nio/T expansion simply by 
expanding in powers of m|,. To obtain the expansion through order g^, we can use bare 
propagators and vertices. The contributions from the three-loop diagrams were first 
calculated by Arnold and Zhai in Ref. ||5|, and later by Braaten and Nieto |[12|. One 
finds 
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(?2?2s) contributions 

All the diagrams except (3f) are infrared finite in the limit — > 0. This implies that the 
contribution is given by using a dressed longitudinal propagator and bare vertices. 
The ring diagram (3f ) is infrared divergent in that limit. The contribution through is 
obtained by expanding in powers of self -energies and vertices and one obtains 

+ |/(.-in.^-n. + 46]/^^^^^^^ 
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(hss) contribution 



For all the diagrams that are infrared safe, the (hss) contribution is of order g'^rn^, i.e. 
g^ and can be ignored. The infrared divergent diagrams contribute as follows 
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(sss) contribution 

The (sss) contribution is given by a three-loop calculation of the free energy of EQCD 
in three dimensions. This calculation was performed in Ref. IIT2II . Alternatively, one can 
isolate the (sss) contributions from the diagrams listed in Ref. HSj. The result is 

1 1 

4 (p2 + ml){q^ + m^,) (r2 + ml)^ 
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The expression for the integrals are given in Appendix|Cj Adding Eqs. ( C.3.1| - ( C.3.13| , 
the final result is 
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Note that all the poles in e cancel. 



5.5 Thermodynamic potentials 

In this section we present the final renormalized thermodynamic potential explicitly 
through order 5^, aka NNLO. The final NNLO expression is completely analytic; 
however, there are some numerically determined constants which remain in the final 
expressions at NLO. 



5.5.1 Leading order 

The leading order thermodynamic potential is given by the contribution from the 
diagrams (la) and (lb) 



Hi-ioop = -^ideal |l - ym2, + 30m|, + ^ + 21og I - 7 + 27 + , 

(5.5.1) 



where Jideai is the free energy of a gas of — 1 massless spin-one bosons and mo and 
p. are dimensionless variables: 

-^ideal = [Nl-l) (-^^') > (5-5.2) 

= (5.5.3) 
H = ■ (5.5.4) 

The complete expression for the leading order thermodynamic potential is given by 
adding the leading vacuum energy counterterm ( |5.3.3| to Eq. (5.5.1 1: 



Olo = -^ideai <! 1 - y m?, + 30^3, + y ("log |-^ + 7 + y)m|,[>, (5.5.5) 



5.5.2 Next-to-leading order 

The renormalization contributions at first order in S are 

AQi = AK^ + Aimlj-^QLo . (5.5.6) 
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Using the results listed in Eqs. (5.3.5 1 and (5.3.61 the complete contribution from the 
counterterm at first order in 3 is 



AQl = J^ideal i I^WJd + 



NcO-s .2 
OK 



495 



^+21og|-21og?«D + 2 







(5.5.7) 



Adding the NLO counterterms ( 5.5. 7[ | to the contributions from the various NLO 
diagrams, we obtain the renormalized NLO thermodynamic potential ||2T]| 
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5.5.3 Next-to-next-to-leading order 

The renormalization contributions at second order in 5 are 
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Using the results listed above, we obtain 
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Adding the NNLO counterterms ( 5.5.10| l to the contributions from the various NNLO 
diagrams we obtain the renormalized NNLO thermodynamic potential. We note that 
at NNLO all numerically determined coefficients of order drop out and we are left 
with a final result which is completely analytic. The resulting NNLO thermodynamic 
potential is 
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We note that the coupling constant counterterm listed in Eq. ( 5.2.5| coincides with 
the known one-loop running of the QCD coupling constant 



llNc^ 

24712 



(5.5.12) 



Finally, note that if we use the weak-coupling value for the Debye mass = 
4i7iNciXsT-^/3, the NNLO HTLpt result ( 5.5.1l| l reduces to the weak-coupling result 
through order and we have checked that this is the case. 



5.6 Thermodynamic functions 
5.6.1 Mass prescriptions 

The mass parameter niu m HTLpt is completely arbitrary. To complete a calculation, it 
is necessary to specify niu as function of g and T. In our case this implies that the free 
energy T is obtained by specifying mo as a function of T and in the thermodynamic 
potential Q. In this section we will discuss several prescriptions for the mass parameter. 



Variational Debye mass 

The variational mass is given by the solution to the variational gap equation which is 
defined by 



dm 
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(5.6.1) 



Applying it to ( 5.5. 11| , the NNLO gap equation reads 
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At leading order in HTLpt, the only solution is the trivial solution, i.e. = 0. In 
that case, it is natural to chose the weak-coupling result for the Debye mass. This was 
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done in Ref. IIZOl . At NLO, the resulting gap equation has a nontrivial solution, which 
is real for all values of the coupling ||2T1 . At NNLO, the solution to the gap equation is 
plagued by imaginary parts for all values of the coupling. The problem with complex 
solutions seems to be generic since it has also been observed in screened perturbation 
theory II30I and QED II22I . In those cases, however, it was complex only for small values 
of the coupling. 

Perturbative Debye mass 

At leading order in the coupling constant g, the Debye mass is given by the static 
longitudinal gluon self-energy at zero three-momentum, = Yli{0,0), i.e. 

^NcOCsT^ . (5.6.3) 

The next-to-leading order correction to the Debye mass is determined by resummation 
of one-loop diagrams with dressed vertices. Furthermore, since it suffices to take into 
account static modes in the loops, the HTL-corrections to the vertices also vanish. The 
result, however, turns out to be logarithmically infrared divergent, which reflects the 
sensitivity to the nonperturbative magnetic mass scale. The result was first obtained by 
Rebhan EH and reads * 

3ml = ml^^l^a^'^ 

where nJmag is the nonperturbative magnetic mass. We will not use this mass prescrip- 
tion since it involves the magnetic mass which would require input from e.g. lattice 
simulations. 

BN mass parameter m\ 

In the previous subsection, we saw that the Debye mass is sensitive to the nonper- 
turbative magnetic mass which is of order g^T. In QED, the situation is much better. 
The Debye mass can be calculated order by order in e using resummed perturbation 
theory. The Debye mass then receives contributions from the scale T and eT. Effective 
field theory methods and dimensional reduction can be conveniently used to calculate 
separately the contributions to from the two scales in the problem. The contribu- 
tions to and other physical quantities from the scale T can be calculated using bare 

*In Ref. ItTI , it was shown that the gauge dependent part of the static gluon self-energy 11^^(0, k) 
vanishes when it is evaluated on shell, i.e. when = —nil. ^^^^ is in accordance with general gauge- 
dependece identities l44l . 
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(5.6.4) 
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propagators and vertices. The contributions from the soft scale can be calculated using 
an effective three-dimensional field theory called electrostatic QED. The parameters 
of this effective theory are obtained by a matching procedure and encode the physics 
from the scale T. The effective field theory contains a massive field Aq that up to 
normalization can be identified with the zeroth component of the gauge field in QED. 
The mass parameter of Aq gives the contribution to the Debye mass from the hard 
scale T and can be written as a power series in e^. For non-Abelian gauge theories, the 
corresponding effective three-dimensional theory is called electrostatic QCD. The mass 
parameter niE for the field Aq (which lives in the adjoint representation) can also be 
calculated as a power series in It has been determined to order by Braaten and 
Nieto II12I . For pure-glue QCD, it reads 



l + ^U + ^T + ^log 



3n V4 2 ' 2 °A7tT 



(5.6.5) 



We will use the mass parameter wje as another prescription for the Debye mass and 
denote it by the Braaten Nieto (BN) mass prescription. 

5.6.2 Pressure 

In this subsection, we present our results for the pressure using the variational mass 
prescription and the BN mass prescription. 

Variational mass 



In Fig. 5.5 we compare the LO, NLO, and NNLO predictions for the real part of the 
pressure normalized to that of an ideal gas using the variational mass and three-loop 
rurming of ocg llTOll . Shaded bands show the result of varying the renormalization scale 
^ by a factor of 2 around }i = 2nT. 



In Fig. 5.6 we show the NNLO result for the imaginary part of the pressure 
normalized by the ideal gas pressure using the variational mass and three-loop ruiming 
of DCs II70I . The imaginary part decreases with increasing temperature and is rather 
small beyond 3 — 4 T^. 

Due to the imaginary parts, we abandon the variational prescription in the remain- 
der of the chapter. 

BN mass 



In Fig. 5.7 we show the HTLpt predictions for the pressure normalized to that of an 



ideal gas using the BN mass prescription and one-loop rurming of Oig in Eq. ( 5.5.12[ | 



The bands are obtained by varying the renormalization scale by a factor of 2 around 
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Figure 5.5: Comparison of LO, NLO, and NNLO predictions for the scaled real part of the 
pressure using the variational mass and three-loop running iTOll . Shaded bands show the 
result of varying the renormalization scale ^ by a factor of 2 around }i = InT. 



}i = InT. In Fig. 5.8 we again plot the normalized pressure, but now with three-loop 
running of <Xs llTOl . The agreement between the lattice data from Boyd et al. 1541 is very 
good down to temperatures of around 3 Tc. Comparing Figs. 5.7 -{5^ we see that using 
the three-loop rurming, the band becomes wider. However, the difference is significant 
only for low T, where the HTLpt results disagrees with the lattice anyway. For T > 3Tc, 
the prescription for the running makes very little difference. 

Until recently, lattice data for thermodynamic variables only existed for tempera- 
tures up to approximately 5 Tc. In the paper by Endrodi et al II72II , the authors calculate 
the pressure on the lattice for pure-glue QCD at very large temperatures. In Fig. |5.9[ we 
show the results of Endrodi et al as well as Boyd et al, together with the HTLpt NLO 
and NNLO predictions for the pressure using the BN mass prescription and three-loop 
running of as II70I . The two points from Ref. I|72ll have large error bars, but data points 
are consistent with the HTLpt predictions. 

It is interesting to make a comparison of convergence between HTLpt and weak- 
coupling expansion. Analyzing the weak-coupling result listed in Eq. (2.3.351, we 
find for the case of SU(3) Yang-Mills theory that in order to make the magnitude 
of the coefficients of each order smaller than that of the previous order, one has to 
require the temperature be higher than 5.36 x 10^ Tc. However HTLpt meets the same 
requirement at temperatures higher than 25.6 Tc, which is an improvement of five 
orders of magnitude. 
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Figure 5.6: The NNLO result for the scaled imaginary part of the pressure using the 
variational mass and three-loop running llTOll . The two curves arises from the two complex 
conjugate solutions to the gap equations. Shaded bands show the result of varying the 
renormalization scale f/ by a factor of 2 around }i = 2nT. 



5.6.3 Pressure at large Nc 

The large Nc limit is achieved by taking Nc to be large while holding g^Nc fixed. The 
large Nc coupling is defined by A = g^Nc- As can be seen from Eqs. (5.5.51, (5.5.81 



and ( 5.5.11[ l that the ratios of the thermodynamic potentials over Jideal solely 



functions of A which have no residual dependence on Nc or g after the substitution 
g^ — > A/Nc, while the same is true for the BN mass (5.6.51. Therefore the scaled 



HTLpt thermodynamics, i.e. the thermodynamics obtained by taking ratio over J-ideai/ 
is independent of the actual number of colors Nc up to three-loop order. This is in 
line with a recent lattice study by Panero II73I showing that the thermodynamics of 
SU(N) Yang-Mills theories has a very mild dependence on Nc, supporting the idea 
that the QCD plasma could be described by models based on the large Nc limit. In 
Fig. ( 5.10[ l we plot the HTLpt predictions for the pressure at large Nc through NNLO 



with three-loop running llTOll at f/ = 2tzT together with the SU(3) prue-glue lattice data 
from Boyd et al. ||54| . The curves in Fig. ( 5.10| are exactly the same as those in Fig. ( |5.8| 



demonstrating the independence of the scaled SU (Nc) pressure with respect to Nc. It 
is unknown at this stage whether higher-order HTLpt contributions would spoil this 
independence, however from the comparison of the NNLO result with the lattice data, 
the Nc dependence from the higher order corrections for T > 3Tc might be tiny in case 
there are any. 
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Figure 5.7: Comparison of LO, NLO, and NNLO predictions for the scaled pressure using 
the BN mass and one-loop running ( 5.5.12| . The points are lattice data for pure-glue with 
Nc = 3 from Boyd et al. 1541 . Shaded bands show the result of varying the renormalization 
scale by a factor of 2 around = 27rT. 



5.6.4 Energy density 

The energy density S is defined by 



£ = T-T^. (5.6.6) 



In Fig. 5.11 we show the LO, NLO, and NNLO predictions for energy density normal- 
ized to that of an ideal gas using the BN mass prescription and three-loop running of 
ft^s llTOll . The bands show the result of varying the renormalization scale by a factor of 
2 around }i = 2nT. Our NNLO predictions are in very good agreement with the lattice 
data down to T ~ 2 T^. 



5.6.5 Entropy 

The entropy density is defined by 



S = . (5.6.7, 



In Fig. 5.12 we show the entropy density normalized to that of an ideal gas using the 
BN mass prescription and three-loop running of ocs II70I . The points are lattice data 
from Boyd et al. 1541 . Our NNLO predictions are in excellent agreement with the lattice 
data down to T ~ 2 Tc. 
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Figure 5.8: Comparison of LO, NLO, and NNLO predictions for the scaled pressure using 
the BN mass and three-loop running llTOl with SU(3) pure-glue lattice data from Boyd et al. 
Il54ll . Shaded bands show the result of varying the renormalization scale }i by a factor of 2 
around = 2nT. 



5.6.6 Trace anomaly 

In pure-glue QCD or in QCD with massless quarks, there is no mass scale in the 
Lagrangian and the theory is scale invariant. At the classical level, this implies that the 
trace of the energy-momentum tensor vanishes. At the quantum level, scale invariance 
is broken by renormalization effects. It is convenient to introduce the scale anomaly 
density £ — 3V, which is proportional to the trace of the energy-momentum tensor. 
The trace anomaly can be written as 



In Fig. 5.13 we show the HTLpt predictions for the trace anomaly divided by using 
the BN mass prescription and three-loop running of ocs II70I . The points are lattice data 
from Boyd et al. Il54l . For temperatures below approximately 2 Tc, there is a large 
discrepancy between the HTLpt predictions and the lattice. At LO and NLO, the curves 
are even bending downwards. 

At temperatures close to the phase transition it has been suggested that the dis- 
crepancy between HTLpt resummed predictions for thermod5mamics functions and, in 
particular, the trace anomaly is due to influence of a dimension two condensate Il74ti76ll 
which is related to confinement. Phenomenological fits of lattice data which include 
such a condensate show that the agreement with lattice data is improved [ZZllZHl- Al- 
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Figure 5.9: Comparison of NLO, and NNLO predictions for the scaled pressure using the 
BN mass and three-loop running llTO] with SU(3) pure-glue lattice data from Boyd et al. 
15411 and Endrodi et al. Il72l . Shaded bands show the result of varying the renormalization 
scale jihy a factor of 2 around }i = 2nT. 




Figure 5.10: Comparison of LO, NLO, and NNLO predictions for the scaled pressure at 
large Nc using the BN mass and three-loop running llTOl with SU(3) pure-glue lattice data 
from Boyd et al. ||54| . }i = 2nT is taken here. 



tematively, others have constructed AdS/ CFT inspired models which break conformal 
in variance "by hand" Il79fl8l1l . These models are also able to fit the thermodynamical 
functions of QCD at temperatures close to the phase transition. 



In Fig. 5.14 we show the HTLpt predictions for the trace anomaly scaled by T^/T^ 
using the BN mass prescription and three-loop rurming of as llTOl . The points are 
lattice data from Boyd et al. ||54| . The most remarkable feature is that lattice data are 
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Figure 5.11: Comparison of LO, NLO, and NNLO predictions for the scaled energy density 
using the BN mass and three-loop rurming ||70] with SU(3) pure-glue lattice data from 
Boyd et al. 1541 . Shaded bands show the result of varying the renormalization scale jihy a 
factor of 2 around }i = 2tzT. 



essentially constant over a very large temperature range. Clearly, HTLpt does not 
reproduce the scaled lattice data precisely; however, the agreement is dramatically 
improved when going from NLO to NNLO. 



5.7 Conclusions 



In this chapter, we have presented results for the LO, NLO, and NNLO thermodynamic 
functions for SU(Nc) Yang-Mills theory using HTLpt. We compared our predictions 
with lattice data for Nc = 3 and found that HTLpt is consistent with available lattice 
data down to approximately T ~ 3 Tc in the case of the pressure and T ~ 2 Tc in the 
case of the energy density and entropy. These results are in line with expectations since 
below T ~ 2 — 3 Tc a simple "electric" quasiparticle approximation breaks down due 
to nonperturbative chromomagnetic effects ll5Tll52ll ^. This is a nontrivial result since, 
in this temperature regime the QCD coupling constant is neither infinitesimally weak 
nor infinitely strong with ^ ~ 2, or equivalently ocs = g^/ (47r) ~ 0.3. Therefore, we 
have a crucial test of the quasiparticle picture in the intermediate coupling regime. 

The mass parameter ?md in HTLpt is arbitrary and we employed two different 
prescriptions for fixing it. Unfortunately, the variational gap equation has four complex 

'f There have been also hints that the Z(N) interface l82l , gauge-fixing ambiguities l83l , and topological 
objects such as quantum instantons |84] and magnetic monopoles l85l might play important roles on the 
thermodynamics at intermediate temperature. 
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Figure 5.12: Comparison of LO, NLO, and NNLO predictions for the scaled entropy using 
the BN mass and three-loop running ||70l with SU(3) pure-glue lattice data from Boyd et al. 
Il54ll . Shaded bands show the result of varying the renormalization scale }i by a factor of 2 
around }i = 2nT. 

conjugate solutions, two with positive real parts. This has also been observed in scalar 
theory and QED. Whether this is a problem of HTLpt as such or is related to our nio/T 
expansion is unknown. Since it is not currently possible to evaluate the NNLO HTLpt 
diagrams in gauge theories exactly, it is impossible to settle the issue at this stage. On 
the other hand, the BN mass prescription is well defined to all orders in perturbation 
theory and does a reasonable job reproducing available lattice data for temperatures 
above T>3Tc. With QED and Yang-Mills results at hand, the NNLO full QCD HTLpt 
thermodynamics will be a routine extension Il86ll87ll . 
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Figure 5.13: Comparison of LO, NLO, and NNLO predictions for the scaled trace anomaly 
using the BN mass and three-loop running llTOl with SU(3) pure-glue lattice data from 
Boyd et al. 1541 . Shaded bands show the result of varying the renormalization scale ^ by a 
factor of 2 around ji = 2nT. 




Figure 5.14: Comparison of LO, NLO, and NNLO predictions for the scaled trace anomaly 
using the BN mass and three-loop running llTOl with SU(3) pure-glue lattice data from 
Boyd et al. 1541 . Shaded bands show the result of varying the renormalization scale }ihy a 
factor of 2 around // = 2nT. 
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Chapter 6 

Summary and Outlook 



This dissertation is devoted to the study of thermodynamics for thermal gauge theories. 
The poor convergence of conventional perturbation theory has been the main obstacle 
in the practical application of thermal QCD for decades. To improve this embarrassing 
situation, a considerable effort has been put into reorganizing the perturbative series 
at phenomenologically relevant temperatures. The application of hard- thermal-loop 
perturbation theory to the problem carried out in this dissertation leads to laudable 
results for both Abelian and non-Abelian theories. 

The success of HTLpt is not totally unexpected since it is essentially just a reorgani- 
zation of perturbation theory which shifts the expansion from around an ideal gas of 
massless particles to that of massive quasiparticles which are the real degrees of free- 
dom at high temperature. The HTL Feynman rules listed in Appendix [A| show clearly 
that the propagators and vertices are dressed systematically by the thermal medium, 
as a result the interactions get screened in the medium which can be seen, for instance, 
from the quark-gluon three-vertex ( A.8.1[ l that the coupling strength gets screened by 
the thermal mass term explicitly. Therefore the expansion in terms of the HTL Feynman 
rules are self-consistently around a gas of thermal quasiparticles. The fact that the 
mass parameter is not arbitrary but a function of g and T determined variationally 
or perturbatively also indicates that HTLpt doesn't modify the original gauge theory 
but just reorganizes its perturbation series. Gauge invariance which is guaranteed 
by construction in HTLpt is useful both as a consistency check in calculations and as 
a way to simplify calculations. Although the renormalizability of HTLpt is not yet 
proven, the fact that it is renormalizable at NNLO using only known counterterms 
shows promising light along the way. 

So far, thermodynamics for quantum fields has been studied intensively in the 
commimity, both perturbatively through higher orders or numerically on the lattice, 
however real-time dynamics is still in its very early stage of development. Transport 
coefficients are of great interest since they are theoretically clean and well defined non- 
equilibrium dynamical quantities. Along the line of perturbative approach to transport 
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coefficients, considerable efforts have been devoted at leading order in weak-coupling 
expansion II88I . However the only known transport coefficients to next-to-leading 
order are shear viscosity in scalar cp'^ theory Il89l , heavy quark diffusion rate in QCD 
and Af = A supersymmetric Yang-Mills theory II90I , and transverse diffusion rate q in 
QCD lU, and all of them exhibit poor convergence as bad as the case of thermodynamic 
quantities, such as the pressure. Since djmamical quantities are still not well described 
by lattice gauge theory, new resummation techniques are urgently needed in order to 
achieve a better understanding of transport coefficients. 

Although the papers written to date have focussed on using HTLpt to compute 
thermodynamic observables, the goal of this work is to create a framework which can 
be applied to both equilibrium and non-equilibrium systems. HTLpt is formulated in 
Minkowski space, so its application to non-equilibrium dynamics is straightforward. 
With the confidence from thermodynamics, HTLpt is ready to enter the domain of 
real-time dynamics and this might be of great help in deepening our knowledge in the 
properties of the quark-gluon plasma. 

abs.tex include ack.tex include intro.tex include sum.tex include 
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Appendix A 

HTL Feynman Rules 



In this appendix, we present Feynman rules for HTL perturbation theory in QCD, from 
which QED Feynman rules can be obtained by simplifying the relevant color structures. 
We give explicit expressions for the propagators and for the quark-gluon three- and 
four-vertices. The Feynman rules are given in Minkowski space to facilitate future 
applications to real-time processes. A Minkowski momentum is denoted p = {po, p), 
and the irmer product isp ■ q = poqo — p • q. The vector that specifies the thermal rest 
frame is n = (1,0). 

A.l Gluon self-energy 

The HTL gluon self -energy tensor for a gluon of momentum p is 



The tensor T^"^{p, q), which is defined only for momenta that satisfy p + q = Q,is 



The angular brackets indicate averaging over the spatial directions of the light-like 
vector y = (l,y). The tensor T^^ is symmetric in }i and v and satisfies the "Ward 
identity" 



The self -energy tensor n'"' is therefore also symmetric in y. and v and satisfies 



W^p) 



ml[T^\p,-p)-n^n'] . 



(A.1.1) 




(A.1.2) 



(A.1.3) 



Pf.n^'ip) = 0, 

g^yW'ip) = -ml . 



(A.1.4) 
(A.1.5) 
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The gluon self -energy tensor can be expressed in terms of two scalar functions, the 
transverse and longitudinal self -energies and TLi, defined by 



(A.1.6) 
(A.1.7) 



where p is the unit vector in the direction of p. In terms of these functions, the self- 
energy tensor is 



(A.1.8) 



where the tensors Tp and Lp are 



The four-vector rip is 







p^p" 


n^pH'p 


yllV _ 

— 


npn^p 
n2 







(A.1.9) 
(A.1.10) 



(A.1.11) 



and satisfies p-tip = and nt, = 1 — {n-pY /p^. ( A.1.5 1 reduces to the identity 



(d-l)nr(p) + -yn,(p) = ml. 



(A.1.12) 



We can express both self-energy functions in terms of the function defined by 
( |AT2| : 



nr(p) 



m 



D 



r\p,-p)-l + nl 



{d-l)nl 
Udp) = ml[l-T'\p,-p)] , 



(A.1.13) 
(A.1.14) 



In the tensor T^'"{p,—p) defined in (A.1.2i, the angular brackets indicate the 
angular average over the unit vector y. In almost all previous work, the angular 
average in ( A.1.2| has been taken in d = 3 dimensions. For consistency of higher order 
corrections, it is essential to take the angular average in d = 3 — 2e dimensions and 
analytically continue to d = 3 only after all poles in e have been cancelled. Expressing 
the angular average as an integral over the cosine of an angle, the expression for the 00 
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component of the tensor is 

rOO(p, -p) = ^ (1 - c')-'^^ , (A.1.15) 



where the weight function w[e) is 



^ ^ r2(i-e) r(|)r(i-e) ^ ' 



The integral in ( A.1.15[ l must be defined so that it is analytic at po = It then has a 



branch cut rurming from po = — I p I to po = + 1 p | . If we take the limit e — ?► 0, it reduces 
to 

r-(p,-p) = ^log^, (A.1.17) 
which is the expression that appears in the usual HTL self-energy functions. 

A,2 Gluon propagator 

The Feynman rule for the gluon propagator is 

i5'''^^y{f) , (A.2.1) 

where the gluon propagator tensor A^,y depends on the choice of gauge fixing. We 
consider two possibilities that introduce an arbitrary gauge parameter ^: general 
covariant gauge and general Coulomb gauge. In both cases, the inverse propagator 
reduces in the limit ^ — ?► oo to 

A-^(p)f"' = + V^'f - n^'(p) • (A.2.2) 

This can also be written 

^-^{^yy - + \ l''^ (A 2 3) 

" At(p) ^n2AL(p)^P ' ^^-^-^^ 

where and Al are the transverse and longitudinal propagators: 
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The inverse propagator for general ^ is 



covariant , (A.2.6) 



AJipY'" -^{pf -p-nn^'){p'' -p-nn'') Coulomb. (A.2.7) 



The propagators obtained by inverting the tensors in ( A.2.7| and ( A.2.6| l are 



A^^ip) = -AT{p)t;' + A,{p)nX-^\^2)2 
= -Ar(p)T;^ + Ai.(p)nfV'-^ P'^' 



covariant , (A.2.8) 
Coulomb . (A.2.9) 



It is convenient to define the following combination of propagators: 



Ax(p) = Al{p) + -^Arip) . 



(A.2.10) 



Using ( A.1.12| ^ (A.2.41, and (A.2.51, it can be expressed in the alternative form 



Ax(p) = [m|,-dnr(p)] Al(p)At(p), 



(A.2.11) 



which shows that it vanishes in the limit mo — t- 0. In the covariant gauge, the propaga- 
tor tensor can be written 



Af-(p) 



-At(p)/'^ + Ax(p)nf'n^] - ^Ax(p) {p^'n' + n^p') 



+ 



Mp) + ^Ax(p)-^ 



pt^p 



(A.2.12) 



This decomposition of the propagator into three terms has proved to be particularly 
convenient for explicit calculations. For example, the first term satisfies the identity 

[-Ar{p)g,. + A^{p)n,n,]A-J{pr = - ^ + ^^p,nj . (A.2.13) 



A.3 Three-gluon vertex 



The three-gluon vertex for gluons with outgoing momenta p, cj, and r, Lorentz indices 
ji, V, and A, and color indices a, b, and c is 



(A.3.1) 
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where f'^'^ are the structure constants and the three-gluon vertex tensor is 

(A.3.2) 

The tensor T^^^ in the HTL correction term is defined only for p + + r = 0: 

T^^\p,q,r) = -iffy'' (-^ \ . (A.3.3) 

\ VP-yj-y r-yq-yj / 

This tensor is totally symmetric in its three indices and traceless in any pair of indices: 
gpwT^^^ = 0. It is odd (even) under odd (even) permutations of the momenta p, q, and 
r. It satisfies the "Ward identity" 

qnT^'^Prqrr) = T'\p + q,r) -T^'^pj + q) . (A.3.4) 

The three-gluon vertex tensor therefore satisfies the Ward identity 

p^I^'\p,q,r) = A-\qy^ - A^\ry^ . (A.3.5) 

A.4 Four-gluon vertex 

The four-gluon vertex for gluons with outgoing momenta p, q, r, and s, Lorentz indices 
}i, V, A, and a, and color indices a, b, c, and d is 

KliliP'^'r^') = -^s'{fa>,.f.cd{g^'Y^-g"Y') 

+ Imltr (t^T'T'' + T'^rr^^^ T^^'"'{p,q,r,s)] 

+ 2 cyclic permutations , (A.4.1) 

where the cyclic permutations are of {q, v, h), {r, A, c), and (s, a, d). The matrices are 
in the fundamental representation of the SU (Nc) algebra with the standard normaliza- 
tion tr(T«T^) = IS"^. The tensor T^"'^'' in the HTL correction term is defined only for 
p + q + r + s = 0: 

r^'^''{p,q,r,s) = ly^fy^y^'i '^''^ , , 

ip + q)-n , (p + q + r)-n \ \ ^^^^^ 



q-yr-y {r + s)-y r-y s-y {s + p)-y 

This tensor is totally symmetric in its four indices and traceless in any pair of indices: 
g}ivT^^^ = 0. It is even under cyclic or anti-cyclic permutations of the momenta p, q, r, 
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and s. It satisfies the "Ward identity" 

qfiV"'^''{p,q,r,s) = T^-{p + q,r,s) - r'^-{p,r + q,s) (A.4.3) 
and the "Bianchi identity" 

rf''^^p,q,r,s) + r''''"'{p,r,s,q) + V"'^^{p,s,q,r) = 0. (A.4.4) 

When its color indices are traced in pairs, the four-gluon vertex becomes particu- 
larly simple: 

^"'^""KlcIiP'I'r^s) = -ig^NM-l)V^''^''{p,q,r,s), (A.4.5) 
where the color-traced four-gluon vertex tensor is 

Tf''^''{p,q,r,s) = Igl^-g^- - f^g-- - gl^-g-^ - mlT^'^-{p,s,q,r) . (A.4.6) 

Note the ordering of the momenta in the arguments of the tensor 'J'H^'^'^ ^ which comes 
from the use of the Bianchi identity ( A.4.4| . The tensor (A.4.6 1 is symmetric under the 



interchange of }i and v, under the interchange of A and a, and under the interchange 
of v) and (A, a). It is also symmetric under the interchange of p and q, under the 
interchange of r and s, and under the interchange of {p,q) and {r,s). It satisfies the 
Ward identity 

p^r'"''^'^(p,^?,r,s) = r^-(^,r + p,s)-r^'^(^j,r,s + p). (A.4.7) 



A.5 HTL gluon counterterm 

The Feynman rule for the insertion of an HTL counterterm into a gluon propagator is 

-i5'^U^^{p), (A.5.1) 
where n''^(p) is the HTL gluon self-energy tensor given in ( A.1.8| |. 



A.6 Quark self-energy 

The HTL self-energy of a quark with momentum p is given by 



E(P) = m]np) , 
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where 

r-V) . . (A.6.2, 

Expressing the angular average as an integral over the cosine of an angle, it becomes 



Po - \p\c 



(A.6.3) 



The integral in ( A.6.3[ l must be defined so that it is analytic at po = It then has a 



branch cut running from po = — I P I to po = + 1 p | • In three dimensions, this reduces to 

2 

m 



UP) = ;^T70log^^^ + ;^7-p{l-;sT^loge^) • (A.6.4) 



Po + 


IpI 


Po- 


IpI 



A.7 Quark propagator 

The Feynman rule for the quark propagator is 

iS'^Sip) . (A.7.1) 
The quark propagator can be written as 



where the quark self-energy is given by ( A.6.1[ |. The inverse quark propagator can be 
written as 

S-i(p) = p-^p). (A.7.3) 

This can be written as 

S-\p) =A{p), (A.7.4) 
where we have organized Ao{p) and As{p) into: 

A^ip) = (Ao(p),As(p)p). (A.7.5) 
The functions Ao{p) and As{p) are defined as 

Mp) = Vo-yTp, (A.7.6) 
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2 

m 

Mp) = lpl + ^[l-rp] . (A.7.7) 

A.8 Quark-gluon three-vertex 

The quark-gluon three-vertex with outgoing gluon momentum p, incoming fermion 
momentum q, and outgoing quark momentum r, Lorentz index ^ and color index a is 

n{p,(1,r) = gt,{7^-m^^n{p,q,r)) . (A.8.1) 

The tensor in the HTL correction term is only defined for p — q + r = 0: 

rnp,q,r) = (y^(—^]) . (A.8.2) 



This tensor is even under the permutation of q and r. It satisfies the "Ward identity" 

p^f^ {p, q, r) = fiq)- (r) . (A.8.3) 
The quark-gluon three-vertex therefore satisfies the Ward identity 

p^I^{p,q,r) = S-\q) - S~\r) . (A.8.4) 

A.9 Quark-gluon four-vertex 

We define the quark-gluon four-point vertex with outgoing gluon momenta p and q, 
incoming fermion momentum r, and outgoing fermion momentum s. Generally this 
vertex has both adjoint and fundamental indices, however, for this calculation we will 
only need the quark-gluon four-point vertex traced over the adjoint color indices. In 
this case 

^VijiP'q'r,s) = -gW^Ci:Sijf^%p,q,r,s) 

= g'cFSijT^^ (A.9.1) 

where Cf = (N^ — l)/(2Nc). There is no tree-level term. The tensor in the HTL 
correction term is only defined forp + q — r + s = 

This tensor is symmetric in ^ and v and is traceless. It satisfies the Ward identity: 

Pj,T^''{p,q,r,s) = r{q,r-p,s)-r{q,r,s + p) . (A.9.3) 
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A.IO HTL quark counterterm 

The Feynman rule for the insertion of an HTL quark counterterm into a quark propaga- 
tor is 

iS^^-Lip) , (A.10.1) 
where is the HTL quark self-energy given in ( A.6.1[ |. 



A.ll Ghost propagator and vertex 

The ghost propagator and the ghost-gluon vertex depend on the gauge. The Fejmman 
rule for the ghost propagator is 

covariant , ( A. 1 1 . 1 ) 

Coulomb . (A. 11. 2) 

The Feynman rule for the vertex in which a gluon with indices }i and a interacts with 
an outgoing ghost with outgoing momentum r and color index c is 

-gf'^'r^ covariant , (A.11.3) 

-gf''' ij^' - ^-n n'') Coulomb . (A.11.4) 

Every closed ghost loop requires a multiplicative factor of —1. 



I gab 



A.12 Imaginary-time formalism 

In the imaginary-time formalism, Minkoswski energies have discrete imaginary val- 
ues po = i2nnT for bosons and po = ^ (2m + l)7rr for fermions, and integrals over 
Minkowski space are replaced by sum-integrals over Euclidean vectors {2nnT, p) or 
((2n + l)7rT,p), respectively. We will use the notation P = {Po,p) for Euclidean mo- 
menta. The magnitude of the spatial momentum will be denoted p = |p|, and should 
not be confused with a Minkowski vector. The inner product of two Euclidean vec- 
tors is P ■ Q = PqQo + P ■ q- The vector that specifies the thermal rest frame remains 
n = (1,0). 

The Feynman rules for Minkowski space given above can be easily adapted to 
Euclidean space. The Euclidean tensor in a given Feynman rule is obtained from the 
corresponding Minkowski tensor with raised indices by replacing each Minkowski 
energy po by zPq, where Pq is the corresponding Euclidean energy, and multipying by 
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—i for every index. This prescription transforms p = (po/ p) into P = {Pq, p), gf^^ into 
-(5''", and p-q into -P-Q. The effect on the HTL tensors defined in ( |A.1.2| , ( |A.8.2| , and 
(A.9.21 is equivalent to substituting p-n —> —P-N where N = {—i,0), p-y —> —P-Y 



where Y = {—i,y), and yf — ?► Y''. For example, the Euclidean tensor corresponding to 
( |AT2| is 



V{P,-P) 



ypiyv 



P-N 

Yy 



(A.12.1) 



The average is taken over the directions of the unit vector y. 

Alternatively, one can calculate a diagram by using the Feynman rules for Minkowski 
momenta, reducing the expressions for diagrams to scalars, and then make the appro- 
priate substitutions, such as — > —P'^, p ■ q ^ —P Q, and n ■ p ^ in ■ P. For example, 
the propagator functions ( A.2.4| and (A.2.5 1 become 



At(P) = 
Al(P) = 



P^ + UriP) ' 
1 



(A.12.2) 
(A.12.3) 



p2 + ni.(p) • 

The expressions for the HTL self-energy functions Y\j{P) and Wi^{P^ are given by 



(A.1.131 and (A.1.14 1 with Hp replaced by n| = p^ /P'^ and T^'^ip, — p) replaced by 



Tp 



2 J-1 ^ ' iPo-pc 



(A.12.4) 



Note that this function differs by a sign from the 00 component of the Euclidean tensor 
corresponding to ( A.1.2| : 



p ■ 



(A.12.5) 



T°0(P,-P) = -T^\p,-p) 
A more convenient form for calculating sum-integrals that involve the function Tp is 

Tp = ( ^\. ) , (A.12.6) 



where the angular brackets represent an average over c defined by 



(/(c)), ^ w{e) / dc(l-c^)-7(c) 

J 



(A.12.7) 



and w{e) is given in ( A.1.16| l. 
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In the imaginary-time formalism for thermal field theory, the four-momentum P = 
(Pq, p) is Euclidean with = Pq + p^. The Euclidean energy Pq has discrete val- 
ues: Pq = 2nnT for bosons and Pq = (2n + l)7rT for fermions, where n is an integer. 
Loop diagrams involve sums over Pq and integrals over p. With dimensional regu- 
larization, the integral is generalized to d = 3 — 2e spatial dimensions. We define the 
dimensionally regularized sum-integrals by 




bosons , (B.0.1) 
fermions , (B.0.2) 



where 3 — 2e is the dimension of space and }i is an arbitrary momentum scale. The factor 
{e^ / AlTcY is introduced so that, after minimal subtraction of the poles in e due to ultravi- 
olet divergences, }i coincides with the renormalization scale of the MS renormalization 
scheme. 



B.l One-loop sum-integrals 

The simple one-loop sum-integrals required in our calculations can be derived from 
the formulas 

4 ^ ( H ^^e 2Til + m-e)T{n-l-m + e) 

fp(p2)« Utzt) r(n)r(2-2e) ^ ^ 

X C{2n -2m -3 + 26)T^+^"'-^" {271^+^"'-^" , (B.1.1) 
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The specific bosonic one-loop sum-integrals needed are 



45 



4 ^ - ^ 

TpP^ ~ VI V47rT 
1 



2e 



V_ 

(47r)2 V47rT 



27 + 0(e) 
27 + 2 + 0(e) 



(47r)2 V47rT 

The specific fermionic one-loop sum-integrals needed are 



i logp' 

T[P} 

i 



360 



24 V47rT 



le 



1+ 2-21og2 + 2 



2e 



(47r)2 V47rr 

T2 / \2e 



+ 27 + 41og2 + 0(e) 



{P} C^'^)^ 



16 V47rT 
3 1 



l+(^-21og2 + 2||-^)e + 0(e2) 



(B.1.3) 
(B.L4) 
(B.L5) 
(B.L6) 



{P} (P2)3 4 (47t:)2 V47rT 
{P} (P'r 



2e 



5 1 



8 (471)2 V47rT 
2 / N2e 



^+27- ^+41og2 + 0(e) 
J+27-^+41og2 + 0(e) 



The number 71 is the first Stieltjes gamma constant defined by the equation 



-2 + 27 + 41og2 + 0(e) 



(B.1.7) 

(B.L8) 
(B.1.9) 

(B.1.10) 
(B.1.11) 

(B.1.12) 

(B.1.13) 



^(1+z) = - + 7-7iZ + 0(z2) _ 



(B.1.14) 



B.2 One-loop HTL sum-integrals 



We also need some more difficult one-loop sum-integrals that involve the HTL function 
defined in ( |A.12.4[ ). 

The specific bosonic sum-integrals needed are 



1 

pP2 



Tp 



24 \AnTJ 



1 
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Tpp^ 



(47r)2 \AnT) 

— (^) 
(47r)2 V47rT/ 



+ 27 + 21og2 + 0(e) 



(B.2.2) 



2e 



21og2( -+27) +21og^2 



71^ 



0{e) 



(B.2.3) 



3 (471)2 UttT/ 



(l+21og2)(^+27')-^ + ^log2 



21og^2 + 0(e) 



(B.2.4) 



The specific fermionic sum-integrals needed are 



Tp 



Tp 



1 1 



2 (47r)2 V47rr 



le 



+ 27 + l+41og2 + 0(e) 



, (B.2.5) 



t{P}P^{{FY)^ 



(47r)2 V47rT 



21og2( ^+27) +101og22 



n 



+ 0{e) 



(B.2.6) 



(47r)2 V47rT 



le 



1 7r2 



^+2(7 + 21og2)- + — +41og"2 



(47r)2 V47rT 



le 



87log2-47i + 0(e) 



log2 ( - +27 ) +51og^2 + 0(e) 



(47r)2 V47rT 



le 



l + 27 + 41og2 + 0(e) 



(B.2.7) 



(B.2.8) 
(B.2.9) 



It is straightforward to calculate the sum-integrals ( B.2.1[ |-( B.2.8[ | using the rep- 
resentation ( A.12.6[ l of the function 7p. For example, the sum-integral ( |B.2.2[ | can be 
written 



ipj'^' - ipV' 



P2 



2c2 



(B.2.10) 



where the angular brackets denote an average over c as defined in ( A.12.7| . Using the 
factor Pq in the numerator to cancel denominators, this becomes 



iph^' = iph 



2 2 



(B.2.11) 



The first term in the square brackets vanishes with dimensional regularization, while 
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after rescaling the momentum by p 



tpp' 



Tp 



p/c, the second term reads 



4+2(r 



(B.2.12) 



Evaluating the average over c, using the expression (B.1.61 for the sum-integral, and 
expanding in powers of e, we obtain the result ( B.2.2[ |. Following the same strategy, 
all the sum-integrals ( B.2.1[ |-( B.2.9[ | can be reduced to linear combinations of simple 
sum-integrals with coefficients that are averages over c. The only difficult integrals are 
the double average over c that arises from ( B.2.8| and ( B.2.4| l: 



,l+2e 



21og2-21og2(2-log2)e + 0(e2) , 



(B.2.13) 



l + 21og2 



^log2-log22)e 



(B.2.14) 



B.3 Two-loop sum-integrals 



The simple bosonic two-loop sum-integrals that are needed are 

^PqP2Q2R2 



0{e), 

1 T2 



12 (47r)2 V47rTV 

1 / 7i \ 4e 



i 



PqP^Q^r^ 6(471) 
1 T2 
9 (4 



PQ P2Q2^2R2 

PQ 

PQP^Q^r^ 



i 



71)2 KAtzT) 

7rj2 [atzt) 
8 (471)2 V47rT/ 



1 C'(-l) 
- + 10 - 121og2 + 4^f— ^ + Oie) 

e ^(-1) 



1 T2 



4e 



7.521 + 0(e) 



1 2 ^ 8 4^'(-l) 

- H |-41oe2 + + -^^-^ 

e 9 ^ S*^ 3C(-1) 



where R = — (P + Q) and r = |p + q|. 

The simple fermionic two-loop sum-integrals that are needed are 

1 



i 



1 

{pQ}PW^ 



0{e), 

1 T2 



6 (47r)2 V47rr 



4e 



(B.3.1) 

,(B.3.2) 

(B.3.3) 

(B.3.4) 

0{e) , 
(B.3.5) 



(B.3.6) 



1 Tf-i) 1 

-+4-21og2 + 4|A_^ + 0(e) 



(B.3.7) 
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'{PQ}P2Q2,4 



1 



12 (471)2 V47rT 



4e 



1 11 r(-i) 



{PQ}P2Q2,2R2 

PQ 

{PQ}P2Q2,4 



{PQ}^2p2Q2R2 
^2 



72 (47r)2 V47rT/ 

^ ( H \ 

36 (47r)2 UttT/ 

5 r2 / 77 N4e 

72 (4 



4e 



4e 



{PQH2P2Q2R2 



7r)2 V47rT; 

-— m 
18 (47r)2 \47zTJ 



4e 



7.00164 + 0(e) 



r'(-i) 

l-6T + 6||3^ + 0(e) 



+ 9.55216 + 0(e) 



8.14234 + 0(e) 



(B.3.8) 
(B.3.9) 
(B.3.10) 
(B.3.11) 
(B.3.12) 



where R = — (P + Q) and r = |p + q|. 

To motivate the integration formula we will use to evaluate the two-loop sum- 
integrals, we first present the analogous integration formula for one-loop sum-integrals. 
In a one-loop sum-integral, the sum over Pq can be replaced by a contour integral in 
Po = iPo- 

i F{P) = lim / [F{-ipo,p)-F{0,p)]e^P»nB{po) , (B.3.13) 

J P t]^Q+ J Z7ZI Jp 

where nB(po) = 1/ {e^P° — 1) is the Bose-Einstein thermal distribution and the contour 
runs from — oo to +oo above the real axis and from +oo to — oo below the real axis. 
This formula can be expressed in a more convenient form by collapsing the contour 
onto the real axis and separating out those terms with the exponential convergence 
factor ng ( I Po I ) • The remaining terms run along contours from — oo ± fe to and have 
the convergence factor e'^P". This allows the contours to be deformed so that they 
run from to ±zoo along the imaginary po axis, which corresponds to real values 
of Po = —ipo- Assuming that F{—ipo,p) is a real function of po, i.e. that it satisfies 
F(— zpQ,p) = f *(— fpo/P)/ the resulting formula for the sum-integral is 

i F(P) = / f(P)+ / e{po)nB{\po\)2lmF{-ipo + e,p) , (B.3.14) 
Tp Jp Jp 

where e(po) is the sign of po. The first integral on the right side is over the {d + 
1) -dimensional Euclidean vector P = (Po, p) and the second is over the {d + 1)- 
dimensional Minkowskian vector p = (po/P)- 

The two-loop bosonic sum-integrals can be evaluated by using a generalization of 
the one-loop formula ( B.3.14| : 



i FiP)GiQ)HiR) = \j F{P)G{Q)H{R) 

TPQ 3 JPQ 
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e(po)nB(|po|)2ImF(-/po + e,p)Re / GiQ)H{R) 

Jo 



e(po)nB(|po|)2Imf(-zpo + e,p) / e{qo)nB{\qo\)2ImG{-iqo + e,q) 



xReH(K) 
(cyclic permutations of F, G, H) . 



(B.3.15) 



The sum over cyclic permutations multiplies the first term on the right side by a factor 
of 3, so there are a total of seven terms. This formula can be derived in three steps. 
First, express the sum over Pq as the sum of two contour integrals over po, one that 
encloses the real axis Impo = and another that encloses the line Impo = —Imqo. 
Second, express the sum over qo as a contour integral that encloses the real-(^o axis. 
Third, symmetrize the resulting expression under the six permutations of F, G, and H. 



The resulting terms can be combined into the expression (B.3.15 1. 
The fermionic generalization of ( B.3.15| reads 

i F{P)G{Q)H{R) = f F{P)G{Q)H{R) 

- |^e(poK(|po|)2ImF(-fpo + e,p)Re j^G{Q)H{R) 

- |^e(poK(|po|)2ImG(-zpo + e,p)Re j^H{Q)F{R) 

+ / e(po)nB(|po|)2ImH(-/po + e,p)Re / F{Q)G{R) 
Jp JQ 

+ / e(po)nF(|po|)2Imf(-fpo + e,p) / e(i?o)"F(ko|) 2ImG(-f(?o + e,q) 

J p Jq 

xReH(R) 

Ro=i{po+qo)+£ 

- / e(po)nF(|po|)2ImG(-zpo + e,p) / e(i?o)"B(|<7o|) 2ImH(-^^?o + e, q) 

Jp Jq 

xRef (R) 

Ro=!(Po+'?o)+£ 

e(po)nB(|po|)2ImH(-/po + e,p) / £{qo)nF{\qo\)2lmF{-iqo + £,q) 

I Jq 



Po=-ipo+£ 



xReG(R) 



(B.3.16) 



Ro=!(po+<?o)+£ 

where ftf (po) = 1/ {e^P" + 1) is the Fermi-Dirac thermal distribution. The integrals of 
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the imaginary parts that enter into our calculation can be reduced to 

n{p) 1 



J^e{po)n{\po\)2lm^ 
[ e(po)n(|po|)2Imrp 



f{-ipo + e,p) = [ '-^tJ2f{±ip + e,p) , 
Jp p z ^ 

(B.3.17) 

f{-ipo + e,p) 

= - I pn{p)W{c-^+^'f{±ip + e,ip/c)\ . 

(B.3.18) 



The latter equation is obtained by inserting the expression ( A.12.6| for Tp, using ( B.3.17| , 
and then making the change of variable p — > p/c to put the thermal integral into a 
standard form. 

As an illustration for calculating the simple two-loop bosonic sum-integrals, we 
apply the formula ( B.3.15| to the sum-integral (B.3.21. The nonvanishing terms are 



i 



2 I nB{\po\)27zS{pl-p^) 



+ [ nB{\po\)2nS{pl-p^) [ nBi\qo\)2nSiql-q^)^ . (B.3.19) 

The delta functions can be used to evaluate the integrals over po arid qo- The integral 
over Q is given in ( E.0.1| . This reduces the sum-integral to 



pqP^Q^v^ (47r)2 



4-21og2 



2e f MpI 2e 



nB{p)nB{q) 1 



.2 • 



(B.3.20) 



/pq pq 

The momentum integrals are evaluated in ( D.1.2[ | and ( D.2.3[ |. Keeping all terms that 



contribute through order e^, we get the result (B.3.2 1. The sum-integral (B.3.3 1 can be 
evaluated in the same way: 



PQP^Q-r 



2A 



(47r) 



21og2 



^2e f MPl^-2e + 

Jp p Jpq 



WB(p)wB(t?) q^ 
pq 



. (B.3.21) 



The sum-integral (B.3.5 1 can be reduced to a linear combination of ( B.3.2[ | and ( B.3.3[ l 
by expressing the numerator in the form P-Q = PqQo + {^^ ~ ~ ^^12 and noting 
that the PqQo term vanishes upon summing over Pq or Qq. The sum-integral ( B.3.4| l is a 
little more difficult. After applying the formula ( B.3.15| and using the delta functions 
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to integrate over po, qo, and tq, it can be reduced to 

-.2 r ^^(^\ r I 



i 



Mp) 

p V 



q^ 



Q Q2R2 V r2 ' r2 



nB{p)nB{q) f , p^ , \ r 



pq 



r2 q^ q^J A{p,q,r) 



, (B.3.22) 



where A{p,q, r) is the triangle function that is negative when p, q, and r are the lengths 
of 3 sides of a triangle: 



A{p,q,r) = p'' + q'' + r''-2{p^q^ + q^r^ + r^p 



(B.3.23) 



After using ( E.0.7| -( E.0.9| to integrate over Q, the first term on the right side of ( B.3.22| 
is evaluated using ( D.1.2| . The two-loop thermal integrals on the right side of (B.3.22 1 
are given in ( D.2.11[ |-( D.2.14[ |. Adding together all the terms, we get the final result 
(|R3^. 

As an illustration for calculating the simple two-loop fermionic sum-integrals, we 
apply the formula ( B.3.16| to the sum-integral (B.3.71. The nonvanishing terms are 



+ I nf{\f„\)lTzS{pl-p^) I nf(|,„|)27r«(,J-,2)i . (B.3.24) 

J p J q I 

The delta functions can be used to evaluate the integrals over pQ and q^. The integral 
over Q is given in ( E.0.1| . This reduces the sum-integral to 



i 



{Pe}P2Q2,2 



(471) 



+ 4 -21og2 



nF{p)nF{q) 1 



2e 



Mp) -2e 
P P 



(B.3.25) 



/pq pq 

The momentum integrals are evaluated in ( D.1.2[ | and ( D.2.3[ |. Keeping all terms that 



contribute through order , we get the result (B.3.7(. The sum-integral (B.3.8 1 can be 
evaluated in the same way: 



i 



{PQ}P2Q2,4 



(471) 



21og2 



p 



nF{p)nf{q) q^ 



pq 



pq 



(B.3.26) 



The sum-integral ( B.3.10| can be reduced to a linear combination of (B.3.71 and (B.3.8 1 



by expressing the numerator in the form P-Q = PqQo + {^^ ~ P^ ~ and noting 



that the PqQo term vanishes upon summing over Pq or Qq. The sum-integral ( B.3.9| l is a 
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little more difficult. After applying the formula ( B.3.16| and using the delta functions 
to integrate over po, qo, and vq, it can be reduced to 



{PQ} P2Q2,2R2 



/ Mp) f if 


f Mp) f 1 P^\ 


/p p yQp2Q2R2 


P^-r, h P JqQ'R'W q'J P 




2 2 2 2 

^ — p^ — q^ 




- Mp^q^r) 



'^f(p)'^b('?) f P^ — p 



2 ~^ n2 



P=-ip 



(B.3.27) 



/pq pq vr ry HP'^'^) 

After using ( E.0.7| -( E.0.9| to integrate over Q, the first term on the right side of ( B.3.27| 



is evaluated using (D.1.2|. The two-loop thermal integrals on the right side of (B.3.271 



are given in ( D.2.11[ |-( D.2.14[ |. Adding together all the terms, we get the final result 
(|R3^. 



B.4 Two-loop HTL sum-integrals 

We also need some more difficult two-loop sum-integrals that involve the functions Tp 
defined in ( |A.12.4[ ). 

The specific bosonic sum-integrals needed are 



Tr 



1 r2 



t 



pqP^Q^t^ 



Tr 



48 (47r) 



1 T2 



7r)2 V47rT; 



4e- 



2-121og2 + 4 



576 (47r)2 \AttTJ 



- 19.83 + 0(e) 
1 /26 24 



+ 



477.7 + 0(e) 



P-Q 

PQP2Q2r4 



Tr 



1 r2 

" 96 (471)^ 



(JL_) 



4e 



4 log 2 + 4 



a-i)y e 

(B.4.1) 

r'(-i)\ 1 

921og2 + 4|^)- 
(B.4.2) 



+ 59.66 + 0(e) 
The specific fermionic sum-integrals needed are 

T(PO}PWr^ ~ ~ 48 (471)2 ) 



(B.4.3) 



t 



{PQ}P2Q2,4 



Tr 



1 T2 



^+ 2 + 121og2 + 4-^^ 



+ 136.362 + 0(e) 



(JL_) 



4e 



576 (471)2 V47rT. 
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-2 + (-+521og2 + 4 
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+ 446.412 + 



(B.4.5) 



t 



P-Q 

{PQ}P2Q2,4 



Tr 



1 T2 



96 (47r)2 V47rT 



4e 



J,,f41og2.4-^'(-)^^ 



+ 69.174 + 0(e) 



2 2 

— 



To 



1 



8 (47r)2 V47rT 



4e 



^ + f2 + 27 + ^log2 + 2 



(B.4.6) 



+ 46.8757 + 



(B.4.7) 



To calculate the sum-integral ( B.4.1| , we begin by using the representation ( A.12.6| l 
of the function Tr: 



(B.4.8) 



PQP^Q^r^ fpQP^Q^r^ TpqP^Q^ \Rl + r^c^ , 

The first sum-integral on the right side is given by ( B.3.2[ |. To evaluate the second 
sum-integral, we apply the sum-integral formula ( B.3.15| |: 
1 



PQP^Q^{Rl + r^c^ 

Mp) 



2Re 



Mp)Mr) 

pq 



,-3+2e 



2 2 2 2 2 7 2 

Re^^ — -4 ^+2c"^+^^Re ' ^ ^ 



P^(-ip,p/c) J 

, (B.4.9) 



p<j V ' ^{p + ie,q,rc) A{p + i£,q,rc 

where Tc = |p + q/c|. In the terms on the right side with a single thermal integral, the 



appropriate averages over c of the integrals over Q are given in (E.0.12 1 and (E.0.5 1 



2Re 



QQ2{R2 + r^c^) 



,-3+2e 



Po=-ip+£ 



P^(-ip,p/c) 



(471) 



|^+(4-^log2)l + 16-l^-81og2 + ^log^2 



(B.4.10) 



The subsequent integral over p is a special case of ( D.1.2| |: 



/ Mp)p 

Jo 



2€ 



(B.4.11) 



where (fl)^ = r(fl + b)/T{a) is Pochhammer's symbol. Combining this with ( B.4.10| |, 
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we obtain 

p ' 

T2 



2 Re 



qQ' \Rl + r^c- 



+ {c 



-l+2e 



(47r)2 V47rTy 48 



Po=-ip+e \ 

+ 18 - 121og2 + 4 t; I - + 173.30233 



P^{-ip,p/c)/^J 

. (B.4.12) 



For the two terms in ( B.4.24| with a double thermal integral, the averages weighted by 



are given in (D.2.20 1 and (D.2.271. Adding them to (B.4.261, the final result is 



f ^ \4<r 1 

(47r)2 V47rT/ 48 



1 / 1 
^+ 6-121og2 + 4 t; J I - + 18.66 



(B.4.13) 



Inserting this into (B.4.8 1, we obtain the final result ( B.4.1| |. 

The sum-integral (B.4.2 1 is evaluated in a similar way to ( B.4.1[ |. Using the repre- 
sentation ( A.12.6| l for Tr, we get 



PQ P^Q^r^ fpQ P^Q^r^ \ R2 + r2c2 



(B.4.14) 



The first sum-integral on the right hand side is given by ( B.3.3| l. To evaluate the second 
sum-integral, we apply the sum-integral formula ( B.3.15| l: 



PQP^Qh^Rl + r^c^) 

Mp) 

P P 



Re 



p2 Jf-q2 



QQ2r2(R2 + ,2,2) 



Po^-ip+e P 



-l+2e 



pq 



,-1+2. P 



r2 A{p + i£,q,rc) 



Re 



P^{-ip,p/c) 



A{p + i£, qjc). 

(B.4.15) 



In the terms on the right side with a single thermal integral, the weighted averages 
over c of the integrals over Q are given in (E.0.15 1, ( E.0.16| , and ( E.0.11| : 



Re 



p^ + q^ 



QQV(R2 + r2c2) 
1 



,-l+2e 



Po=-ip+£ 



qQ^R^ 



P^(-ip,p/c), 



{471 



1 /35 31 , ^\ 1 
+ (^-^log2 - 



48e2 ^36 24 

313 2477r2 17 
^ 108 ~ ~576 



18log2 + |log^2 



(B.4.16) 
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After using ( B.4.11| to evaluate the thermal integral, we obtain 



V P \ -h Q ^ 

(47r)2 V47rT; 576 



P^(-ip,p/c) I 



l+('l!^-601og2 + 4f'<-') 



V 3 



84.72308 



, (B.4.17) 



For the two terms in ( B.4.29| with a double thermal integral, the averages weighted by 
are given in ( |D.2.22| |, ( |D.2.29| , and ( |D.2.30| . Adding them to ( |B.4.30[ ), the final result 



IS 



i 



pqP^Q 



2^2 



R2 



^2^2 



T2 



(47r)2 KAttTJ 576 



/314 

71- 



24 

-^-92 log 2 + 4^^^ ^ 



+ 270.2 



(B.4.18) 



Inserting this into ( B.4.28| , we obtain the final result (B.4.21. 

To evaluate ( B.4.3[ l, we use the expression ( A.12.6[ l for 7r and the identity P Q 
(R2 - p2 - Q2) /2 to write it in the form 



t 



P-Q 



Tr 



i 



P-Q 

pqP^Q^t^ 



tpP^^Rr^ 



Tr 



PQ pqP^Q^ \ Rq + r2c2 



C2(l-c2) 



(B.4.19) 



The sum-integrals in the first three terms on the right side of ( B.4.32| are given in ( B.1.4| , 
(B.2.2 1, ( B.3.2[ |, and ( B.3.5[ |. The last sum-integral before the average weighted by c is 



given in (B.4.8 The average weighted by is given in (B.4.27 1. The average weighted 
by can be computed in the same way. In the integrand of the single thermal integral, 
the weighted averages over c of the integrals over Q are given in ( E.0.14| and ( E.0.6| |: 



2Re 



2e^-2e 



-3+2€ 



P^{-!>,p/c), 



{Any 



23 , , ^\ 1 104 
^-41og2 - + - 



7r2-31og2 + 81og^2 



(B.4.20) 



After using ( B.4.11| to evaluate the thermal integral, we obtain 



"b(p) 



2Re 



qQ2 \R§ + rV 



+ c 



l+le 



Po^-ip+£ 



(47r)2 V47rr; 



23 
72 



^log2 ) - + 1.28872 
3 e 



P^(-ip,p/c) fj 

(B.4.21) 
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For the two terms with a double thermal integral, the averages weighted by are 
given in p.2.21| l and p.2.28| l. Adding them to ( |B.4.33| , we obtain 

TpoP^Q^ \ + r2c2 / (471)2 UtzTJ 



17 



1 1 \ 1 
^ - ^ - ^^ogl] - - 0.1917 
72 67r2 3 ^ e 



(B.4.22) 



Inserting this into ( B.4.32| along with ( B.4.27| , we get the final result ( B.4.3| . 

To calculate the sum-integral ( B.4.4| , we begin by using the representation ( A.12.6| l 
of the function 7r: 



{PQ}P^Q^ 



(B.4.23) 



{PQ}P^Q^ \Rl + r^c\ 

The first sum-integral on the right hand side is given by ( B.3.7[ l. To evaluate the second 
sum-integral, we apply the sum-integral formula ( B.3.16| l: 

1 

{pq}pWM+^ 



V V 



2Re 



QQ2(R2 + r2c2 



+ c 



-3+2e 



P^{-ip,p/c) 



pq 



pq 



'A(p + z£, q,rc) 



pq 



pq 



A{p + ie,q,rc) ' 
(B.4.24) 



where = |p + q/c|. In the terms on the right side with a single thermal integral, the 
appropriate averages over c of the integrals over Q are given in ( E.0.5| and (E.0.121. 
The subsequent integrals over p are special cases of ( D.1.2| l and ( D.2.3| : 



"b(p)p 
nf(p)p 



-l-2e 



-l-2e 



(i)-4.(^)2. a-i+4e)^ 

n^{p)p-^-^' 



(i)-2.(i)-. a-1 



r2 

12 ' 



1-2 



-l+4e 



(B.4.25) 



This yields 



-2/ ^Re 



P P JQQ^\Rl + r^c^ 



+ / MpI /^-1+2. 



(47r)2 V47rT/ 48 



^ - 6- 121og2-4 



- + 70.122 



P^(-ip,p/c)/^ 

(B.4.26) 



For the two terms in ( B.4.23| l with a double thermal integral, the averages weighted 
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by are given in (D.2.23 and (D.2.31i. Adding them to (BA26 1, the final result is 



7r)2 \AttTJ Us 



(47 



^ - 6-121og2-4 t; J I -+51.9306 



. (B.4.27) 



Inserting this into ( B.4.23| ^ we obtain the final result (B.4.41. 

The sum-integral (B.4.5 1 is evaluated in a similar way to ( B.4.4[ |. Using the repre- 
sentation ( A.12.6| l for Tr, we get 



{PQ}P2Q2,4 



(B.4.28) 



The first sum-integral on the right hand side is given by ( B.3.8| l. To evaluate the second 
sum-integral, we apply the sum-integral formula ( B.3.16| l: 



i 



{PQ}P2QV(R2 + r2c2 
P P 



Re 



p2 + 



QQ2r2(R2 + r2c2) 



+ C 



-l+2e 



Mp) -2 
P P 



P^(-ip,p/c) 



+ 



pq pq r2 A{p + i£,q,rc 

-i+2e f Mp)Mq) p 

Jpq pq 



2 _L^2 

^Re 



2 2 2 

rj. — p'^ — q^ 



A{p + ie,q,rc 



(B.4.29) 



In the terms on the right side with a single thermal integral, the weighted averages 
over c of the integrals over Q are given in ( E.0.11[ |, ( E.0.15[ l, and ( E.0.16[ |. After using 
( B.4.25| to evaluate the thermal integral, we obtain 

^2 



nf(p) 



Re 



p^ + q^ 



Q Q^r^ \R^ + r2c2 



Mp) 1 Ll+2e 

p p p^\ 

j2 



(471)2 



Q Q'K2 
\4nTJ \576 



1 /34 



3 361og2-4^^_j 



t V M _ _^ 229.354 
e 



(B.4.30) 



For the two terms in ( B.4.29| l with a double thermal integral, the averages weighted 
by are given in ( |D.2.25| |, ( |D.2.33| |, and ( |D.2.34| |. Adding them to ( |B.4.30| |, the final 
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result is 



(47r)2 KAttTJ \576 



/118 



52 log2-4 



+ 90.9762 



(B.4.31) 



To evaluate ( B.4.6[ >, we use the expression ( A.12.6[ > for 7r and the identity P-Q 
{R^ - p2 - Q2) /2 to write it in the form 



P-Q 

{PQ}P2Q2,4 



P-Q 



{PQ}P2Q2,4 

V 1 



C2(l-c2) 



{PQ}P2Q2,2 

(B.4.32) 



The sum-integrals in the first three terms on the right side of ( B.4.32| are given in ( B.1.8| , 
( [RZ2| , ^3j7\ , and ( |B.3.10f . The last sum-integral before the average weighted by c is 



given in (B.4.24 1. The average weighted by is given in (B.4.27 1. The average weighted 
by can be computed in the same way. In the integrand of the single thermal integral, 
the weighted averages over c of the integrals over Q are given in ( E.0.6[ l and ( E.0.14[ |. 
After using ( B.4.25| to evaluate the thermal integral, we obtain 



"f(p) 



Re 



(4 



V JQQ^\Rl + r^c^ 



,l+2e 



/ 7 1 \ 1 



P^(-ip,p/c) 

(B.4.33) 



For the two terms with a double thermal integral, the averages weighted by are 
given in p.2.24| | and ( |D.2.32| |. Adding them to ( |B.4.33| , we obtain 



t 



{PQ}P^Q^ \R§ + r2c2 



(4 



— m 



4e 



0.1359 



(B.4.34) 



We finally need to evaluate (B.4.71. Applying ( |B.3.16| | gives 



r{PQ}PWoR'^ 



Re 



p V Jp P 

pq pq 

nB{p)nF{q) 

pq pq 



Re 



2 2 

p — q 



Q\Q^r^Rl + r^c^) 



2 Jl 



P 



c P -q ^-l+2e 



q^ A{p + ie,q,rc 



Po=-ip 

I 



2 2 2 2 2 

Re( l^ji^P_Zi_c-i+2. 
q^ A{p + ie,q,rc) 



111 



APPENDIX B. FOUR-DIMENSIONAL SUM-INTEGRALS 



MP>^BMj^g / p2 - ^2^2 _ p2 _ ^2 



pq 



A{p + ie,q,rc) 



(B.4.35) 



In the terms on the right side, with a single thermal factor, the weighted average is 
given m Eq. ( |E.0.17| , After using Eq. ( |B.4.25| l to evaluate the thermal integral, we obtain 



2 2 

p -r 



Q\Q^r^{Rl + r^c^) 



(471) 



24 ) 



(B.4.36) 



The terms with two thermal factors are given in Eqs. p.2.26) , p.2.35[ > and ( |D.2.36| |. 
Adding them to (|B.4.36[), we finally obtain dBATb. 



B.5 Three-loop sum-integrals 



The three-loop sum-integrals needed are 

1 



i 



pqrP^Q^R^{P + Q + R)^ 



1 T"* / U N6e 



24 (471)2 V47rr 



1 91 r'(-i) 2^'^-^Ko(,) 



(B.5.1) 



PQRP2Q2R4(Q_R)2(R_p)2 

216 (47r)2 V47rT 



' + '^ + ^7 + ^^-^^ + 0(e) 
e 22 ir 11 ^(-1) 11 ^(-3) 



(B.5.2) 



{PQR}P^Q^RHP + Q + RV 



1 



96 (47r)2 V47rT 



be 



(B.5.3) 



1 



PQ{R}P2Q2R2(p + Q + R)2 



1 / \6e 



192 (47r)2 V47rT 
QrR 

{P}QrP2q2r2(P + Q)2(P + R)2 
1 / ^ \6e 



l + ^-|iog2 + 85(z2)_,n:^^o(.) 



a-i) a-3) 



(B.5.4) 



384 (47r)2 V47rT 



1 361 , 76, ^ .r(-3) : 



(B.5.5) 



P{QR}P2Q2R2(p + Q)2(p + R)2 
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3456 {4ny 



\4tzTJ 



6e 



1 23 6 192 28rM)_4r(z3: 
e^5^5^ 25'°S^^5a-l) 5^-3) 



+ 0{e) 
(B.5.6) 



The three-loop sum-integrals were first calculated by Arnold and Zhai and calculational 
details can be found in Ref. ||5|. 
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Three-Dimensional Integrals 



Dimensional regularization can be used to regularize both the ultraviolet divergences 
and infrared divergences in three-dimensional integrals over momenta. The spacial 
dimension is generalized to = 3 — 2e dimensions. Integrals are evaluated at a value 
of d for which they converge and then analytically continued to d = 3. We use the 
integration measure 

We require one integral that does not involve the thermal distribution function. The 
momentum scale in these integrals is set by the mass m = mjj. 



C.l One-loop integrals 



The one-loop integral is given by 



Specifically, we need 



/ 1- 



/p (p2 -I- m2)" 
1 

8n 



r(i)r(n) 



-m 



= [ log{p^ + m^) 
Jv 



6n \2m) 



h 



4:71 \2mJ 



Sum \2mJ ^ ^ ' 



(C.1.1) 



(C.1.2) 

(C.1.3) 
(C.1.4) 
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C.2 Two-loop integrals 

We also need a few two-loop integrals on the form 

Kn 



1 



pq p2 + {q^ + m^Y (p + q)2 ' 
111 



/pq + + [(p + q)^]" ' 

Specifically, we need }i, }2, and which were calculated in Refs. lUTllTZl : 



h 



K2 



1 (Ji_ 
4(47r)2 \2mJ 



4e 



2 + 0{e) 



4(47r)2m2 \2mJ 



8m2(47r)2 V2m/ 



(C.2.1) 
(C.2.2) 



(C.2.3) 
(C.2.4) 
(C.2.5) 



C.3 Three-loop integrals 

We also need a number of three-loop integrals. The specific integrals we need are listed 
below and were calculated in Refs. IUTI [121 . They are special cases of more general 
integrals defined in Ref. 



i 



pqr (p2 + m2)(<j2 + ^2)^2(p + q + r)2 



m 



(JL) 

2(47r)3 \2mJ 



6e 



8 + 0{e) 



pqr (p2 + m2)(<?2 + ni^) (p - q)2(q - r)2(r - p)2 



m f }i 
4(47r)3 V2m; 
(r2 + m2)2 



8 + 0(e) 



pqr (p2 + m2)(^j2 + ^2) (p - q)4(q _ r)2(r - p)2 



m 



L 



4(47r)3 \lm) 



6 + 0(e) 



pqr (p2 + jtt2)(^2 + ^2)(^2 + ^2) (q_r)2(r-p)2 



(JL\ 



6€ 



12 



1 1 
pqr (p2 + m2)(^j2 + m2) (p - q)2(q - r)2(r - p)2 



8m(47r)3 V2m/ 



6e 



2 + 0(e) 



(C.3.1) 



(C.3.2) 



(C.3.3) 



(C.3.4) 



(C.3.5) 
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pqr {p^ + m^){q^ + m^){r^ + m^)^ (q-r)2(r-p)2 



i (JL-) 

'(47r)3 \2mJ 



6e 



[ ^ 

7pqr + m2) + m2) [(q - r)2 + m^] [(r - p)2 + m^] 



" ~ (47r)3 V2mJ 

X 



6e 



8 -41og2 + £>(£) 



(P - q)^ 



pqr (p2 + m2) ((^2 + ^2) [(q _ r)2 + f„2] [(^ - p)2 + m^] r2 

2m / jM \ 6e 



(47r)3 V2m 



[l-21og2 + C>(e)] , 
1 



(C.3.6) 



(C.3.7) 



(C.3.8) 



(p - q)^ 



pqr + m2) ((^2 + ^2) [(q _ i.)2 + ^2] [(f _ p)2 + „2] ^4 



3m / f( 



(47r)3 V2m 



6e 



l--log2 + 0(e) 



(C.3.9) 



Ipqr (p2 + m2) {q^ + m2) [(q - r)2 + m^] [(r - p)2 + m^] 

^ (P - q)^ 

pqr (p2 + m2) ((j2 + ^2) [(q - r)2 + m2] [(r - p)2 + m2] 



i 



i 



3m(47r)3 v2m 



4m3(47r)' 



6e 



[l-log2 + C>(e)] , (C.3.11) 

1 1 

pqr (p2 + m2) {q2 + m2) [(q - r)2 + m2] [(r - p)2 + m2] r2(p-q)2 

1 1 



Xqr (jf^^ + m^) if + m2) [(q - r)2 + m2] [(r - p)2 + m2] 



(C.3.13) 



Finally, we need the combination 

1 



i 



(p-q)^ 



pqr (p2 + m2)(^2 + ^2)(^2+^2) (q _ r)2(r - p)2 

{q^ + m^) 1 



+ 



L 



pqr (p2 + m2) [(r - p)2 + m2] [(q - r)2 + m2] r2(p - q)2 
2m / }i 



(C.3.14) 
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Three-Dimensional Thermal 
Integrals 



The three-dimensional thermal integrals involve the Bose-Einstein distribution nsip) = 
1/ (e^P — 1) and the Fermi-Dirac distribution nf{p) = 1/ (e^P + 1). 



D.l One-loop integrals 



The one-loop integrals can all be obtained from the general formulae 

Mp)^2a _ a2 + 2a-2e)r(2 + 2a-2e)r(l) 



P P 



1-2 



-l-2a+2e 



r(i-e) 

V V 



(eV)'T2+2«-2^ (D.1.1) 
(D.l. 2) 



D.2 Two-loop integrals 



The simple two-loop thermal integrals are 



pq 



nB{p)nB{q) j_ 



r nBip)nB{q) 
Jpq pq r4 

r nF{p)nF{q) 1 
Jpq pq r2 

nF{p)nF{q) f 

pq 



1 



4e 



, 4e 



4 (47r)2 V47rT 

— m 

(47r)2 ^4711/ 

T2 / II 



1 14 C'(-l) 

_ + _+41og2 + 4|i-J+0(.)^ 



1 1 ITf-l) 



(471)2 V47rTy 3 

1 T2 / ^ \4e 

~ 36 (471)2 V47rTJ 



[l-log2 + C>(e)] , 



(D.2.1) 

I , 
(D.2.2) 
(D.2.3) 
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pq 



5 + 67 + 61og2-6 

1 / ^( \4e 



pq 



nB{p)nF{q) 
pq 



36 (47r)2 V47rT 
X 7-67-181( 
1 / ^ 



4e 



18 (471)2 V47rT 
X 1-67 -12 log 2 + 



a- 


1) 




-1) 




-1) 


^.^'^ 


-1) 




-1) 



(D.2.4) 



+ 29.509e + O(e) 



+ 31.098e + O(e) 



, (D.2.5) 



. (D.2.6) 



We also need some more complicated two-loop thermal integrals that involve the 



triangle function defined in Eq. (B.3.231: 

nB{p)nB{q) 



pq P'? q^A{p,q,r) 



L 



nB{p)nB{q) 



pq P"? Mp^q^r) 



nBip)nB{q) p^ 



pq P'? q^A{p,q,r) 



nB{p)nB{q) p^{p^ + q^^ 



pq P^ r^Aip,q,r) 



nF{p)nf{q) 



pq pq q^A{p,q,r) 



L 



nf{p)nF{q) 



pq pq A(p,i?,r) 



7 



4e 



48 (47r)2 V47rT/ 



1 /22 ^ 



35 /e 



V 7 

40.3896 + 0(e 



■a-i) 



(D.2.7) 



1 r2 / ^ \4e 



24 (47r)2 V47rf/ 

t2 



-2+2 1 + 7 + ^^^ 



71"- 



4 + 47 + Y - + 4(1 + 7) 
2S^ + 0(e) 



^3) r2 / ^ 

240 (47r)2 Katct) 
+ 2||l + 0(.) 

1 r2 / \4e 

48 (47r)2 Vi/rT/ 
- 86.46 + 0(e) 

7 r2 / 7i \4e 



(D.2.8) 



1 ^ -3 

e ^ -3 



1 /14 

-2 + 1^ + 107-6 



(D.2.9) 

r(-l)\ 1 



V 3 



(D.2.10) 



96 (47r) 
f2^ 

f 47.2406 + 0(e) 

1 T2 / \4e 
"48 (47r)2 Vi/rT/ 
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pq 



pq 



nj:{p)nj:{q) p^{p^ + q^) 
pq r^A{p,q,r) 



nF{p)nj:{q) p^ 



pq pq Mp^q^r) 



nF{p)nB{q) p^ 



pq pq Mp.q.r) 



nF{p)nB{q) 



pq 



pq q^A{p,q,r) 



pq 



I 



Mp>^iq) 
pq 



q^A{p, q,r) 
nF{p)nB{q) 

pq pq ^{p,q,r) 



4 log 2 - 471 + 4(1 + 7 + log 2) 



a-i) 



49^(3) T2 / ^ \4e 



(D.2.12) 



1920 (47r)2 V47rT; 

,r(-3) , ,r(3) 

~96 (47r)2 Vi/rT/ 



2 + 21og2 



26 



(D.2.13) 



+ ^ + 107 + 101og2 



96 (471)2 V47rT; 



41.1580 + 0(e) 



(D.2.14) 



4e 



+ 2 l + 7 + log2 



r2 , „ ^ 4e 



96 (4tiV \4ttTJ 



37.0573 + 0(e 

" 1 



(D.2.15) 



96 (4; 



+ 2 l + 7-log2 



7^3) T2 



+ 19.2257 + 0(e) 



(D.2.16) 



7r)2 V47rT/ 



4e 



+ 2 



1920 (47r): 



1 2 

-+2--log2 



f(-3) 

'2 



+ 2 



+ 0(e) 



24 (47r)2 V47rT/ 



4e 



+ 4 + 27-51og2 



+ - + 18.1551 + 0(e) 



■-— m 
96 (47r)2 V47rTy 



(D.2.17) 

703) 
80 

(D.2.18) 



4e 



+ 2 l + 7 + 51og2 



+ 84.2513 + 0(e) 



(D.2.19) 



The most difficult thermal integrals to evaluate involve both the triangle function 
and the HTL average defined in ( A.12.7| . There are two sets of these integrals. The first 
set is 



pq 



pq 



nB{p)nBiq) 

pq 

nB{p)nB{q) 

pq 



Re( 



Re 



22 2 2 
r c — p — q 

A{p + ie, q, rc) 

21 2 2 
r^c — p^ — q^ 

A{p + ie, q, rc) 



^ [0.138727+ 0(e)] , 



(D.2.20) 



j2 



4e 



67r2 (47r)2 V47rr 
^1 



6.8343 + 0(e) 



(D.2.21) 
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nB{ p)nBiq) 

pq 



Re( 



2 2 2 2 

r^c^ — — q^ 



71' 



-1 T2 



A{p + ie,q,rc) 247r2 (47r)2 V47rTy 



4e 



2 2 2 2 

r c — p — q 



nFip)nF{q) ^2 

pq pq \ A{p + i£,q,rc 



X 
j2 



1 



+ 15.3782 + 0(e) 



2 2 2 2 

r c — p — q 
A{p + ie, q, rc) 

Re/ 



pq 

nj:{p)nj:{q) ^ 
pq 

nBip)nF{q) P^-q^ / r^c^ - p^ - q^ 



(47r)2 
T2 

(47r: 



2 2 2 2 

r c — p — q 
A(p + i£,q,rc) 



[0.01458 + C(e)] , 
2 [0.017715 + C(e)] , 

T-2 

-0.011578 + 0(e)] 



(47r)2 ' 



[0.17811 + 0(e)] 



pq \A.{p + ie,q,rc) / (47r)2 

The second set of these integrals involve the variable Vc = |p + q/c|: 



(D.2.22) 
(D.2.23) 
(D.2.24) 
(D.2.25) 
(D.2.26) 



nB{p)nB{q) 



2 7 

p^ — q'^ 



Jpq pq \ A{p + i£,q,rc) 



1 T2 



8 {4ny 



\AttTJ 



4e 



nB{ p)nB{q) 

pq 



A{p + ie,q,rc) 



1 T2 



24 (4/1)2 KAttTJ 



4e 



^ +13.442 + 0(e) 



^+ 16.381 + 0(e) 



Ipq 



nB{ p)nB{q) f 

pq q^ 



Re/ci+2^4f 



2 2 

p"^ — q^ 



A{p + i£,q,rc 



4e 



Ipq 



pq 



1 t2 

48 (4/1)2 V47rT 



V4/rr/ 



^ +6.1227 + 0(e) 



(^2 A(p + ze,(^,rc) 



5-81og2 T2 / n \4e 
144 (4/r)2 \4/rTJ 



^ +100.73 + 0(e) 



nF{ p)nB{q) / ^-i+2e 'c 



pq 



2 2 2 

rt — p — q 



A{p + i£,q,rc) 



^ [0.19678 + 0(e)] , 



Mp)nB{q) 



2 2 2 

n — p^ — q'^ 



f -FKFJ'^BK<1J ^^fl+2e^ 

Jpq pq \ A{p + te,q,rc) 



(Any 



[4.8368 X 10-^ + O{e)] , 



Ipq 



2 2 2 

n — p^ — q^ 



nFip)nB{q) f_ / c^+^e ^ 

pq q^ \ A{p + i£,q,rc) 



(D.2.27) 



(D.2.28) 



(D.2.29) 



(D.2.30) 



(D.2.31) 



(D.2.32) 
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1 T2 



nV KAttTJ 



4e 



96 (47r) 



-7.77235 + 0(e) 



(D.2.33) 



pq 



A{p + i£,q,rc 



ll-81og2 T2 



288 

nf{p)nF{p) 



(47r)2 V47rT 



4e 



+ 7.79813 + 0(e) 



(D.2.34) 



pq 



pq 



Re 



2 2 2 2 2 

A{p + i£,q,rc) 



1 r2 



24 (47r)2 V47rT 



4e 



^+ 2 + 27 + 21og2 + 2^^ 



40.316 + 0(e) 



e 

(D.2.35) 



nB{p)nF{p) / r. 



pq 



pq 



12 2 12 



q^ A{p + ie,q,rc 



1 I }l \4e 

12 (IttP V47rT; 



1 / C'(-l)\ 1 

^+ 2 + 27 + 41og2 + 2^^ 



52.953 + 0(e) 



e 

(D.2.36) 



The simplest way to evaluate integrals like ( D.2.1| l-( D.2.6| l whose integrands factor 
into separate functions of p, q, and r is to Fourier transform to coordinate space where 
they reduce to an integral over a single coordinate R: 



/ f{v)8m{r) = I f{R)g{R)m 

Jpq JR 



(D.2.37) 



'pq 

The Fourier transform is 

f{R) = I .'P V(P) , 
Jp 

and the dimensionally regularized coordinate integral is 

R V 471 ' ' 

The Fourier transforms we need are 

'3 , „, ^\ / o \ 3+2a-2e 



(D.2.38) 



(D.2.39) 



p2^g/pR 



1 r(f + a-e) 2 e (2 



n ia)l{-oc] 



2\' 



R 



(D.2.40) 



1 1 



P P 47rr(2 



I J I dvV^'^-^-'n{p)]^^_,{pR) 



(D.2.41) 



If a is an even integer, the Fourier transform ( D.2.41| is particularly simple in the limit 
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3: 



p p 

MpI 2 ^tp.R 
P 

Mp) gip.R 
p p 

Jp p 



( 



COth X 



— — — coth X — coth X ^ 



2R V 

^ - — cschx ) , 



AtzR \x 

7TT^ 



I csch^x 



r cschx 



(D.2.42) 
(D.2.43) 
(D.2.44) 
(D.2.45) 



where x = nRT. We can use these simple expressions only if the integral over the coor- 
dinate R in ip.2.37} converges for d = 3. Otherwise, we must first make subtractions 
inside of the integrand to make the integral convergent. 

The integrals p.2.2^ - ^DT^ can be evaluated directly by applying the Fourier 
transform formula ( D.2.37[ | in the limit e — ?► 0. The integral ( D.2.1[ | however requires 
subtractions. It can be written 



nB{p)nB{q) 1 



Mp) f Mq) _ ^ 
pq p \ q q 



+ T 



Mp) 



qq r 



(D.2.46) 



'pq pq Jpq p \ q r/ Jp p 

In the second term on the right side, the integral over q is proportional to p^^^^^ 



, so 



the integral over p can be evaluated using ( D.l.l[ |. This first term on the right side is 
convergent for d = 3 so it can be evaluated easily using the Fourier transform formula 



(D.2.371. The integral over R reduces to a sum of integrals of the form dx x"' coth" x. 
Although the sum of the integrals converges, each of the individual integrals diverges 
either as x — )■ or as x — ?► oo. A convenient way to evaluate these integrals is to use 
the strategy in Appendix C of Ref. ||5l. The integrals are regularized by using the 
substitution 



dx x"' coth" X 



r(i + ^) 

2S 



dx x"'^^ coth" X . 



(D.2.47) 



The divergences appear as poles in 3 that cancel upon adding a convergent combination 
of these integrals. 

The integrals ( D.2.7[ |-(D.2.9 1 and ( D.2.11| -( |D.2.13| can be evaluated by first averag- 
ing over angles. The triangle function can be expressed as 



4pY(l -cos^^) , 



(D.2.48) 



where 9 is the angle between p and q. For example, the angle average for ( D.2.7| and 
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p.2.11| | 



IS 



A{p,q,r) 



w{e) 



p q 



7-1 p^q^ 



After integrating over x and inserting the result into (D.2.71 and ( D.2.11[ |^ the integral 
reduces to 



n{p)n{q) 



pq 



pq q^A{p,q,r) 



n{p)n{q) 



1 - 2e »2 7 - 6e 1 , _ 



/pq pt? 

The integrals over p and q factor into separate integrals that can be evaluated using 
( |D.1.1| | and p.l.2[ |. After averaging over angles, the integrals p.2.8| , ( |D.2.9| , ( |D.2.12[ ) 
and ( D.2.13| l reduce to 

^2 



n{p)n{q) 



pq 



pq A{v,q,r) 



1 — 2e f n{p) f n{q) 1 



q ^ r 



n{p)n{q) 



4e 
1 -2e 
8e 



q q q* 



/pq q^A{p,q,r) 

which again can be evaluated using ( D.l.l| and ( D.1.2| . 

The integral ( D.2.10| and ( D.2.14| can be evaluated by using the identity 



(D.2.51) 
(D.2.52) 



_ _ _ 

\r^A{p,q,r)/.,^ 2e\r^ 



+ 



l-2e 1 



p q 



(D.2.53) 



The identity can be proved by expressing the angular averages in terms of integrals 
over the cosine of the angle between p and q as in ( |D.2.49[ |, and then integrating by 
parts. Inserting the identity ( |D.2.53| into (D.2.101 and (D.2.141, the integrals reduce to 



n{p)n{q) p^{p^ + q 



pq 



pq 



r^A{p,q,r) 2e 



n{p)n{q) p^ 

pq PR ''^ 
1 —2e f n{p)n{q) 1 



8e 



pq q' 



(D.2.54) 



The integral in the first term on the right is given in (D.2.2 1 and (D.2.4 1, while the second 
term can be evaluated using (D.1.1 1 and ( D.1.2[ |. 

The integral ( |D.2.23| | can be evaluated directly in three dimensions by first aver- 
aging over c and x, and then integrate the resulting functions numerically over p and 
q. 

To evaluate the weighted averages over c of the thermal integrals in Eqs. ( D.2.20| - 
( D.2.22[ | and Eqs. ( D.2.24[ |-( D.2.26[ l, we first isolate the divergent parts, which come 
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from the region p — — > 0. We write the product of thermal functions in the form 

nip)niq) = (n{p)niq) - + , (D.2.55) 

where s = (p + (^)/2. In the difference term, the HTL average over c and the angular 
average over x = p • q can be calculated in three dimensions: 

\ A(p + /e,^,rc)4 4p^ ^|p-<j| 2(p2 - 6^ 



A{p + ie,q,rc)/^^^ 3{p^ - q^Y llpq ^ \p - q\ 

(3p2+(?2)(p2 + 3tj2) p 



6{p^ — q^)^ q 



log^^, (D.2.57) 



Re 



r2 A(p + ie, q, rc) f^,, 3(p2 - q^f \ 1 ^ pq 

- 4^ log - 4^ log P ) , (D.2.58) 

Apq \p-q\ p^-q^ ^ q ' 

p2 — q^ r^c^ — p2 _ ^2 \ ^ 2 



A(p + ze, (7, rc) / 4:pq{p^ — q^) 

- 2p^log 



c,x 

|p2_„2n 



(D.2.59) 



pq 

The remaining two-dimensional integral over p and can be evaluated numerically: 
/ /^. ^V-P^-)^ ^ [5.292X 10-] , 



V pq p'^q^ J \ A{p + ie,q,rc) {Any 

(D.2.60) 



r / nB{p)nB{q) s^n^^js) \ /, r^c^ -p'-q' \ ^ [3292 x10-^1 

7pq V pt? / \ A{p + ie,q,rc) (47r)2 L J 



nB{p)nB{q) 


sM(s) 


pq 


p^q^ 


nBip)nB{q) 


s2n|(s) 


pq 


pY 


Mp)M<i) 


s2n|(s) 


pq 




nF{p)nF{q) 




pq 


p2(^2 


nF{p)nF{q) 


S^Hp (s) 


pq 


p2^2 



(D.2.61) 



ipqV P*? pY y \ A(p + ze,t7,rc)/^ (47t)2 l J 

(D.2.62) 

/^4 r2c2-p2-^2 \ ^ T2 ^^^^^^ 

\ A(p + ze,(?,rc)/^ (47r)2 L J 

(D.2.63) 

r / nr{p)n,{q) s^nUs) \ q^ /, r^c^ -p'-q' \ ^ ^7 792 x IQ-'] 

ipqV pq pV J \ A(p + ze,^,rc)/^ (47r)2 L J' 

(D.2.64) 

f ( nB{p)n,{q) _ s'nB{s)nF{s) \^J f-ci' rh^ - P' - \ = ^ [0.178II] . 
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(D.2.65) 



The integrals involving the n^(s) term in (D.2.55 ' are divergent, so the HTL average 
over c and the angular average over ^ = p • q must be calculated in 3 — 2e dimensions. 
The first step in the calculation of the n^(s) term is to change variables from p and q to 

s = {p + q)/2,^ = Apcj/ (p + cj)^, and x = p • q: 



2^2 fiP'^'^) 



64 



pq p^q 



{Any 



X ^'rf/5r''(l - i6)-'/'(/(s+,s_,r) +/(s_,s+,r))^ , 

(D.2.66) 

where s± = s[l ± ^/l — fi] and r = s[4 — 2/3(1 — x)]^''^. The two terms inside the 
average over x come from the regions p > q and p < q, respectively. The integral over 
s is easily evaluated: 



dss^-*'nl(s 



r(2-4e) [ai-4e)-a2-4e)]r 



2-4e 



dss^-*'^n^(s) = r(2-4e) -(1 - 2*^)^(1 - 4e) 



fl -2 



-l+4e 



)a2-4e) 



.2-4e 



/ dss^-^'nFis)nBis) = l-^+'^'T (2 - Ae)a2 - Ae)T 
Jo 

It remains only to evaluate the averages over c and x and the integral over j6. 



(D.2.67) 

(D.2.68) 
(D.2.69) 



The first step in the calculation of the n^{s) term of ( |D.2.20| |, ( |D.2.21[ ) and ( |D.2.24[ ) is 
to decompose the integrand into two terms: 



2 2 2 2 

r c — p — q 



A{p + ie,q,rc) 2 ^ [p + i£± q)^ — r^c^ ' 

The weighted averages over c gives a hypergeometric function: 

1 1 il 



(D.2.70) 



(p + ie ± qY — f^c^ 



3-2e (p + fe±ij)2 



(p + f£±ij)2 



(D.2.71) 



5 1 

2' ^ 



(p + Z£ ± - r2c2 (3-2e)(5-2e) {p + k±qf 



{p + i£±qYj 

(D.2.72) 



In the -\-q case of ( D.2.71[ | and ( D.2.72[ |, the /e prescription is unnecessary. The 
argument of the hypergeometric function can be written 1 — ^y, where y = (1 — x) 12. 



125 



APPENDIX D. THREE-DIMENSIONAL THERMAL INTEGRALS 



After using a transformation formula to change the argument to ^y, we can evaluate 
the angular average over x to obtain hypergeometric functions with argument j6. The 
averages over x of ( D.2.71| l and ( D.2.72| l are 



^ 1 



3-2e 
2e 



2-2£,l + £ 

'3 



/5 



(l),(l)_2.(2)-2.(i)-. ^_,j. (l-2e,l-e 



(l)-.(2)-3. 



2-3e 



(D.2.73) 



5-2e 
2e 



I 2 - 2e, 1 + e 



(l),(l)_2.(2)-2.(|)-. ^_,^ 1 - 2e, I - e 



(l)-.(2)-3. 



2-3e 



(D.2.74) 



where {a)}, is Pochhammer's symbol which is defined in (F.0.3 1. Integrating over j6, we 
obtain hypergeometric functions with argument 1: 

^2 



Jo 



1/2 



(p + (?)^ — r^c^ 



1 (l).(2)-2. 
4e (1)-. 



(1)-2.(1)-. ^( \-2e,\-e,\,\ 
(2)-2.(2)-2.(l). \\-2e,2-2e,\^e 



(l)-3.(l)-2.(i)-. ^ (\-?>e,\-2e,\-e 



(i)_3<r(2)-3. ' V 2-3e,2-3e 



(1-/3) -1/2 

JO 



-3<r 
^4 



(D.2.75) 



(p + q^- — r^c^ 



c,x 



4e(3 - 2e) 



2e 



1 - 2e, 1 - e, f , 1 
(2)-2/ U -2e,2-2e,l + e 



(1)- 



(l),(l)_2.(2)-2.(i)-. (l)-3. ^ (l-3e,l-2e,l-6 



{l)-e{2)-3e ilUe 



\ -3e,2 -3e 



(D.2.76) 



Expanding in powers of e, we obtain 



Jo 

/^/3r'^(l-/3) 

^0 



-1/2 



-1/2 



(p + q)'^ — r^c^ 

(p + q)'^ — r^c^ 
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c,x 



w- 
24 

712 

72 



(1+ 3.54518 e) , (D.2.77) 



(1 + 10.8408 e) . (D.2.78) 
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In the —q case of ( D.2.71| and ( D.2.72| |, the argument of the hypergeometric func- 
tions can be written (1 — fiy) / (1 — j6 ± ie), where y = (1 — x) /2 and the prescriptions 
+/£ and —is correspond to the regions p > q and p < q, respectively. These regions 
correspond to the two terms inside the average over x in ( D.2.66[ |. In order to obtain 
an analytic result in terms of hypergeometric functions, it is necessary to integrate 
over j6 before averaging over x. The integrals over j6 can be evaluated by first using 
a transformation formula to change the argument of the hypergeometric function to 
— — — j6) and then using the integration formula ( F.0.10| to obtain hypergeo- 

metric functions with arguments y or 1 — y: 



-3/2, 




{l)-2e [ 



e 
3 



2e (1), 



3 



5-2£ Ol^f 

^ {l)-2e 

2{5-2e) (i; 



3^ ill 
4(5_e)(3_e)(i 




3e(3e-l)(3e-2) 



(D.2.79) 




After averaging over x, we obtain hypergeometric functions with argument 1: 



/ dBR 

lo 



-2e, 



/3) 



-1/2 



1 



2£ 



{p + ie — qY — r^c^ 
\-e,\-2e,\ 



4e (1)^2, \^ 2-2e,l + e 
1 (2)_2.(l).(i)-. 



1 



i — e - — 2e 1 



2e (i)_,(i),(i)_2. 



2e i 
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,(2)_2.(l)-2.(l).(i)- 



free \-i-^t \^ ;-^t \^jt \ 2)-€ p ( ^ 3e, 2 



8e(l-3e) (l)-e(2)-3. 



2-3e,2-3e 



(D.2.81) 



(p + ie — qY - r^c^ j 
i),(l)-,(l),(|)-, ^/l-e.l-2g.l.f 



4e(l-2e)(3-2e) {\)-re{\)ie ^ \ l-le,\-, e,^ 

1 ^1 



2 "-'2 



e4 -2e,2 



.!7re V 2 



(i),(i)_2.(i)Ki)-.(f)- 



16e(l -3e)(3-2e) 



(l)2.(3)-3. 



P ( 1 — 1 — 3e, 2 



3-3e,2-3e 



(D.2.82) 



Expanding in powers of e and then taking the real parts, we obtain 



(p + ie — qY ~ ''^'^^ 



Res^ / d|6j6-2e(i_ fl)-i/2 
Jo 

Res2 /^d^j6-2^(l-^)-i/2 ^ ^, , , 

Jo \ (p + ze - (j)2 - r2c2 



c,x 



12 + 



(1 + 0.34275 e) , 

(D.2.83) 
(1 + 1.10518 e) . 

(D.2.84) 
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Inserting the sum of the integrals ( D.2.77| l and ( D.2.83| into the thermal integral ( D.2.66| 
and similarly for the integrals weighted by c^, we obtain 

T2 









p2(j2 






1 


p2q2 


/ 




7pq 


p2q2 



Re 



2 2 2 

r c — p — 




A(p + /£, 


re) 


2 2 2 

r c — p — 




A(p + f£, t?. 


rc) 


2 2 2 

r c — p — 




A{p + i£,q, 


rc) 



(47r)2 

1 



[0.133434] , 
7r)2 UttT/ 



67r2 (4 

[5.53165 X 10"^ 



+ 7.0292 



(D.2.85) 
, (D.2.86) 
(D.2.87) 



(471)2 

Adding these integrals to the subtracted integrals in ( D.2.60| , (D.2.61 1 and (D.2.631, we 



(D.2.20 


l,( 


D.2.21 


1 and ( 


D.2.24 



To evaluate the subtraction in the integrals ( D.2.62| and ( D.2.64| , we use the identity 
q'^ = [r^ + q^ — p^ — 2p ■ q) /2. The integral with q^ — p^ irv the numerator is purely 
imaginary. Thus the real part of the integral can be expressed as 



s-^n^(s) q^ 
pq 



Re( 



2 2 2 2 

r c — p — q 
A(p + k, q,rc) 
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2n2 



pq p^q 



A{p + ie,q,rc) 



(D.2.88) 



The evaluation of the first term in Eq. ( |D.2.88| l follows the same procedure as for ( D.2.20| , 
but just with Tip now instead of ng. The result reads 



pq p^q^ 



0{e). 



(D.2.89) 



A(p + i£,q,rc) 

It remains only to evaluate the integral in Eq. ( D.2.88| with p • q in the numerator. We 
begin by using the identity 



2 2 2 2 

^2 p • q y — — q"^ 



A{p + ie,q,rc) 



P + r / 2. / p q 

{p^-q^ + ie)^^ r2 



pqc" 



+ is ± q)^ \{p + z£ ± q)^ — ^^c^ 



(D.2.90) 



In the first term on the right side, the average over c is a simple multiplicative factor: 
(c^)c = 1/ (3 — 2e). The average over x gives hypergeometric functions of argument /6: 



p q 



\-e,\ 
3-2e 



2- e,l 

3- 2e 



(D.2.91) 



The integral over j6 gives hypergeometric functions of argument 1: 



,-1/2 p'^ + q^ /p-q 



1 f2)_ 



2e 



{l)-2e 
1 (3) -2. 



12 



2 - 2e, 1 - 1 
I -2e,3-2e 

Jl-e,l 



.l)-2e 

Expanding in powers of e, we obtain 



I -2e 



/2-2e,2-e,l 
I -2e,3-2e 

'2-e,l 
I -2e 



(D.2.92) 



JO 



-1/2 p^+ii /p-q 



(p2 _ ^2)2 \ ^2 16 



[1 - 1.02148 e] . (D.2.93) 



In the second term of ( D.2.90| l, the average over c is given by ( D.2.72[ |. In the -\-q 
term, the average over x = p ■ q is 



1 7'l 



1/ 



{p + qf 



5-2e 
4e 



J 2-eX\ 
\^3 -2e,l + e 
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5 (l),(l)-2.(3)-2.(D-. 

4e (l)-e(3)-3e 



1 -2e, 



5 _ 
2 

3-3e 



l-2g / 2 - 2e, I - 
1 - e i 3 - 3e 



(D.2.94) 



Integrating over we obtain hypergeometric functions of argument 1: 



Jo 



-1/2 



p qc" 



(2)- 



2e 



(p + (?)2-r2cVc,:, 

Y 2-2e,2-e,l,f 
\ |-2e,3-2e,l + e 



2 - 2e, 1 - e, 1, f 
; - 2e, 3 - 2e, 1 + e 



(l),(l)_2.(3)-2.(i)-. 



6e{2-3e) (l)_,(f) 



2J-3e 



2 - 3e, 1 - le, 



3e, 3 — 3e 



l-2e ^ / 2-3e,2-2e,§ - 
1-e I |-3e,3-3e 



(D.2.95) 



Expanding in powers of e, we obtain 



-1/2 



pqc" 



(p + — r^c^ 



c,x 



18 



(1 - 0.0728428 e) .(D.2.96) 



In the — term in the integral of the second term of ( D.2.90| l, we integrate over j6 be- 
fore averaging over x. The integral over /S can be expressed in terms of hypergeometric 
functions of type 2^1: 



Jo 



1/2 



4pq 

(p-(?)2 \{p + ie-q) 



2 _ ^2^2 



2(3-2e)e (i)_2. 



1-y 



(1) 



f^(l-2y) (l-y)-3/2F 



4(3-2e)e (-i), 
1 



2e,l 



1-y 



2 - 3e, : 



(D.2.97) 



The phase in the last term is e^^"^ for the f{s-^-, S-, r) term of (D.2.66 v which comes from 
the p > cj region of the integral, and e'"^ for the /(s-, s+, r) term, which comes from the 
p < q region. The average over x = p ■ q can be expressed in terms of hypergeometric 
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functions of type 3F2 evaluated at 1: 



Jo 



-1/2 



1 (2)-2. 

4(3-2e)e {\)_,, 



4p ■ q 



1 -e,2-2e,l 
3 - 2e, 1 + e 



1 (l),(3)_2.(-i)-. 



(3-2e)e (!)_,(- 1),(3)_2, 



- - — £ - —2e 1 
2 2 ' 

2 26", 2 ~1~ ^ 



1 + 



c,x 

2-e,2-2e,l 
3 - 2e, 1 + e 



e, ^ - 2e, 1 



l + 2e 



2(l-e) \3-26,-l+e 



16(2 -3e)e 



(l)e(2)-2.(2)-2.(i)-. 



(l)-.(3)-3. 



1 -e,2-3e, : 



3 — 3e, 3 — 3e 



1-g ^ /2-e,2-3e,| - 
l-2e I 3-3e,3-3e 



(D.2.98) 



The expansion of the real part of the integral in powers of e is 
Jo "^"^ XIP-'?) 



2 (p + fe — q)^ — r^c^ 



18 



(1 - 0.796858 e) . (D.2.99) 



Inserting p.2.93| , p.2.96| , and p.2.99| into the thermal integral of p.2.90| , we obtain 



r ^!|^P_qj, /,2 



2 2 2 2 

r^c — — q"^ 
A(p + ie,q,rc) 



^M^lp^j^ /^2 



pq P^'?^ ''^ 



2 2 2 2 

r c — p — 



(47r)2 V47rT/ 247r2 



+ 13.52098 



12 



log 2 



A(p + z£, ij,rc) (47r)2 6n^ 
Inserting these along with ( D.2.85| and ( D.2.89| into ( D.2.88| ^ we obtain 

1 



(D.2.100) 
(D.2.101) 



22 2 2 
r c — p — q 



— fV^Re ( c^— ^ 

pq p^r^ \ A(p + Z£,(j,rc 



(47r)2 l47rT/ 247r2 



15.302796 



pq P^''^ 



Re( 



22 2 2 
r c — p — q 



T2 1 - tt^ 



A{p + ie,q,rc) (An)^ 6n^ 



12 



log 2 



(D.2.102) 
(D.2.103) 



Adding this integral to the subtracted integral in (D.2.621 and ( D.2.64[ |, we obtain 
the final result in ( |D.2.22[ | and ( D.2.25[ |. The subtracted integral appearing in (D.2.65 1 



vanishes due to antisymmetry of the integrand. Thus the final result (D.2.26 1 is given 
by p.2.65|). 
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To evaluate the weighted averages over c of the thermal integrals in (D.2.271- 
( D.2.36[ l, we first isolate the divergent parts, which arise from the region q ^ 0. The 
integrals p.2.31| l and ( |D.2.32[ ) can be computed directly in three dimensions without 
any isolation of divergence, as described above. For the integrals ( D.2.27| and ( D.2.28| l, 
a single subtraction of the thermal distribution Mb suffices to remove the divergences: 



nsiq) = (nB{q) --)+-. (D.2.104) 
For the rest, a second subtraction is also needed to remove the divergences: 

nB{q) = (nsiq) - l + D+ l-^- (D.2.105) 



In the integral (D.2.301 and (D.2.34i, it is convenient to first use the identity r 



+ 2p • q/c + (j^/c^ to expand them into three integrals, two of which are (D.2.27 1 and 



(D.2.291, and ( D.2.31| and ( D.2.33| l, respectively. In the third integrals, the subtraction 
( D.2.105| is needed to remove the divergences. 

For the convergent terms, the HTL average over c and the angular average over 
X = p ■ q can be calculated in three dimensions: 



Re ( c 



-1 'c 



A{p + ie,q,rc 



c,x 



1 , 2p 



v + ^ log P±i - log til 



Re ( c 



' A{p + ie,q,rc 



Apq \2p + q p '^P ~ 1 

1 q^{Ap^ + 3q^) 2p 



(D.2.106) 



6(4p2-(j2) 3(4p2-(j2)3 



log 



{p + q){Ap^ + 2pq + q^ 
llpq{lp + qY 



log 



p + q {p - q){Ap'^ -2pq + q^ 



P 



llpq{lp-qf 



log 



P 

(D.2.107) 



Re ( p ■ q 



A(p + ^■£,^J,rc 



q{llp^-q^) Ap 



c,x 



6pq 6p(4p2 — (^2^2 



log 



{p + q){2p^-2pq-q^) p + q {p - q){2p^ +2pq - q^) \p - q\ 



12p^q{2p + qy 



2 2 2 2 2 

Re( ^'-P rj-p'-q^ 
q^ A{p + ie,q,rc) 



log 



,-1 



4p l2p'^q{2p-qY 
log 2' 



log 



4p 
(D.2.108) 



4pij2 



(J log 



r , c,x 

p + R 
\p-(^\ 



plog 



I 2 2 1 

\p — q \ 



(D.2.109) 
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The remaining two-dimensional integral over p and q can be evaluated numerically: 
J^<l p \ q q^J \ A{p + ie,q,rc) /, (4/1)2 L J 

(D.2.11( 
(D.2.11] 

Jpd p \ q q^ 2q) q^ \ A{p + ie,q,rc) (4/1)2 L J' 

(D.2.11: 

/ ^f^-^ + l)P^Re/ ;;-^.^-^\ \ = ^[-3.729 x 10-3] 
Jpq p \ q q^ 2qJ q^ \A{p + ie,q,rc)/, (4/r)2 L J 

(D.2.11; 

r ^fMi)_T IW^ 7^,,, rg-p2-,2 . ^ r2^[i.4S2xlO-], 
ipq p V 'J 2qJ q^ \ A{p + i8,q,rc) /, (4/r)2 L J 

(D.2.11^ 

/ ^(^-^ + i)P^Re(.-#:P^\ =^[-2.832 X 10-] 



(D.2.115) 

pq P q '^^\ q^ A{p + ie,q,rc) q^ q^ 

— 2 [4-134 X 10-2] , (D.2.116) 



(471) 



Mp)nFiq) ^^^/ rt-p^ r^-p^-q^ ^_i^2, _ 1^-1+2^ + 1^^2^ 
pq P \ ^2 A{p + ie,q,rc) q^ q^ 



j2 



(4/r)2 



2.530 X 10"^ 



(D.2.117) 



The integrals involving the terms subtracted from nB{q) in (D.2.1041 and (D.2.105 1 are 
divergent, so the HTL average over c and the angular average over x = p ■ q must be 
calculated in 3 — 2e dimensions. The first step in the calculation of the subtracted terms 
is to replace the average over c of the integral over q by an average over c and x: 

^ / f(A ^c-p^-q^ \ _ ^_TXn-l JLi^^2e(l)^/.7,,2Ne^7„Nl-«-2e 



X ^/(c) c3-"-2^(l - c2)"-2+2^ £(x ^ c - k)^-"-^'^ . (D.2.118) 
The integral over p can now be evaluated easily using either ( B.4.11| l and ( B.4.25| or 

/ nB{p)p-'-'^ = ^^^(l-46)(.V')'^T'-'% (D-2.n9) 

133 



APPENDIX D. THREE-DIMENSIONAL THERMAL INTEGRALS 



Jn Jn 



2-2e 



(D.2.120) 



It remains only to calculate the averages over c and x. The averages over x give 2F1 
hypergeometric functions with argument [(1 ^ c)/2 — ie]^^: 



{xTc-ie)-"-^'l = (ITC) 



-n-2er: l-e,n + 2e 
2-2e 



{x{x^ c — ie) 



-n-le\ 



:(lTc) 



l-e,n + 2e 
3-2e 



[(lTc)/2 j , 

(D.2.121) 

[{lTc)/2-ie]-'] 



2-£,n + 2e 
3-2e 



[(lTc)/2-fe]-i 



(D.2.122) 



Using a transformation formula, the arguments can be changed to (1 ^ c)/2 — ie. If the 
expressions ( D.2.12l| and ( D.2.122| | are averaged over c with a weight that is an even 
function of c, the + and — terms combine to give 3F2 hypergeometric functions with 
argument 1. For example. 



-l-2e 



c,x 



^-,7r. (l)3.(l)-2. j, l-2£,l-£,e 
e I (1)2.(2)-. I 2-£,l-3£ 



1 (2)-2.(l).(i)-. 

3e (l)-.(l 

, J2ne (1)-3.(1). ^ A + 1 + 2e,4e 

(l)_4e(2)2. I 2 + 2(r,l+3e 



(D.2.123) 



Upon expanding the hypergeometric functions in powers of e and taking the real parts, 
we obtain 



Re((l-c2)2-^(xTc 



ze 



-l-2e 



[-6 + 2(1- log 2)e2] , 



C,.T 



Re ^0^(1 - c^)^^ ^(x T c - ie) 
Re({l-c^)^+^'J^{xTc-iey 



-l-2e 



c,x 



3-2e 



c,x 



Re/x{l-c^)^+^'J^{xTc-iey 



2-2e 



c,x 



n 



n 



+ ^(2 -3 log 2)^2 



n 



|e + ^(l-61og2)e2 



(D.2.124) 
, (D.2.125) 
(D.2.126) 
.(D.2.127) 



If the expressions ( |D.2.12l| | and ( D.2.122[ | are averaged over c with a weight that 
is an odd function of c, they reduce to integrals of 2^1 hypergeometric functions with 
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argument y. For example. 



c{l-c^)^+^'J^{xTc-ie) 



-2-2e 



(2)-2.(i)-. J^2.-'-(l) 



(l)-.(l)- 
8 



c,x 



(2)2. Jo 



dyy-'^{l-yY+^\l-2y\F 



y 



,2!7Te(l)-3e 



3(1 +3e) (1)_4. 



j\yy^+^{l-yY+^\l-ly\F(^ 



-3e 
2 + 2e,l+4e 



2 + 36 



(D.2.128) 



The expansions of the integrals of the hypergeometric functions in powers of e are 
given in (F.0.23 l-( F.0.24[ |. The resulting expansions for the real parts of the averages 
over c and x are 



Re(c(l-c2)i+2^^(x^c 



le 



-2-2e 



c,x 



Re ( xc{l - c^f T c - ^^y^' 



le 



_^^14(l-log2)^^ 
2(1 -log 2) 



(D.2.129) 



c,x 



^ + ^log2-^log22+^)e. (D.2.130) 



Multiplying each of these expansions by the appropriate factors from the integral over 
q in p.2.118[ ) and the integral over p in p.2.119| l and ( |D.2.120| , or ( |B.4.11| and ( |B.4.25| l, 
we obtain 



222 
n — p — q 



"b(p) 1 Re /c-^+^^^ 

pq p \ h{p + k,q,rc, 



l7r)2 V 



(47r)2 V47rT 



4e 



A{p + i£,q,rc) 
1 

24 



e 
T 



■2 + 41og(27r) 



(D.2.131) 



4e 



nsip) P^ 

pq P 

Mp) p q 

pq P 



Re(ci+2^^ 



A{p + ie,q,rc 



Re ( c 



,2e- 'c 



^ (47r)2 V47rT 
^ + ^+41og(27r) 

r /_ 1 

~ (47r)2 V 12 



(D.2.132) 
(D.2.133) 



(471)2 V47rT 



24 



1 11 



nsip) P^ 

pq P 



Re/ci+2-4f 



e 3 

T2 / N4e 



+ 41og(27r) 



XD.2.134) 



A{p + ie,q,rcj 
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1 

24 



(D.2.135) 



pq 



2 2 2 

— p — q 



1^ \ 



4e 



1 

18 



(l-log2) 



(47r)2 V47rT 



14 , 
3 ^ 



12 



(D.2.136) 



pq p 



2 2 2 

— p — q 



--3 \ A{p + i£,q,r, 



(4 



7r)2 V47rT/ 



4e 



1 

48 



12 4,^ 

loe2 + 4^-^ 

6 3 3'°8^+^a-i) 



pq p r 



2 2 2 

— p — 



pq P 

X 



1 

36 

^f(p) P q 

pq P 



A{p + k,q,rc 



2 2 2 
p — 



A(p + /e,ij,rc 



(D.2.137) 
(D.2.138) 



4e 



(4 



1 14 

(l-log2)(- + --41og2 + 4:^^_jj 



12 



Re ( c 



,2e '^c 



(D.2.139) 
(D.2.140) 



A{p + i£,q,rc) (47r)2 

Adding these integrals to the subtracted integrals in (D.2.110 l-( D.2.112] l, we obtain 
the final results in p.2.27[ )-( [a229| . Combining ( |D.2.115| l with p.2.134| and ( |D.2.136| l, 
we obtain 



/" nB{p)nB{q) p q ^^^^2e 
Jpq pq q^ 



2 2 2 

n-p - r 



A{p + ie,q,rc 
ji \4e5-21og2 



72 



+ 11.6689 



. (D.2.141) 



(47r)2 V47rT. 

The final integral p.2.30| is obtained from p.2.27| , p.2.29| , and p.2.141| by using the 
identity = p^ + 2p • q/c + q'^/c^. 

Adding Eqs. ( |D.2.137[ ) and p.2.138| to the subtracted integral ( |D.2.114| we obtain 
the final result in Eq. ( |D.2.33| >. Combining ( |D.2.115| | with ( |D.2.139| | and p.l.UO} , we 
obtain 



^F(p)nB(^) p-q, 

pq P'? 



Re(c2^^ 



2 2 2 

rf — p — q 



(47r)2 V47rr 



Aip + ie,q,rc) 







- - 15.2566 




e 



(D.2.142) 



The integral ( |D.2.34[ | is obtained from ( |D.2.31[ ), ( |D.2.33| | and ( |D.2.142[ |. Finally con- 
sider ( D.2.35| l and (D.2.361. In order to evaluate them we need two subtractions for 
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each integral 



nf(p)nf(p) 1_ 

pq P'? 



(4 

1 

"12 



7r)2 V47rT; 



4e 



-+2 + 21og2 + 27 + 2|^ 



nF{p)nF{q) 1 ^^-i+2e^^ 



pq 



1 

~24 



7r)2 V47rTy 



4e 



(471) 



nB{p)nF{q) ]_ ,2e^ 

pq P'? 



^+ (2 + 27 + 41og2 + 2? ^ 

j2 



a-i: 



53.1065 



(471) 



V47rT/ 



4e 



■2 + 41og2 + 27 + 2 



pq 



nB{p)nF{q) J_ 
P'? 



(4 



— m 
7r)2 V47rTy 



4e 



1 

12 



+ 2 + 27 + 61og2 + 2 



a-i) 



+ 69.7096 



(D.2.143) 



(D.2.144) 



(D.2.145) 



(D.2.146) 



The subtractions can be evaluated directly in three dimensions and the results are 
given in Eqs. (D.2.116|-(D.2.117| The integrals ( D.2.35[ l and ( D.2.36[ l are then given 
by the by the sum of the difference terms (D.2.1161 and ( D.2.117| and the subtraction 
terms p.2.143| l-( |D-2.146| |. 
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Four-Dimensional Integrals 



In the sum-integral formula (B.3.151, the second term on the right side involves an 
integral over four-dimensional Euclidean momenta. The integrands are functions of 
the integration variable Q and R = — (P + Q). The simplest integrals to evaluate are 
those whose integrands are independent of Pq: 





1 




(471) 


1 


1 




(471) 




1 




(471) 



1 



-F4-21og2 



- + l-21og2 



-2) [l + (-2-21og2)e] 



(E.0.1) 
(E.0.2) 
(E.0.3) 



Another simple integral that is needed depends only on = Pq -\- jp-. 



(E.0.4) 



where (fl)^ is Pochhammer's symbol which is defined in (F.0.3 1. We need the following 
weighted averages over c of this function evaluated at P = (— /p, p/c): 



-1+2? 



P^(-ip,p/c) 



P^{—ip,p/ c) 



(471) 



(471, 



^+21og^2 ' 



21og2 



4 

(E.0.5) 
(E.0.6) 



The remaining integrals are functions of Pq that must be analytically continued to the 
point Po = —ip + e. Several of these integrals are straightforward to evaluate: 



= 0, 



(E.0.7) 



Po=-ip 
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Q Q^r^R^ 



Po=-ip 



{Any 



1 l-21og2 

^ + — 



+ 10-21og2 + 21og^2 



7n^ 



Po=-ip 



(471) 



2 -21og2 



(E.0.8) 
(E.0.9) 



We also need two weighted average over c of the integral in ( E.0.7[ | evaluated at 
P = {—ip,p/c). The integral itself is 



qQ^R^ 



(47r 



1 ^^,^2)y-2.(l).l(l)-.(l)- 



3-2e 



e 4 (2)_2. 



(E.0.10) 



The weighted average is 



,l+2e 



qQ^R^ 



1 ,2e^2-2e 1 



P^{-ip,p/c) I 



(47r 



48 



1 2(10 + 3 log 2) 



3e 



4 40, ^ ^, 2. 37r2 
- + ylog2 + 2V2+ — 



. (E.0.11) 



The most difficult four-dimensional integrals to evaluate involve an HTL average 
of an integral with denominator Rq + r'^c^: 



Re 



Re 



Re 



1 /c2(l-c2 



QQ^Rg + J-'cV, 3(47r) 



Re 



Re 



2 2 



QQ2r2(R2 + r2c2) 



(471) 



2 -21og2 



8-41og2 + 41og^2 



71 



1 1 



(47r: 



1 20 

^ + y- 

5-61og2 
3e 



61og2 



|-21og2 + 41og22-^ 



1 1 



2e„-2-2e 



4 (47r)2 
1 



14 2 



(4/r: 



e 3 3 

13 -16 log 2 
12i 



log 2 



29 19, ^ 8, 2. 4 2 

log2+ -log'^2 TT^ 

9 18 * 3 9 



(471; 



(E.0.12) 
(E.0.13) 



(E.0.14) 
(E.0.15) 



(E.0.16) 
(E.0.17) 
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The analytic continuation to Pq = — ip + c is implied in these integrals and in all the 
four-dimensional integrals in the remainder of this subsection. 

We proceed to describe the evaluation of the integrals ( E.0.12[ | and ( E.0.14[ |. The 
integral over Qo can be evaluated by introducing a Feynman parameter to combine 
and Rq + r-^c^ into a single denominator: 



QQHRi + rV) 



-3/2 



x + xc^y + 2(1 -x)r-p + {l- xYp^- is] , (E.0.18) 



where we have carried out the analytic continuation to Pq = + ^- Integrating over 
r and then over the Feynman parameter, we get a j¥\ hypergeometric function with 
argument 1 — c^: 



qQ2(K2 + ,2c2) (47,)2 

x(l-c2)-^F I 2 



e (2) -3. 

3 le,l-e 



2-3e 



(E.0.19) 



The subsequent weighted averages over c give hypergeometric functions with 
argument 1: 



j_ / c2 



c^(l-c^) 
Q \ R§ + r2c2 



1 

(47r)2 



(l).l.,..(2) -^(l)-2.a)-2. 

5^ 



e 3 (i)-2.(2)-3. 



xf 



' 1 - 2e, I - 2e, 1 - e 
^ - 2e, 2 - 3e 



(E.0.20) 



(471 



1 (g7^2^ep-2e Ae'^e (2)-e(l)-2e(2)-2e 



2-2e 



e 15 
2e, 1 - e 



2e, 2 - 3e 



(D-2.(2)-3. 

1) . (E.0.21) 



After expanding in powers of e, the real part is (E.0.141. 

The integral (E.0.15 1 has a factor of l/r^ in the integrand. After using (E.0.18 1, it is 



convenient to use a second Feynman parameter to combine (1 — x + xc^y with the 
other denominator before integrating over r: 

X j [{1-X + xc^y + 2y{l - x)r-p + i/(l - x)^p^ - ie] '^^^ . 

(E.0.22) 
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After integrating over r and then y, we obtain 2^1 hypergeometric functions with 
arguments x(l — c^). The integral over x gives a 2^1 hypergeometric function with 
argument 1 — c^: 



QQ^r^{Rl + r^c^) (4/1)2 



7 2\e -2-2e 



2(1 + 2e) 



3e 



{\)-le 

-{l-c^)-'F 



-3e 



(E.0.23) 



After averaging over c, we get a hypergeometric functions with argument 1: 



.(e>2)^p-2-2. 



(1), / 1 (-^)-.(l)- 



3-2e 



(D- 



2e 



1^,,, (-|)-.(l)-2.(2)-2. j. / 1 - 2e, i - 2e, -e 
2 (i)-2.(l)-3. V 5-2e,-3e 



(E.0.24) 



After expanding in powers of e, the real part is (E.0.15 1. 



To evaluate the integral (E.O.I6I, it is convenient to first express it as the sum of 
three integrals by expanding the factor of in the numerator as = + 2Tp ■ r + r^: 



t = / fr +2^:^ + 1 



2^2 ^ 



(E.0.25) 



To evaluate the integral with p • r in the numerator, we first combine the denominators 
using Feynman parameters as in ( E.0.22| . After integrating over r and then y, we obtain 
2F1 hypergeometric functions with arguments x(l — c^). The integral over x gives 2^1 
hypergeometric functions with arguments 1 — c^: 



pr 



(1)- 



3e- 



2^' t ^5^-2^ 
,7re (l)-2g(l)-e ^-^ _ ^^ye^ [ 2 ~ 2e, -£ 



(D- 



1 -3e 



l-c2 



(E.0.26) 



After averaging over c, we get a hypergeometric function with argument 1: 



(4/1): 



7 2\e -2e^ 



-(.V)^p 



2e2 



3-2e (I) 



2;-2e 
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+ 



l^rne (l)-e{-^)-le{l)-2e ^ / 1 " | - le, -6 

3 (|)-2e(l)-3. V 2-2e,l-3e 



(E.0.27) 



After expanding in powers of e, the real part is 

^2 



Re 



L 



p r 



1 



-l+log2 
3e 



20 14, ^ 2, 2^ 

\ log 2 log^ 2 H 

9 9^ 3^ 36 



. (E.0.28) 



Combining this with ( E.0.12[ | and ( E.0.14| , we obtain the integral ( E.0.16| |. 

To evaluate the integral ( E.0.17[ |, we first express the numerator as a sum of two 
integrals whose averages have been calculated: 



2 2 

q — p 



2p ■ r + 



QQ2r2(R2 + r2c2) 



(1)-.(1)-2. 1 
(l)-3. e 

(l)-.(l)-2. 



^ (i)-2. 



3e 



(2)- 



3e 



(E.0.29) 



The two hypergeometric functions are now combined into a single hypergeometric 
functions, which yields 



I — 
Jo 



QQ2r2(R2 + ,2,2) 



Averaging over c, yields 
2p ■ r + 



QQ^r^{Rl + r^c^) 



c,x 



(471) 



(i) 



(l)-.(2)-2. 



(2)- 



3e 



2-3e 



1-c^ 



(E.0.30) 



(471)^ 



(i) 



2)-2e 



-1+e' 



(1)- 



2e 



(E.0.31) 



Expansion in powers of e, yields Eq. ( E.0.17| . 
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Hypergeometric Functions 



The generalized hypergeometric function of type pF^ is an analytic function of one 
variable with p + q parameters. In our case, the parameters are functions of e, so the 
list of parameters sometimes gets lengthy and the standard notation for these functions 
becomes cumbersome. We therefore introduce a more concise notation: 



[ h N 



(F.0.1) 



The generalized hypergeometric function has a power series representation: 



/ OCi, 0t2, . . . , (Xp 
\ ^ ^ 



where {a)i, is Pochhammer's symbol: 



^ A («l)n(«2)n- • • («p)n „ 



r(fl + b) 



(F.0.2) 



r(«) ■ 

The power series converges for |z| < 1. For z = 1, it converges if Res > 0, where 



(F.0.3) 



p-i p 



(F.0.4) 



i=l i=l 

The hypergeometric function of type p+if^+i has an integral representation in terms of 
the hypergeometric function of type pF^: 



Jo y fii,...,pci 



tz = 



I{H + v) I j6i,...,j3^,^ + v 



(F.0.5) 
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If a hypergeometric function has an upper and lower parameter that are equal, both 
parameters can be deleted: 



0Ci,0i2, ...,Oip,V 



^ I Oil, 0i2, . . . , dp 

\ h N 



(E0.6) 



The simplest hypergeometric function is the one of type iFq- It can be expressed in 
an analytic form: 



iFo(«;;z) = (i-z) 



(E0.7) 



The next simplest hypergeometric functions are those of type 2F1. They satisfy trans- 
formation formulas that allow an 2F1 with argument z to be expressed in terms of an 
2F1 with argument z/ (z — 1) or as a sum of two 2Fi's with arguments 1 — z or 1/z 
or 1/(1 — z). The hypergeometric functions of type 2F1 with argument z = 1 can be 
evaluated analytically in terms of gamma functions: 



0ii,Ci2 



r(|6i)r(|6i-«i-^2) 
r(^i-ai)r(^i-a2) 



(E0.8) 



The hypergeometric function of type 3F2 with argument z = 1 can be expressed as a 3F2 
with argument z = 1 and different parameters | 



Cii,Ci2,0iz 



r(/3i)r(/32)r(s) 



r(ai+s)r(a2 + s)r(a3) 



|6i - aa, j62 - ^3/ s 

(X.I+ 5,0.2 + S 



(E0.9) 



where s = j6i + /52 — ai — ^2 — *3- If all the parameters of a 3F2 are integers and half-odd 
integers, this identity can be used to obtain equal numbers of half-odd integers among 
the upper and lower parameters. If the parameters of a 3F2 reduce to integers and 
half-odd integers in the limit e — > 0, the use of this identity simplifies the expansion of 
the hypergeometric functions in powers of e . 

The most important integration formulas involving 2F1 hypergeometric functions 
is ( F.0.5| with p = 2 and q = 1. Another useful integration formula is 



t 



i-t I 

T{al + ^)Tia2 + ^)Ti^l)T{-^) 

r(ai)r(a2)r(/3i + ^) 



T{li)T{v) ^ f CLi,Ci2,V 



(-z)''F 



OL1 + H,0L2 + n,v + n 
^1 + jM, 1 + f/ 



(EO.IO) 



This is derived by first inserting the integral representation for 2F1 in ( F.0.5[ l with 
integration variable t' and then evaluating the integral over t to get a 2F1 with argument 
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1 + t'z. After using a transformation formula to change the argument to —t'z, the 



remaining integrals over t' are evaluated using (F.0.5 1 to get sFa's with arguments — z. 

For the calculation of two-loop thermal integrals involving HTL averages, we 
require the expansion in powers of e for hypergeometric functions of type pfp-i with 
argument 1 and parameters that are linear in e. If the power series representation 
( F.0.2[ | of the hypergeometric function is convergent at z = 1 for e = 0, this can be 
accomplished simply by expanding the summand in powers of e and then evaluating 
the sums. If the power series is divergent, we must make subtractions on the sum 
before expanding in powers of e. The convergence properties of the power series at 
z = 1 is determined by the variable s defined in ( F.0.4[ |. If s > 0, the power series 
converges. If s — )■ in the limit e — ?► 0, only one subtraction is necessary to make the 
sum convergent: 



Ui h-i 




r(/3i)---r(/3p_i: 



r(ai)r(a2) ■ ■ ■r(ap) 

( *1 jn 



v-^' (n + 1)-'-^ 

(/3i)„ ■ ■ ■ (/3g)„n! T{a^)T{cc2)---T{ocpy ^ > 

(RO.ll) 



If s ^ 
gent: 



- 1 in the limit e — > 0, two subtractions are necessary to make the sum conver- 



Ul 



r(/3i)---r(/3p_i) 
r(ai)r(a2) • • •r(ap) 



[as + l)+tas + 2)] 



r(/5i)---r(/3p_i^ 



h V (i6i)„---(/5,)„n! r(«i)r(«2) • • ■ r(ap) 



{n + l)-'-^ + t{n + l) 



-s-2 



(R0.12) 



where t is given by 



^ ^ ^ K - 2) g (^,-l)(/3,-2) 
i=l ^ !=1 ^ 



(R0.13) 



The expansion of a pFp_i hypergeometric function in powers of e is particularly 
simple if in the limit e — > all its parameters are integers or half-odd-integers, with 
equal numbers of half-odd-integers among the upper and lower parameters. If the 
power series representation for such a hypergeometric function is expanded in powers 
of e, the terms in the summand will be rational functions of n, possibly multiplied by 
factors of the polylogarithm function xp{n + a) or its derivatives. The terms in the sums 
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can often be simplified by using the obvious identity 



k-l 



J^lf(n)-f{n+k)] = J^fii) 

n=0 !=0 



(E0.14) 



The sums over n of rational functions of n can be evaluated by applying the partial 
fraction decomposition and then using identities such as 



E 

n=0 



n + a n + h 



ip{h)-xp{a), 
rp'ia) . 



(E0.15) 
(E0.16) 



The sums of polygamma functions of n + 1 orn + ^ divided by n + 1 orn + ^ can be 
evaluated using 



E 

oo 

E 

oo 

E 

oo 

E 

)i=0 



xp{n + l) log(n + l) 



n + 1 


n + 1 


rp{n + l) 


log(n + 1) 


n + \ 


n + 1 




log(n + l) 


n + 1 


n + 1 


tp{n + l) 


log(n + l) 


n + \ 


n + 1 



i 2 ^ 

22 '^^ ' 



-^(7 + 21og2)2 + ^-7i, 



-7^ -41og2 + 21og^2 ^ 



12 



-^(7 + 21og2)2-^-7i, 



(E0.17) 
(E0.18) 
71 , (F.0.19) 
(E0.20) 



where 71 is Stieltje's first gamma constant defined in ( B.1.14| . The sums of polygamma 
functions of n + 1 or n + j can be evaluated using 



i:^(^^(n + l)-log(n + l) + ^^ 



^ + ^7-^log(27r), (F.0.21) 



00 / 1 \ 1 1 

(^t/'(n + l)-log(n + l) + — J = -7-log2--log(27r). (F.0.22) 

We also need the expansions in e of some integrals of 2^1 hypergeometric functions 
of y that have a factor of 1 1 — 2i/ 1 . For example, the following two integrals are needed 
to obtain ( |D.2.129| >: 







j\yy'+%l-y)'+^\l-2y\F(^ 



2 + 2e,l + e 
2 + 3e 



^+0 + ^og2)^^ (F-0-23) 
|+(^ + ^log2)e. (F.0.24) 



146 



APPENDIX F. HYPERGEOMETRIC FUNCTIONS 



These integrals can be evaluated by expressing them in the form 



2jjdyf-Hl-yr-Hl-2y)Fiy^' 



y 



(R0.25) 



The evaluation of the first integral on the right side gives 3F2 hypergeometric functions 
with argument 1. The integrals from to ^ can be evaluated by expanding the power 
series representation ( F.0.2| of the hypergeometric function in powers of e. The resulting 
series can be summed analytically and then the integral over y can be evaluated. 
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